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quantitative  assessment  of 

ENVIRONMENTAL  IMPACTS  IN 
THE  AQUATIC  ENVIRONMENT 


1  INTRODUCTION 


Background 

Ouring  the  first  decade  of  legislated  environmental  assessment  of  govern¬ 
mental  activities,  there  have  been  changes  in  approaches  and  attitudes  toward 
the  National  Environmental  Policy  Act  (NEPA),  as  well  as  implementation  of  new 
legislation  and  administrative  guidelines  requiring  new  forms  of  environmental 
planning.  New  terminology,  such  as  "scoping,"  "fate  and  effect,"  "hazard 
evaluation,"  etc.,  has  appeared,  and  new  methodologies  have  been  developed  for 
environmental  planning  and  management.  Unfortunately,  legislative  acts  and 
rhetoric  tend  to  grow  faster  than  do  the  tools  to  carry  them  out.  As  a 
result,  quantifying  and  measuring  the  significance  of  environmental  impacts 
still  remains  an  plusive  target  during  planning,  and  flexible  and  efficient 
tools  are  needed  to  predict,  evaluate,  and  mitigate  these  impacts.  To  meet 
this  need,  the  U.S.  Army  Construction  Engineering  Research  Laboratory  (CERL) 
has  developed  a  prototype  package  of  computerized  impact  evaluation  procedures 
called  the  Rational  Impact  Assessment  System  ( R IAS)  which  use  site-specific 
quantification  routines  to  answer  such  questions  as,  "How  bad  will  the  impact 
be?"  or  "Will  the  impact  be  significant?" 


Objective 


The  objective  of  this  report  is  to  document  the  development  of  RIAS  as  an 
impact  evaluation  tool. 


Approach 

Because  the  mechanisms  of  environmental  impacts  can  be  extremely  complex 
and  varied,  it  is  difficult  to  construct  one  comprehensive  simulation  tool  for 
predicting  them.  For  this  reason,  the  computer  software  developed  to  support 
RIAS  consists  of  a  series  of  independent  modules  which  can  be  used  either  as 
separate  programs  or  together  as  subroutines  within  a  larger  control  program. 
The  number  of  modules  used  depends  entirely  on  the  decision-maker's  needs  and 
the  types  of  impacts  identified  through  initial  scoping. 

Data  collection,  processing,  and  impact  simulation  carried  out  in  RIAS 
proceed  through  the  use  of  three  general  procedures:  Filter  Questions 
(FQUES),  Water  Quality  Simulation  (SIMWQ),  and  Rational  Threshold  Value  Test 
(RTVTEST)  (see  Figure  1).  Appendix  C  provides  program  listings  for  these 
three  systems.  FQUES  collects  and  organizes  relevant  environmental  setting 
and  project  information  data  through  a  computerized  format  of  filter  ques¬ 
tions.  SIMWQ  simulates  primary  impacts  on  the  physical /chemical  attributes  of 
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the  aquatic  receiving  system  over  both  temporal  and  spatial  dimensions. 

RTVTEST  Is  a  set  of  rational  threshold  values  (RTV)  models  which  predict  the 
significance  of  the  primary  and  secondary  impacts  listed  by  SIMWQ.  This 
series  of  analyses  simulates  the  Impact  chain  of  events  and  provides  a  uniform 
method  for  environmental  data  handling. 

Chapters  2  and  3  document  the  development  of  RIAS  as  an  impact  signifi¬ 
cance  evaluation  tool,  describing  these  three  supporting  modules  and  providing 
user  information  for  application  of  RIAS  to  the  quantitative  assessment  of 
impacts  on  aquatic  ecosystems.  Four  example  applications  of  RIAS  (Chapter  4) 
and  user  information  for  data  collection  and  organization  are  presented 
(Chapter  5). 


Mode  of  Technology  Transfer 

The  information  in  this  report  will  be  Issued  as  a  DA  Pamphlet  in  the  200 
series  and  as  the  module  called  RIAS  in  the  remote  terminal  ADP  system  enti¬ 
tled  Environmental  Technical  Information  System  (ETIS). 


from  scoping 
procedures 


Potential 

Impacts 
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Figure  1.  General  design  and  information  flow  for  the  RIAS. 


2  FILTER  QUESTIONS 

FQUES  has  been  developed  to  execute  a  program  which  asks  a  series  of 
questions  about  a  specific  project's  environmental  setting  and  organizes  the 
information  obtained  into  a  data  file.  The  data  file  is  then  used  as  input 
for  simulations  and  other  evaluation  protocols. 

FQUES  can  be  used  to  set  up  a  new  data  file  or  to  revise  an  existing  data 
file.  Output  from  FQUES  stores  data  related  to  such  things  as  number  of 
stream  reaches,  number  of  conservative  and  nonconservative  water  quality 
attributes,  tributary  inputs,  point  source  discharges,  hydraulic  rating  param¬ 
eters,  boundary  conditions,  and  biological  parameters. 


3  WATER  QUALITY  SIMULATIONS  (SIMWQ) 


The  development  of  SIMWQ  has  been  restricted  to  one-dimensional,  non- 
dispersive,  steady-state,  plug  flow  models  for  receiving  streams.  Whenever 
possible,  rate  equations  have  been  limited  to  first-order  reaction  kinetics. 
This  results  in  linear  rate  equations  which  provide  analytical  solutions.  In 
addition,  it  avoids  the  necessity  of  using  complicated  numerical  solution 
techniques  for  sets  of  differential  equations  and  produces  a  much  more  usable 
model . 

SIMWQ  considers  two  types  of  sources  or  sinks  for  water  quality  consti¬ 
tuents.  First-order  decay  or  accumulation  terms,  similar  to  the  form  of  the 
Streeter-Phel ps  equation,  are  used  to  represent  most  biological  activity  in 
the  stream.  These  terms  can  also  be  used  for  distributed  sources  or  sinks 
which  affect  the  stream  equally  along  a  longitudinal  gradient.  Examples  of 
these  terms  are  benthic  oxygen  demand  or  nonpoint  source  runoff.  The  receiv¬ 
ing  watershed  is  represented  by  a  series  of  stream  reaches  within  which  all 
model  parameters  are  constant.  At  the  end  of  each  reach,  new  model  parameters 
are  calculated,  based  on  local  environmental  data  and  point  source  inputs. 
These  results  are  combined  with  upstream  values  based  on  conservation  of  mass 
and  assuming  complete  mix.  This  modeling  approach  is  not  new;  however,  it  is 
an  efficient,  flexible  system  for  tracing  changes  in  water  quality. 

One  major  addition  to  the  computer  program  which  is  the  basis  for  SIMWQ 
is  the  capability  to  handle  branched  watersheds.  This  required  coding  the 
boundary  conditions  for  reaches  rather  than  a  modification  of  the  analytical 
solutions. 

Figure  2  shows  the  general  conceptual  design  of  SIMWQ  and  the  variables 
it  can  analyze.  Tables  1  and  2  define  the  model's  variables  and  parameters, 
respectively.  Table  3  presents  the  rate  equations  which  represent  the  heart 
of  the  model.  Table  4  lists  the  analytical  solutions  derived  from  the  equa¬ 
tions  provided  in  Table  3.  These  solutions  assume  time-constant  model  parame¬ 
ters  and  apply  only  within  reaches.  Tributaries  and  any  applicable  point 
source  inputs  between  reaches  are  accounted  for  by  Eq  1: 


QlCi  +  Q2C2 

Ql  +  Q2 


[Eq  1*] 


where  C  denotes  concentration,  Q  is  flow,  and  subscripts  1  and  2  refer  to  the 
different  flows  being  combined. 


Several  simplifying  assumptions  were  used  to  reduce  the  system  of  rate 
equations  to  the  desired  forms,  including: 


1.  Algae  concentrations  will  be  relatively  constant  for  a  given  reach 
and  a  given  season  of  the  year. 

2.  Higher  organisms,  such  as  fish  and  invertebrates,  do  not  signifi¬ 
cantly  affect  the  rate  of  concentration  change  for  any  of  the  attri¬ 
butes  considered. 


^Variables  for  all  equations  in  text  are  defined  in  Table  2. 


10 


3.  For  the  purpose  of  SIMWQ,  certain  water  quality  attributes  can  effec¬ 
tively  be  considered  conservative  substances  (e.g.,  TDS ,  hardness, 
pH,  total  alkalinity). 

These,  and  other  assumptions  involved  with  steady-state  and  nondispersl ve 
models,  must  be  reevaluated  in  each  application  of  SIMWQ. 


Distributed 

Sources 


Point  Source 
ft  Tributary 


a  Sinks  lnPuts 


Figure  2.  Conceptual  organization  of  interactions  between  water 
quality  attributes  within  each  stream  reach  in  SIMWQ. 
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Table  1 


SIMWQ  Variable  Listing 


Symbol 

Definition 

Units 

ALK 

Total  alkalinity  as  CaC03 

mg/d 

C 

Total  dissolved  inorganic  carbon 

mg/d 

CDET 

Carbon  content  of  detritus  in  water  column 

mg/d 

C02 

Oissolved  carbon  dioxide 

mg/d 

COLI 

Col i form  bacteria  concentration 

#/100  md 

CON  Si 

ith  additional  conservative  water  quality 
constituent 

mg/d 

D 

Dissolved  oxygen  saturation  deficit 

mg/d 

HARD 

Water  hardness  as  CaC03 

mg/d 

L1 

Biochemical  oxygen  demand  (BOD5) 

mg/d 

NCONi 

1th  additional  nonconservative  water  quality 
constituent 

mg/d 

NDET 

Nitrogen  content  of  detritus 

mg/d 

NH3 

Total  ammonia  nitrogen 

mg/d 

N02 

Nitrite  nitrogen 

mg/d 

N03 

Nitrate  nitrogen 

mg/d 

02 

Dissolved  oxygen 

mg/d 

POET 

Phosphorus  content  of  detritus 

mg/d 

pH 

pH  of  surface  water  (-logio  CH+]) 

— 

P04 

Phosphate- phosphorus 

mg/d 

TW 

Temperature  of  water 

°C 

TSS 

Total  suspended  solids 

mg/d 

TOS 

Total  dissolved  solids 

'mg/d 
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Table  2 


SIMWQ  Parameter  Listing 

★ 

Symbol  Definition  Values  Sources 

1.  Distributed  Sources/Sinks 

SL  Scour/runoff  of  BOD^ 

SNrn(,.  Monpoint  sources  of  non- 

'  “  "  conservative  pol lutants 

SNH3  Surface  runoff  of  NH^ 

SN03  Surface  runoff  of  NO^ 

Sq2  Daily  mean  net  1°  production 

Sp04  Surface  runoff  of  P04 

SS0D  Sediment  oxygen  demand 

Sss  Scour/erosional  inputs  of 

suspended  solids 

ST  Natural  heat  inputs  from 

atmosphere 

2.  Reaction  Rates  (k[  =  gJ“2  •  K?°°C) 

kANH3  Algal  uptake  of  NHj  0.1  to  4.0 

kAN03  Algal  uptake  of  NOj  0.1  to  4.0  White  &  Dracup 

kAP04  Algal  uptake  of  P04  0.005  to  0.5  White  &  Dracup 

k^  Decomposition  of  BODg  0.01  to  2.5  HEC,  Zison 

k£s  Bottom  exchange  of  B0D5  0.0  to  2.0 

*See  References,  pp  32-33. 
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Table  2  (Cont'd) 


Symbol  Definition  Values 


Sources 


* 


kNC0Ni 

Biodegradation/ decay  of 

attribute 

ith  nonconservative 
pollutant 

dependent 

kT 

kndet 

Mineral izati on/dissolu¬ 

0.001  to  0.02 

HEC 

tion  of  nitrogen  portion 
of  detritus 

J 

nH3 

Microbial  conversion  of 

0.01  to  2.5 

Zison,  pp  188-197 

NH3 

HEC 

kN02 

Microbial  conversion  of 

0.020  to  0.5 

Miller  4  Jennings 

no2 

0  to  10.0 

White  4  Dracup  HEC 

kT 

*N03 

Microbial  conversion 

«  0.001 

Miller  4  Jennings 

of  N03 

kT 

kpdet 

Mineral ization/dissolu- 
tion  of  phosphorus 
portion  of  detritus 

0.001  to  0.02 

HEC 

(0.01) 

*P04 

Decay  of  phosphates  due 

0.005  to  0.5 

White  4  Dracup 

to  microbial  uptake/ 
conversion 

* 

Reaeration  of  dissolved 
oxygen 

kjj°  =  a  Ub  HC 

Covar 

kT 

*ss 

Decay  of  suspended 
solids  (mineralization/ 
biotic  breakdown) 

0.001  to  0.02 

HEC 

kT 

Ksss 

Settling  of  suspended 
sol  ids 

0  to  2.0 

HEC 

*See  References  pp  32,  33. 
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Table  2  (Cont'd) 


Symbol 

Defini tion 

Values 

Heat  exchange  with 

(See  Temp.  Model 

atmosphere 

in 

the  Computer 
Li  sting) 

kT 

ktss 

Overall  decay  of 
suspended  solids 

^kTSS  =  kSS  +  kSSS/dePth) 

kT 

KSS 

+  ksss/dePth 

★ 

Sources 


*See  References  pp  32,  33. 
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Table  3 


SIMWQ  Rate  Equations* 


Rate  Equations  (within  reaches,  excluding  point  sources) 


1.  Temperature 


2.  BOD 


£'ST-kT  T‘ 


£'-<kI*kIs>  Lt  +  $L 


3,  Suspended  Solids 


-  -^ss  •  ISSt  ♦  sss 


dt 


..  T  *$SS  ,  TCC  .  c 

(kSS  flepflT  TSSt  Sss 


4.  Phosphate-Phosphorus 


TT  *  'kP04  •  ro«t  *  kJoET  •  P0ETt  *  SW>4  -  kjp04  •  A 


5.  Ammonia-Nitrogen 


^  *  -“JhS  •  NH3t  *  kN0ET  •  "DETt  *  S"«  *  kANH3  •  A 


6.  Nitrite-Nitrogen 


=  -kJQ2  •  N02t  +  kJH3  *  NH3t 


*  Equation  variables  are  defined  In  Table  2. 
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Table  3  (Cont’d) 


7.  Nitrate-Nitrogen 


dN03 

dt 


*  ^02  •  N02t  -  k 


N03 


N03  +  S(403  +  k 


AN03 


8.  Dissolved  Oxygen  Deficit 


dD  _  .T  ,  .  T  .  NH3  A  .T  ,  N02  .  T  _ 

—  -  -kL  «  Lt  kNH3  •  Lt  kHQ2  Lt  -  kR  Ot  -  Sq2 


9.  Conservative  Constituents 


dCOHi  =  Q 
dt 

10.  Nonconservative  Constituents 

dNCONI  _  .  T  urnu.  ,  c 

— dt —  "  “kNCONi  '  NC0Nl  SNC0N1 


sSOD 
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Analytical  Solutions  to  Table  3* 

(These  equations  apply  only  within  stream  reaches 
which  have  constant  model  parameters.) 


1.  Biochemical  Oxygen  Demand  (L) 

L  T  T  L 

Lt  =  (L0  -  ITj)  expMk^  +  k^s)t)  +  Fj 

-L  SL 

where  J<j  =  - 

^kL  +  kLS^ 

2.  Total  Dissolved  Solids  (TSS) 

TSSt  -  (TSS0  -  exp(-(kjs  +  ^_)t)  +  ™ 


_TSS 
where  Kj 


_ $SS _ 

*kSS  +  kSS  +  kSSSdepth) 


3.  Suspended  Solids  Portioning 


CDET  t  =  PC(TSSt) 
NDETt  =  PN(TSSt) 
PDETt  =  PP(TSSt) 

4.  NH3 


n  n 

NH3t  *  (NH30  -  Ft  -  %) 


T  n  T  n 

exp(-  kNH3tJ  +  h.  exp("kSSt)  +  ^2 


*  Equation  variables  are  defined  in  Table  2. 


Table  4  (Cont'd) 


-p-n  ^kNDET  ’  PM(TSSo  ~  -l 
where  =  - - - ~ - 

(kNH3  "  kSS) 

"  T  TSS  T 

^2  =  (kNDET  PN  *  -1  +  SnH3  "  kANH3A)/kNH3 

5.  N02 

N02t  =  (N02o  -  5l  “  ^2  '  53  ^  expf-kj^t)  +  5l  exp^ “kNH3t^ 


+  Ng  exp(-  kSS^  +  N.3 


where  ” 


t  _n  _n 

kNH3  (NH3o  '  h  ~  h] 


N02  '  *NH3 


T  ^ 
k  TT 
*NH3  -1 


k  _  k 

kN02  *SS 


53  - 


6.  N03 


T  ^ 
k  IT 
*NH3  -2 


'N02 


no  no  T  no  T  no  T 

N03t  =  N03o  +  TCj  -  ^2  exp(-k(J|02t)  -  expf-k^t)  exp(-k^t) 
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Table  4  (Cont'd) 


no  j  -r 

where  K4  =  ^ ^no2  53  +  $N03  “  ^AN03  ^ ^ 


no 

^  =  N02o  -  5l  -  52  "  53 


1°  =  ^h 


7.  P04 


NH3 


no  k 


J  -  **02  xr 

E,  -  —52 


kSS 


P04t  =  (P04o  -  -  1^)  exp(-kj04t)  +  kJ  exp(-kjst)  + 


TSS 


where  ^  = 


P  kPDET  <PPlTSS0  ‘Hi  ^ 


kT  -  kT 
kP04  kSS 


TSS 


T  i 

P  kP04  ’  PP  ’  5i  +  SP04  '  *AP04  •  A 


lP04 


8.  Dissolved  Oxygen  (0,) 

02t  =  02SAt[  -  Dt 


lNH3  =  3 >43  .  NH3t 
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Table  4  (Cont'd) 


L 


D 


1402 

t 


=  (0 


+*1 


+*3 


1.14  •  N02t 

_0  _D  _0  _D  _D 
o  _  1S_1  "  1S-2  ~  —3  ~  ^4  ~  exp(“^^T) 

IT  _! P  T 

exp(-(K^  +  k^lt)  +  Kg  exp(-kfJ1H3tl 

y  _P  j 

exp(-kN02  t)  exp(-  kS$° 


where 


k[  (L0  -  lx) 
kR  -  (kj  +  kLS) 


n  n 


3.43  k(|(H3  (NH30  -  K,  -  Kg)  +  1.14  k 


1  =2' 


kT  -  kT 
kR  kNH3 


1.14  kJQ2  (N02o 


kT  -  kT 
*R  *N02 


3,43  kl4H3  -1  +  1,14  kN02  52 


kT  -  kT 

kr  kss 


T 

1402 
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Table  4  (Cont’d) 


V 


*  3. 


43  *  1.14  kJ02j3  - 


02 


+  SS0D)/kR 


9.  Conservative  Attributes 
CONS1t  =  C0NSio 


10.  Nonconservative  Attributes 


NC0Nit  =  (NC0Nio  - 


sNC0Ni 


kNC0Ni*  exp^"kNC0Nit) 


SKCOMi 

*NC0Ni 


11.  Temperature 

Tt  =  (T0  -  KT )  exp(-kjt)  +  KT 
where 


kj  =  1.17  x  10-3  +  pL(a  +  bV)(gj  +  6.1  x  10"*  P) 

KT  =  (q$N  +  qat  -  7.36  x  10-2  -  pL(a  +  bV)(aj  -  ea 
-  6.1  x  10-4  p  .  AT)/kj 
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4  RATIONAL  THRESHOLD  VALUE  TEST  MODELS 


An  analytical  approach  such  as  that  provided  by  RIAS  requires  measurable 
indicators  of  impact  signif icance.  To  determine  such  significance,  threshold 
values  must  be  established.  Therefore,  it  is  necessary  to  develop  concepts 
for  using  RTVs  to  measure  the  significance  of  impacts  within  the  aquatic 
environment.  The  RTVTEST  models  used  to  develop  RIAS  are  a  subset  of  avail¬ 
able  model s.l  Impacts  can  be  analyzed  at  three  levels  of  effect  by  applying 
one  or  more  of  the  following  tests: 

1.  WQRTV  —  Assessment  of  the  extent  of  predicted  violation  of  existing 

ambient  water  quality  standards. 

2.  SIRTV  —  Assessment  of  the  effect  of  organic  pollution  on  the  micro¬ 

bial  community. 

3.  TURTV  —  Assessment  of  the  expected  concentrations  of  toxic  compounds 

on  overall  environmental  toxicity  of  receiving  system 

( species-speci fic) . 

The  flexibility  of  RIAS  allows  these  RTVTESTs  to  be  used  singly  or  in  conjunc¬ 
tion  with  each  other. 


Water  Quality  Standard  Assessments  (WQRTV) 

This  RTVTEST  quantifies  the  magnitude  of  water  quality  violations  which 
will  be  caused  by  the  impacts  of  the  project  being  evaluated.  This  test  is 
relatively  straightforward,  comparing  existing  stream  standards  to  the  output 
of  pertinent  SIMWQ  attributes.  Impact  is  quantified  in  terms  of  the  degree  of 
violation  (mg/*)  at  specific  points  in  the  WQ  (i,j,k)  profile,  the  spatial 
extent  of  violations  (mg/si/miles) ,  and  the  temporal  extent  of  violations.  The 
RTV  level  in  this  case  is  the  existing  water  quality  standard. 


Saprobic  Indices  (SIRTV) 


This  system2  was  developed  as  an  empirical  relationship  between  aquatic 
organisms  and  organic  water  pollution.  This  relationship  has  been  termed  the 
Saprobian  system  and  uses  the  concept  of  an  indicator  species.  The  applica¬ 
tion  of  the  empirical  relationships  of  the  Saprobian  system  in  a  quantitative 
index  was  introduced  by  Pantle  and  Buck3  and  expanded  by  Sladecek.4  They 

T  E.  E.  Herricks  and  M.  J.  Sale,  Development  of  Rational  Threshold  Values  for 
Aquatic  Ecosystems  (University  of  Illinois,  1978). 

2  R.  Kolwitz  and  M.  Marsson,  "Okologie  der  Pflanzlichen  Saprobien,"  Ber.  Dt. 
Bot.  Ges.,  Vol  26A  (1908),  pp  505-519. 

3  (1.  Pantle  and  H.  Buck,  "Die  Biologische  Uberwachung  der  Gewasser  und  die 
Darstellung  der  Ergebnisse,"  Gas.  Wass.  Fach.,  Vol  96,  No.  604  (1955). 

4  V.  Sladecek,  "The  Measures  of  Saprobfty,11  Verb.  Int.  Ver.  Limnol.,  Vol  17 
(1969),  pp  546-559;  and  V.  Sladecek  ."System  of  Water  Ouality  from  the  Bio¬ 
logical  Point  of  View,"  Erg.  Limnol.,  Vol  7  (1973),  pp  1-218. 
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describe  the  Saprobic  index  S  as  ranging  between  1  and  4  and  1  ani  3,  respec¬ 
tively.  Other  literature  has  shown  the  relationship  between  the  Saprobic 
index  and  BOD5  in  the  stream.5  Table  5  gives  the  relative  values  of  S  and 
BOD5.5  Saprobic  indices  have  generally  been  used  as  a  classification  scheme  in 
Europe,  but  have  not  been  widely  used  in  the  United  States. 

The  Saprobic  index  can  easily  be  calculated?  using  B11D5  concentrations 
provided  by  SIMWQ: 


,,  ^  =  l.U75(L( i  ,k) )  -  0.473 
*  "  0.218(L( i  ,k) )  +  0.904 

if  0  <  L( i ,k)  <  50  mg/ a 


[Eq  3] 


_  0 .0189 ( L < i  ,k ) )  -  7.938 
"  0.0021(L{ i  ,k) )  -  1.882 

if  L(i  ,k)  >  50  mg/  8. 


where  SI(i,k)  is  the  Saprobic  index  in  ith  point  in  space  and  kth  point  in 
time. 


The  values  of  SI(i,k)  can  then  be  used  to  interpret  the  impact  of  organic 
effluents  on  the  community  structure  of  the  receiving  stream.  For  example, 

RTV  levels  can  be  set  at  SI(i,k)  <  2.0  for  no  significant  impact  and  at  2  < 
SI(i,k)  _<  3  for  minimal  impact;  tTJen  the  output  from  the  SIRTV  routine  can  be 
used  to  quantify  the  extent  of  temporal  and  spatial  impacts  within  the  aquatic 
environment. 


Environmental  Toxicity  (TURTV) 

The  toxicity  unit  concept8  has  proven  to  be  a  useful  tool  for  Integrating 
biological  response  to  both  primary  toxicants  and  modifying  factors  (e.g., 
dissolved  oxygen  [DO],  temperature,  pH,  etc.).  This  index  has  been  used  suc¬ 
cessfully9  to  assess  the  biological  significance  of  water  pollution  Impacts. 

A  toxic  unit  of  a  specific  pollutant  is  simply  a  concentration  equal  to  the 

®  J .  Rothscheln,  "Saprobity  as  a  Criterion  of  Oxygen  Regime"  (in  Slovakian 
with  English  summary),  Pr.  Stud.  VUVH  Bratislave,  Vol  63  (1972),  pp  1-134; 
and  V.  Sladecek  and  F.  fucek,  "Relation  of  the  Saprobic  Index  to  BOD5,"  Wa~ 
ter  Res . ,  Vol  9  (1975),  pp  791-794. 

5  V.  Sladecek,  "The  Measures  of  Saprobity";  and  "System  of  Water  Quality  from 
the  Biological  Point  of  View." 

?  V.  Sladecek  and  F.  Tucek,  "Relation  of  the  Saprobic  Index  to  BOD5." 

8  K.  S.  Lublnski ,  R.  E.  Sparks,  and  l.  A.  Jahn,  Development  of  Toxicity  In- 
dices  for  Assessing  the  Quality  of  the  Illinois  River,  Research  Report  No. 

TJ IlU-VRr-7T-UOTb  (Uni versUy  of  I'll i'rioTr,' TT/4T." 

9  W.  V.  Brigham,  0.  A.  McCormick,  and  M.  J.  Wetzel,  The  Watersheds  of 
Northeastern  Illinois:  Quality  of  Aquatic  Environment  Based  (Jpon  tfater 
ffuaTity  and  Fishery  Data,  Staff  Paper  Wo.  31,  HIPC  (Illinois  Natural  History 
Survey,  1971]). 
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%-hour  LC5u  fo-  a  target  organism.*  Toxic  units  are  calculated  as  the  ratio 
of  simulated  ambient  concentrations  of  an  attribute  divided  bv  its  LC5Q.  This 
ratio  is  also  analogous  to  a  pollutant's  application  factor,11*  and  a  threshold 
level  can  be  specified  to  ensure  protection  for  target  organisms  (e.g., 
requiring  a  toxic  unit  _<  0.01  would  be  equivalent  to  an  application  factor  of 
100,  which  is  used  for  many  chlorinated  hydrocarbons).  Toxic  units  can  also 
be  accumulated  (summed)  for  all  potential  toxicants  to  obtain  an  overall  index 
of  environmental  toxicity.  Experience  has  indicated  acceptable  levels  of 
total  toxicity  units,11  but  one  must  remember  that  this  index  is  just  a  first 
approximation  of  biological  response. 

The  problems  involved  with  measuring  environmental  toxicity  include 
adjusting  for  the  effects  of  environmental  modifiers  and  predicting  synergis¬ 
tic  effects  of  various  pollutant  combinations.  Information  on  the  effects  of 
combinations  of  toxicants  is  available  for  only  a  few  species,  but  more  data 
are  being  collected  daily.  However,  despite  these  limitations,  the  toxic  unit 
model  is  still  the  best  general  model  now  available  for  impact  assessment. 

When  toxicity  data  are  available  at  multiple  trophic  levels,  the  toxic  units 
can  also  be  used  at  various  levels  to  measure  the  sensitivity  of  aquatic  com¬ 
munities. 


*  LC50  is  a  measure  of  the  concentration  level  of  the  toxic  material  that  will 
kill  50  percent  of  the  species  being  used  in  the  test  within  a  given  time 
interval  (e.g.,  96  hours). 

10Quality  Criteria  for  Water  (U.S.  Environmental  Protection  Agency  [USEPAj, 

imr —  - 

^R.  Lloyd  and  0.  H.  M.  Jordan,  "Predicted  and  Observed  Toxicities  of  Several 
Sewage  Effluents  to  Rainbow  Trout,"  J.  and  Proc.  Inst.  Sew.  Purif.  (Brit.), 
Vol  2  (1963),  pp  183-186. 


Table  5 


Values  of  S  and  BOD5  (=L)  for  Upper  Limits 
of  Individual  Saprobic  Degrees* 


Degree 

S 

L 

Note 

Katharobity 

-0.5 

0.0 

Purest  water 

Zernosaprobity 

0.5 

1.0 

Very  clean 

01 igosaprobity 

1.5 

2.5 

Clean 

Beta-mesosaprobi ty 

2.5 

5.0 

Mild  pollution 

A1 pha-mesosaprobi ty 

3.5 

10.0 

Pollution 

Polysaprobity 

4.5 

50.0 

Heavy  pollution 

Isosaprobity 

6.5 

400.0 

Sewage 

Meta saprobity 

6.5 

700.0 

Septic 

Hypersaprobity 

7.5 

2,000.0 

Putrefaction 

Ultrasaprobity 

8.5 

120,000.0 

Lifeless  liquors 

♦From  V.  Sladecek,  "The  Measures  of  Saprobity,"  Verb.  Int.  Ver.  Limnol., 

Vol  17  (1969),  pp  546-559;  and  V.  Sladecek  "System  of  Water  Quality  from  the 
Biological  Point  of  View,"  Erg.  Limnol.,  Vol  7  (1973),  pp  1-218 . 


5  APPLICATIONS 


This  chapter  provides  a  series  of  application  scenarios  which  demonstrate 
the  utility  of  RIAS.  Appendices  A  and  S  provide  sample  output  of  the  RIAS 
system  application.  These  examples  concentrate  on  problems  pertinent  to 
current  Army  activities,  but  are  not  meant  to  be  exhaustive.  The  role  of  RIAS 
or  a  similarly  designed  methodology  in  improving  impact  planning  will  become 
apparent  in  this  discussion. 

The  general  framework  for  decision-making  using  RIAS  should  be  envisioned 
as  an  iterative  process  using  the  three  computerized  procedures  --  FQUES, 
SIMWQ,  and  RTVTF.ST  --  as  the  central  tools.  Adequate  methods  for  impact  scop¬ 
ing  presently  exist  in  the  form  of  matrix  methodology  such  as  the  Environmen¬ 
tal  Impact  Computer  System  (EICS),*2  a  computerized  system  to  help  environmen¬ 
tal  planners  identify  and  mitigate  impacts  of  proposed  Army  projects  or 
activities.  However,  a  shortcoming  of  matrix  methodology  is  that  it  is  not  a 
data-handl  ing  tool  and  provides  little,  if  any,  quantification  potential. 

This  is  the  purpose  of  developing  secondary  algorithms  like  RIAS  for  impact 
assessment.  While  computer-based  matrix  methodologies  work  well  for  specify¬ 
ing  a  project's  potential  impacts,  they  cannot  detect  small  changes  in  project 
specifications  which  may  mean  the  difference  between  significant  or  nonsigni¬ 
ficant  impacts.  RIAS,  which  has  this  capability,  uses  quantitative  procedures 
to  integrate  impacts  according  to  changes  in  project  specifications.  Environ¬ 
mental  setting  data  are  used  to  predict  impact  magnitudes.  Subsequent  to  ini¬ 
tial  use  of  these  procedures,  project  specifications  can  be  changed,  allowing 
impact  assessment  of  various  alternatives  to  be  done  easily  and  inexpensively. 
In  this  way,  the  RIAS  methodology  provides  a  truly  quantitative  tool  for 
impact  management. 

Besides  quantification,  another  major  advantage  of  RIAS  is  standardiza¬ 
tion.  The  computer-based  algorithms  define  the  organization  of  environmental 
information  and  specify  its  use  in  a  consistent,  repeatable  protocol.  At  this 
time,  the  links  between  primary-level  impact  assessment  (i.e.,  scoping  activi¬ 
ties  like  EICS)  and  secondary  procedures  such  as  RIAS  should  be  via  paper  ties 
only.  This  could  be  in  the  form  of  RAMIT  (ramification/mitigation)  statements 
output  from  EICS,  descriptor  package  writeups,  and  user  manuals  such  as  the 
appendices  to  this  report.  Ultimately,  if  a  system  such  as  RIAS  received 
broad-based  support,  data  sets  specifying  environmental  setting  could  be 
assembled  for  all  Army  bases,  providing  "off  the  shelf"  assessment  capability 
whenever  new  missions  altered  base  operations.  Since  the  goals  of  impact 
modeling  should  be  flexibility  and  easily  usable  prediction  tools  coupled  to 
readily  available  data  sets,  RIAS  can  be  seen  as  a  design  prototype  for  this 
type  of  methodology. 

The  following  sections  provide  examples  of  how  RIAS  can  be  applied  to 
assessing  environmental  impacts  resulting  from  Army  projects. 


12r.  Baran  and  R.  0.  Webster,  Interactive  Environmental  Impact  Computer  System 
(EICS)  User  Manual ,  Technical  Report  N-BO/AUAUMOyU  (U.5.  Army  Construction 
Engineering  Research  Laboratory  [CERL],  September  1979). 


Example  1 


The  impacts  occurring  in  aquatic  receiving  systems  which  are  the  most 
frequently  studied  originate  from  point  source  discharges  of  domestic  sewage 
effluents.  This  type  of  effect  must  be  considered  in  environmental  impact 
statements  and/or  assessments  of  Army  activities.  Receiving  stream  impacts 
caused  by  organic  pollution  occur  when  mission  changes  alter  sewage  treatment 
plant  loadings  or  treatment  efficiencies  or  when  old  sewage  treatment  plants 
are  upgraded.  RIAS  provides  an  assessment  tool  for  analyzing  these  situa¬ 
tions.  The  BOD/DO  models  within  RIAS  are  only  a  subset  of  a  more  general 
model.  The  outputs  from  both  WQRTV  and  SIRTV  quantify  these  types  of  impacts 
and  place  impact  predictions  in  easily  understandable  terms.  (See  Appendix  B, 
Part  3. ) 


Example  2 

Waste  discharges  from  many  Army  industrial  or  laboratory  activities  con¬ 
tain  toxic  components  which  can  adversely  affect  aquatic  biota.  The  magnitude 
and  spatial  extent  of  these  impacts  are  largely  a  function  of  environmental 
setting,  such  as  watershed  dilution  capacity,  ambient  water  quality,  and  local 
target  species.  The  TURTV  routine  in  RIAS  provides  a  consistent  method  of 
considering  such  information  and  quantifying  impact  magnitude  using  the  toxic 
unit  concept.  (See  Appendix  A,  Part  3,  Section  C.) 


Example  3 


Vehicle  maintenance  activities  account  for  point  source  discharges  of 
many  potentially  harmful  water  quality  attributes,  including  suspended  solids, 
detergents,  oils  and  greases,  and  general  BOO.  While  washrack  facilities  are 
being  redesigned  and  relocated,  the  impacts  of  these  facilities  within  a 
watershed  context  could  be  evaluated  using  the  RTVTESTs  in  RIAS.  Thus,  site- 
specific  design  activities  could  be  made  more  efficient  by  considering  the 
watershed  assimilation  capacity  of  sensitivity  (Appendix  B,  Part  3)  as  shown 
by  RIAS. 


Example  4 

Impacts  from  landfill  leachate  are  increasingly  affecting  the  aquatic 
environment.  The  significance  of  these  impacts  depends  on  many  factors, 
including  the  water  body's  assimilative  capacity  and  the  sensitivity  of  biota 
in  the  receiving  watershed.  RIAS  is  an  excellent  tool  for  evaluating  the 
severity  of  these  impacts.  If  leachate  rates  and  initial  concentrations  can 
be  estimated  and  isolated  in  a  watershed,  a  simulated  point  source  discharge 
can  be  created  to  represent  leachate  inputs.  (See  Appendix  B,  Part  3.) 
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6  USER'S  GUIDE  TO  RIAS 


Organization  of  Information  Inputs 

The  first  step  in  using  the  RIAS  computer  routines  for  an  impact  assess¬ 
ment  is  organizing  the  data  sources  and  describing  the  problem.  This  requires 
several  stages  of  data  collection  and  organization,  including:  (1)  identifi¬ 
cation  of  the  control  parameters  for  simulation  and  assessments,  (2)  collec¬ 
tion  of  data  on  boundary  conditions  for  the  analysis,  and  (3)  estimation  of 
kinetic  rate  coefficients  and  source/ sink  terms  for  the  water  quality  simula- 
ti ons . 

The  first  step  is  determining  the  impact  types  to  be  considered  and  the 
geometric  description  of  the  watershed  to  be  analyzed.  Potential  impacts  must 
have  already  been  identified  by  some  type  of  scoping  procedure  (e.g.,  EICS). 
Physical  and  chemical  water  quality  attributes  are  required  at  this  point,  as 
well  as  environmental  modifiers  which  might  be  important  in  weighing  impact 
significance.  This  information  is  used  to  specify  the  attributes  to  be 
modeled  in  SIMWQ.  The  watershed  description  consists  of  identifying  reach 
lengths,  drainage  areas  upstream  from  the  top  of  each  stream  reach,  tributary 
and  effluent  locations,  and  bifurcation  structure.  (Chapter  2  provides  cri¬ 
teria  for  specifying  reaches.)  Drainage  areas  are  calculated  using  standard 
U.S.  Geological  Survey  maps. 

The  technique  to  be  used  for  numbering  reaches,  bifurcations,  tribu¬ 
taries,  and  effluents  is: 

1.  Numbering  reaches:  Number  reaches  beginning  at  the  top  of  the  most 
upstream  minor  branch  of  the  stream  to  be  modeled.  Proceed  downstream  until  a 
confluence  is  encountered.  Skip  to  the  top  of  the  next  most  upstream  minor 
branch  and  continue  downstream  to  the  next  confluence.  When  no  more  minor 
branches  remain,  proceed  down  the  main  branch  from  upstream  to  downstream. 

2.  Numbering  bifurcations:  Each  channel  bifurcation  (confluence  loca¬ 
tion)  is  designated  by  a  real  number  consisting  of  digits  in  the  tens,  ones, 
tenths,  and  hundredths  places.  For  example,  for  "10.05,"  the  whole  number 
part  of  the  indicator  ("10"  in  this  example)  represents  the  receiving  reach 
downstream  of  the  confluence.  The  fractional  part  of  the  indicator  ("5"  in 
this  example)  represents  the  last  reach  of  the  minor  branch  which  is  entering 
a  higher-order  stream. 

3.  Numbering  tributary  Inputs:  Using  the  numbering  system  described 
above,  number  upstream  tributaries  first,  doing  the  more  minor  branches  first. 
This  convention  is  not  critical,  but  will  provide  more  consistency.  Remember 
that  the  inputs  at  the  top  reach  of  each  branch  must  be  designated  as  a  tribu¬ 
tary  in  order  to  set  boundary  conditions  for  the  simulation  model.  The  index 
of  each  tributary  designates  the  stream  reach  into  which  it  empties,  not  sim¬ 
ply  its  number. 
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4.  Numbering  effluent  inputs:  Using  the  numbering  system  described 
above,  number  effluents  in  an  upstream  to  downstream  manner.  As  with  tribu¬ 
taries,  the  index  number  of  an  effluent  represents  the  number  of  the  reach 
into  which  it  empties,  not  the  number  of  the  effluent. 

The  second  stage  of  data  organization  is  specifying  boundary  conditions 
of  water  quality  attributes  at  tributaries  and  effluents  within  the  watershed. 
Techniques  for  doing  mass  balances  on  Army  installations  are  available  which 
are  adequate  for  describing  effluents.13  Another  source  of  effluent  informa¬ 
tion  is  the  National  Pollutant  Discharge  Elimination  System  (NPDES)  permits 
for  point  source  discharges.  Many  sources  are  available  that  describe  infor¬ 
mation  on  tributary  inputs.1^  However,  when  no  information  on  boundary  condi¬ 
tions  is  available,  field  data  must  be  collected. 

The  last  stage  of  data  organization  is  specifying  the  kinetic  terms  in 
SIMWQ  equations.  This  is  the  most  important  requirement  for  ensuring  accurate 
impact  assessment.  The  structure  used  in  SIMWQ  takes  advantage  of  a  class  of 
widely  used  simulation  models  (steady-state,  plug  flow)  whose  parameters  are 
well  understood.  The  review  of  Zison,  et  al . ,  is  a  good  primary  source  of 
information  on  the  state  of  the  art  of  estimating  these  model  parameters.1^ 
Table  2  provides  a  range  of  values  experienced  for  all  model  terms  and  refer¬ 
ences  to  previous  modeling  work  in  which  they  were  used. 


Example  Problem 

To  illustrate  the  use  of  all  the  RIAS  routines,  Appendix  B  provides  an 
example  impact  assessment.  Figure  3  shows  the  layout  of  a  hypothetical  Army 
post  where  the  potential  impacts  from  four  point  sources  of  pollutants  on 
aquatic  receiving  systems  will  be  analyzed.  Pre-analysis  information  has 
already  been  organized,  and  the  example  begins  as  FQUES  is  executed  to  build 
up  the  project  specification/environmental  setting  data  base.  All  user 
responses  in  the  computer  output  have  been  underlined. 


13G.  W.  Schanche,  et  al.,  Water/Wastewater  Survey  Guidelines,  Technical  Report 
N-U/ADA033223  (CERL,  1975X1 

^E.  E.  Herricks  and  M.  J.  Sale,  Development  of  Rational  Threshold  Values  for 
Aquatic  Ecosystems  (University  of  Illinois, 

15S.  W.  Zison,  et  al,.  Rates,  Constants  and  Kinetic  Formulations  In  Surface  Wa- 
ter  Quality  Modeling,  EPA-t>uu/3-/8-Un?  (USEPA,  Iv/d). 
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*▼11111  U I  / 

Reservation 

Boundaries 


/ 


Point  Source  Effluents 
■  Sewege  Treatment  Plants 

▼  Vehicle  Maintenance  Area 

•  Solid  Waste  Disposal  Site 


Reach  5  / 

Reach  6 

Reach  7^ 


Figure  3.  Physical  layout  for  example  problem. 


7  CONCLUSION 


This  report  has  described  the  framework  and  demonstrated  the  utility  of 
RIAS,  a  computerized  technique  which  uses  the  concept  of  rational  threshold 
values  to  determine  Impact  significance.  This  system  uses  a  rigorous  set  of 
analysis  procedures  to  identify  significant  effects  resulting  from  Army 
activities  on  aquatic  ecosystems. 
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APPENDIX  A: 


ANALYTICAL  SOLUTIONS  USED  IN  SIMWQ 


The  general  solution  used  in  the  SIMWQ  equations  is  demonstrated  here  for 
a  simple  steady-state  situation  of  coupled  water  quality  attributes.*  The 
equations  developed  here  are  used  only  within  a  stream  reach  in  which  all 
model  parameters  are  constant.  Implicit  assumptions  are: 

1.  Nondispersive,  plug  flow 

2.  Temporal  steady-state  conditions  at  all  upstream  loading  points, 
source/sink  terms,  and  stream  discharge 

3.  First-order  reaction  kinetics. 


The  initial  system  of  coupled  reaction  is  given  by: 


dAt 

~dt 


-kiAt  +  SA 


[Eq  Al] 


dBt 

dt 


-k2Bt  +  k3At  +  SB 


where: 

At  =  concentration  of  attribute  A  at  time  t 
B^  =  concentration  of  attribute  B  at  time  t 
The  general  solution  is: 

SA 

At  -  sAk!  +  (Ao  “ 


[Eq  A2] 


Bt  -  C 


k3sA 

klk2 


+ 


SB 

kiT 


+ 


[B( 


k3sA 


kl<k2  "  kl)  - 


k3Ao 
k2  -  1 


k3SA 

klk2 


SB 

*2 


3  exp(-k2t) 


+  [ - — - - 3  exp{-kit) 

lk2  -  ki»  -  mf^rr 


*  A  more  detailed  discussion  of  this  type  of  equation  development  can  be  found 
in  Chapters  4  and  5  of  R.  V.  Thomann,  Systems  Analysis  and  Water  Quality 
Management  (McGraw-Hill,  1972). 
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Generally,  this  type  of  solution  takes  on  the  form: 

n 

Ct  «  a  +  r  Bi  exp(-Yjt)  [Eq  A3] 

i=l 

where  a  ,  &j,  and  Xi  are  constants  within  a  reach  which  can  be  calculated  from 
data  inputs,  and  C*  is  the  concentration  of  any  attribute  at  time  t. 

In  some  cases,  this  analytical  solution  degenerates;  i.e.,  if  k2  is  input 
as  zero  or  if  k£  =  k^,  the  result  will  be  an  illegal  arithmetic  operation. 
Therefore,  contingencies  have  been  made  in  the  FORTRAN  code  to  avoid  this  type 
of  solution  breakdown.  Using  the  coupled  system  defined  initially,  four  con¬ 
ditions  can  be  identified  which  cause  a  breakdown: 

1.  k^  =  0,  k2  nonzero 

2.  k2  =  0,  kj  nonzero 

3.  ki  =  k2  =  0 

4.  k2  ■  ki  i  0. 

The  solutions  used  to  avoid  blowup  in  SIMWQ  can  be  derived  as  follows: 
Condition  1  (kj  =  0,  k2  i  0): 

At  =  a0  +  SAt  [Eq  A4] 

♦  jjfi  *  <*o  - k-^  ♦  hv2 

Condition  2  (k2  =  0,  k^  1  0): 

At  ■  ^  *  (Ad  -  j£>  Mp(-klt)  [Eq  A5] 

Bt  *  -  ^3  exp(-klt)  ♦  [SB  ♦ 

Condition  3  (kj  =  k2  »  D): 

At  •  A0  +  SAt  [Eq  A6] 

Bt  =  B0  +  SQt 
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[Eq  A7] 


Condition  4  (kj  =  k£  i  0): 


At  =  ^  +  (Aq  -  exp(-kit) 

k^S/v  Sr  k^Si  Sr 

Bt  *  ♦  -i)  +  (B0  -  J-p.  -  J-  exp(-k2t) 

kjk2  K2  k2kj  k2 


+  k3  ( A0  -  -l)t  exp(-k].t) 
ki 


This  same  type  of  logic  and  equation  formulation  is  applied  to  all  sets 
of  coupled  equations  in  SIMWQ.  Because  of  the  additional  terms  compounded  in 
t,  the  form  of  the  general  solution  is  expanded  to 

n 

C  =  «  +  E  expi-y^t)  +  fit  +  et  exp(pt)  [Eq  A8] 

i=l 


As  before,  all  the  constant  terms,  a,6,Y,5,e,p,  can  be  calculated  based  on 
input  for  each  reach  and  C  could  be  the  concentration  of  any  attribute.  These 
terms  are  all  calculated  in  the  subprogram  APARAM,  which  also  includes  the 
logic  to  avoid  unnecessary  breakdowns  in  the  computer  software  codes.  This 
logic  increases  the  usability  of  the  final  product  and  requires  a  minimum  of 
user  expertise. 
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APPENDIX  B 


RIAS  USE  EXAMPLE 

Part  1:  New  Data  File  Creation 


TERMINAL:  257 

79 / 1 2/ Id .  10.32.93.  ' 

UNIVERSITY  Of'  ILLINOIS  CYBER  175.  NOS  1  .  3  -  485/985. 

SIGNON :  341447562 
PASSWORD 

MIHMI 

TERMINAL:  257,  TTY 

RECOVER/  CHARGE:  bill ■ eeusa ■ os7770 
LAST  RECORDED  SIGNON  AT  10:23  12/18/79 
/-fquea 


A ###»#«##»#«»##***#**# A###*#*#####**### 

THIS  PROGRAM  ALLOWS 
THE  USER  TO  BUILD  UP  A  NEW 
DATA  FILE  OR  TO  REVISE  AN  OLD  DATA 
FILE  FOR  SUBSEQUENT  CONTROL  AND 
INPUT  FOR  WATER  QUALITY  SIM- 
ULATONS  UNDER  'RIAS' 

t 


DO  YOU  WISH  TO  BEGIN  CREATING  A  NEW  DATA  FILE 
(ANS:  YES  OR  NO)  7  X 

TIME  INVARIANT  PARAMETERS 

I)  TYPE  IN  THE  NAME  OF  THE  DATA  SET 
7  example  no.  1 
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XXI) 


ANSWER  THE  rai-I.JWX.VO  JUESTIJ.'IS  A30UT  S IMULA 71  ON 
CONTROL  PARAMETERS. 

HOW  MANX  REACHES  (AMS: 1-20)  ?  3_ 

HOW  MANX  TIME  PERIODS  (ANS  :'-'2T  ?  2_ 

HOW  MANX  ADDITIONAL  CONSERVATIVE  WATER  OUALITX 
ATTRIBUTES  3EX0ND  THE  CORE  ATTRIBUTES 
(ANS:  0-’2>  ?  j. 

HOW  MANX  ADDITIONAL  NOMCONSER VATI VE  ATTRIBUTES 
(ANS:  0-4)  ?  2_ 

INPUT  THE  NAMES  OF  THESE  ADDITIONAL  PARAMETERS. 
CONSERVATIVE  ATTRIBUTES: 


ATT  . 

NO  . 

9)7 

hard  . 

ATT  . 

NO  . 

10) 

?  ££ 

ATT  . 

NO  . 

’  ’  ) 

?  z  n 

NONCONSERVAT 17 

£  A 

:  a  _  w  -  r 

ATT  . 

NO  . 

1  2  ) 

?  cl  2 

ATT  . 

NO  . 

1  3 ) 

?  To  »  ; 

INPUT  THE  MODELING  CODES  FOR  THE  EIGHT  CORE 
ATTRIBUTES  (ANS:  0  OR  1  ) 


INPUT  THE  PARAMETERS  DESCRIBING  THE  WATERSHED 
STRUCTURE  FOR  THIS  SIMULATION. 


HOW  MANX  TRIBUTARX  INPUTS  7  5 
INPUT  CODES: 


£  1  :  £ 

3) ?  i 

4) ?  7 

5) ?  I 

HOW  MANX  POINT  SOURCE  DISCHARGES  ?  £. 
INPUT  CODES: 

1  )?  2 
2)?  “ 

3)?  T 
4  )?  2. 

HOW  MANX  MfTlN  BRANCH  BIFURCATIONS  ? 
BIFURCATION  CODES: 

1 )  7  5.0’ 

INPUT  THE  LENGTH  OF  EACH  REACH  AND  THE 
DRAINAGE  AREA  UPSTREAM  FROM  THE 
TO?  OF  EACH  REACH  (ANS:  MILES  AND 
SQUARE  MILES). 

1  )  7  3  .  2  .  1  .  3 

3) 7  2.5. 4.9 

4) 7  2 ■ 2  .  6  .  3 

5) 7  " 

6>?  2.’, '3-0 
7)7  5.7,33- 

8  )  7  20 .  , 99 J  ■ 

INDICATE  HOW  XOU  WISH  TO  MODEL  STREAM 
HXORAULICS. 

0)  MEAN  DEPTH  AND  VELOCITX  SPECIFIED 
FOR  EACH  REACH  AND  TIME  PERIOD. 

1)  HXDRAULIC  HATING  PARAMETERS  USED  FOR 
EACH  REACH. 

(ANS:  USE  EITHER  0  OR  ’ )?  0 
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ZERO  TIME  VARIANT  PARAMETERS  FIRST  (Y  OR  N  )  ?  n 


TIME  VARIANT  PARAMETERS  FOR  TIME  PERIOD  NO.  1 


I)  HIDRA’JLIC  PARAMTERS . 

INPUT  MEAN  VELOCITIES  FOR  EACH  REACH. 
?  .  2» ,  ■ 1 7 ,  . 23 .  • 25 ,  . 29 ,  ■ 3 ,  . 3d,  ■ 93  • 

INPUT  MEAN  DEPTHS  FOR  EACH  REACH. 

?  . 16, . 14, .2,  .22. .25 . 

,.25,..  <  ERROR,  RETYPE  RECORD  AT  THIS  FIELD 
25,  -33  ,  -3 

II)  BOUNDARY  CONDITIONS  AT  TRIBUTARIES. 


INPUT  AMBIENT  RATER  OJALITY  CONDITIONS  FOR 
TRIBUTARIES  1  THROUGH  5 


TEMP  . 

,20 ■ 1 , 20  .  , 25 ■ 2 , 27 
B0D5 
TSS 
NH3 
N02 
NO  3 

P04 
-D.O  . 

HARD  . 

PH 
ZN 
CL2 
COLI  . 


20 . 1 , 20 . 1 , 20 . 1 , 25 . 2, 25. ; 


•DEL* 


?  .5, • 5,  .5,  ’ .0,  ’  .2 
?  IQ  ■  ■  1  0  ■  ■  1  l'.  p  ■  2i  25  ■ 
?  0 . 0  i.O  ,  J  ,  0 

7  0 . 0 . 0  ■  0  ■  0 
7  1.3. 1.3. .5, 2. 6. 2.  8 
?  ^0Sju.05,-9?.-J.^. 

7  3. 25. 8.25 , i ■ 25 , 7. 6, 

7  iso',  220 , 220',  295 , 295 
7  7.9,7. j. 7. 75,7. 9,779 
7  0 ■ , 0 ■ , 0,0,0 
7  0 , 0  ,  0 , 0 , 3 
?  0  ■  0 . 0  ■  0  ,  0 


7.9 


MEAN  DISCHARGE?  . 0 3 ■  ■ 09 .  . 05  ■  ■  5 , 28 


III)  BOUNDARY  CONDITIONS  AT  EFFLUENT  DISCHARGES. 


INPUT  AMBIENT  RATER  QUALITY  CONDITIONS  FOR 


EFFLUENTS 

THROUGH  4 

TEMP. 

o 

27. ,25. 4.22. ,25 

B0D5 

? 

5.9,136, 200 , . 5 

TSS 

o 

3.  '24,250,  ’2 

NH  3 

? 

.5,22.6, '.5, 10. 

N02 

7 

0. 0.0.0. 

7  9 . 8 , 0 . 2 , 0 


MEAN  DISCHARGE?  2..3..’.,.2 


REACTION  RATE  COS? FI  Cl  ENTS . 


INPUT  THE  INDICATED  RATS  COEFFICIENT  FOR 
REACHES  1  THROUGH  3 


?  .05,  .05.  .05,  . 05 • .  • 05 .  .05.  .03. -05 


,.05. .05.  <  ERROR,  RETIRE  RECORD  Al 
.05 


THIS  FIELD 


KANO  3 

KAP04 

KL 

KLS 

KNDST 

KN  H  3 

KN02 

KNO  3 

KPDET 

KP04 

KR 

KS3 

KSSS 

KT 

KTS3 
XtICA  1 
KNCA2 


.  05 .  ■ 05, . 05.  . 05 ,  . 05,  . 05 ,  ■ 05 .  ■ 05 


.  1  6  .  .  ’  6  ,  .  '  6  ,  .  1  a  ,  .  16  ,  .  ’  6  ,  ,  ’  6  ,  6 

.O'  ,  ,01  ,  .0-  ■  .0*2,0  1  .  .0'  ■  -0’  ■  -.O’ 

■  02.  .02.  . 02,  .02,  .02,  .02,  .02,  .02 

.  65  .  .  5  5  ■  .  65  6  5  .  ■  0.5  .  .  65  ,  .  65  ,  .  o_j_ 

3.5.2  -.5  ■  2  .  5  ,  2  ■  5  ■  2,-  5  ■  2  ■  5 , 2  .  5  ■  2,-  3 
.001 . .OP" , , 3’ ' . .00* . .00’ . .0' i . . 30 ^ 
.  02.  .  02,  .  02^  .  02,  J2,  .  02  ■  ■  02  ■  ■  02 


?  ,  001  ,  ■  00  -1  ,  ,  0  ’  1  .  .  0  0  '  .  .  0  0  ’  .  ■  0 

?  .02. .02. .02. .02. ,02. .02. .02 
?  -2,  •2.-2,  ,2,  .2,  •  2^.-2.  ,2- 

1  0,0 ,0,0, 0.0, 0.0,0. 

^  .  0  1  ,  .  0  1  .  .  0  1  ,  .  0  ' ,  .  0  1  ,  .  0  *  ,  .  0  ’ 

?  0.0, 0,0, 0.0, 0,0 


?  .  02  .  .02.  .02,  .02,  .02,  ,02,  .0  2,  .02 

?  0,0, 0,0, 0,0, 0,0 

?  0.0,0, 0,0, 0,0,0 


V)  DISTRIBUTED  SOURCE/SINK  PARAMETERS. 

INPUT  THE  INDICATED  SOURCE/SINK  TERM  FOR 
REACHES  1  THROUGH  3 


Part  2:  Data  File  Revision 


TERMINAL:  46 

79/  '  2/ ’ d .  ^  3  •  “2 . 04 . 

UNIVERSITY  OF  ILLINOIS  CYBER  175. 

SIGNON :  3« ’447562 , San? 

TERMINAL:  46,  TTY 

RECOVER/  CHARGE:  61  U  .  causa , ps 7770 
LAST  RECORDED  SIGNON  AT  1’:55  l2/'B/79 
/net , tapa9 
/-fauas 


JOB  ACTIVE. 


NOS  ’ . 3  -  485/485 . 


THIS  PROGRAM  ALLOWS 
THE  USER  TO  3UILD  UP  A  NEW 
DATA  FILE  OR  TO  REVISE  AN  OLD  DATA 
FILE  FOR  SUBSEQUENT  CONTROL  AND 
INPUT  FOR  WATER  QUALITY  SIM- 
ULATOMS  UNDER  'RIA3' 


DO  YOU  WISH  TO  BEGIN  CREATING  A  NEW  DATA  FILE 
'ANS:  YES  OR  NO)  7  n 


DO  YOU  WISH  TO  REVIEW  AND/OR  REVISE  THE  EXISTING 
DATA  FILE 

(ANS:  YES  OR  NO)  ?  y 


REVIEW  OF  CONTENTS  OF  CURRENT  DATA  FILE 


THE  CURRENT  CONTENTS  OF  YOUR  DATA  FILE  HAVE  THE 
•  TITLE: 

EXAMPLE  NO.  1 


THIS  DATA  SET  SPECIFIES  SIMULATION  OF  THE 
FOLLOWING  WATER  QUALITY  ATTRIBUTES: 
TEMP.  B0D5  TSS  NH  3  N02 

NO  3  P04  D.O.  HARD.  PH 

ZN  CL2  COLI. 

SIMULATIONS  WILL  BE  RUN  FOR  2  TIME  PERIODS 
FOR  A  TOTAL  OF  0  STREAM  REACHES. 
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THS 

FIED  WATER  C- , 

(ED  STRUCTURE 

15  A3 

FOLLOWS : 

H£ACH 

LENG  TH 

DRAINAGE 

INPUTS 

SO. 

(MI.) 

(SO. MI . ) 

EFF 

TRIB 

1 

3 . 200 

'  .  30  0 

4 

3 

2 

.  dOO 

1  .  300 

i 

2 

3 

2 .500 

u  .  90  0 

•  t 

’ 

A 

2. 200 

6  .  300 

2 

•  • 

5 

3  •  500 

11.100 

i  • 

•  • 

6 

2.100 

1 3 .000 

3 

*  • 

7 

5.700 

32.000 

•  • 

4 

3 

20.000 

990 . 000 

a  a 

5 

UMBER  OF  MAJOR  BIFURCATIONS  OF  THE  MAIN 
CHANNEL  OF  THIS  RECIEVING  STREAM  IS  1 


BRANCH  NO.  ’  INCLUDES  REACHES  1  TO  ' 

THE  MAIN  CHANNEL  INCLUDES  REACHES  2  TO  8 


HYDRAULIC  MODELING  WILL  8E  DONE  USING  MEAN 

VELOCITIES  AND  DEPTHS  FOR  EACH  REACH  AND  TIME 
PERIOD. 

NONE  OF  THE  ABOVE  PARAMTERS  CAN  BE  ALTERED  WITHOUT 
CREATING  A  TOTALLY  NEW  DATA  SET  (I.E.,  BY 
STARTING  OVER  WITtf  "  FQUES  * ) . 

DO  YOU  WANT  TO  CONTINUE  (ANS:  YES  OR  NO)?_y_ 

WHICH  TIME  PERIOD  DO  YOU  WANT  TO  REVIEW  ?  2 


INDICATE  WHICH  OF  THE  FOLLOWING  PARAMETERS  YOU 


WANT 

TO 

REVIEW  (RESPOND  WITH  THE  PROPER  NUMBER 

1  ) 

HYDRAULIC  PARAMETERS 

2) 

INITIAL  CONDITIONS  IN  TRIBUTARIES 

3) 

INITIAL  CONDITIONS  III  EFFLUENTS 

4) 

KINETIC  PARAMETERS 

5) 

DISTRIBUTED  SOURCE/SINK  PARAMETERS 

?  4 

b) 

BIOLOGICAL  PARAMETERS 

INDICATE 

THE 

INDEX  NUMBER  OF  THE  KINETIC  PARAMETER 

YOU 

ARE 

INTERESTED  IN  REVIEWING 

. 

1  ) 

KANH  3 

2) 

KAN03 

3) 

KAP04 

4) 

KL 

5) 

KLS 

6) 

KNDET 

7) 

KH  H3 

d) 

KN02 

9) 

KN0  3 

10) 

KPDET 

1  ’  ) 

KP04 

12) 

KR 

13) 

KSS 

14) 

KSSS 

15) 

KT 

16) 

KTSS 

17) 

KNCA  ’ 

19) 

KNCA2 

WHICH  ONE 

7  ’ 

a 

THE  VALUES  CURRENTLY  SPECIFIED  FOR  YACA2  WILL 

BE  LISTED  3EL0W  31  REACH.  TO  CHANCE  A  VALUE 
RESPOND  TO  THE  TRAILING  '?'  WITH  THE  NEW 
VALUE. 


1 ) 

0  .  ? 

.  002 

2) 

0  .  ? 

.  002 

3) 

0.  ? 

.  002 

A  ) 

0  .  ? 

.  002 

5) 

0.  ? 

.  002 

b> 

0  .  ? 

.002 

7) 

0.  ? 

TocST 

S) 

0.  ? 

.  002 

REVIEW  ANOTHER  RATE  COEFFICIENT  CANS:  YES  OR  NO)?_n_ 

CONTINUE  FOR  THIS  TIME  PERIOD  (ANS:  YES  OR  NO)?  n 

REVIEW  ANOTHER  TIME  PERIOD  (ANS:  YES  OR  NO)  ?  y 

WHICH  TIME  PERIOD  DO  YOU  WANT  TO  REVIEW  ? 

INDICATE  WHICH  OF  THE  FOLLOWING  PARAMETERS  YOU 

WANT  TO  REVIEW  (RESPOND  WITH  THE  PROPER  NUMBER) 

1)  HYDRAULIC  PARAMETERS 

2)  INITIAL  CONDITIONS  IN  TRIBUTARIES 

3)  INITIAL  CONDITIONS  IN  EFFLUENTS 

4)  KINETIC  PARAMETERS 

5)  DISTRIBUTED  SOURCE/SINK  PARAMETERS 

6)  BIOLOGICAL  PARAMETERS 


...DICATS  THE  INDEX  NUMBER  OF  THE  DISTRIBUTED 
SOURCE/SINK  TERM  YOU  WISH  TO  REVIEW 


1  ) 

SL 

2) 

SNH3 

3) 

SNO  3 

«) 

S02 

5) 

SP04 

t>) 

SSOD 

7) 

SSS 

B) 

ST 

9) 

SAT  1 

10) 

3AT2 

1  1  ) 

SAT3 

12) 

SA  T 4 

13) 

SAT5 

? 

THE  VALUES  CURRENTLY  SPECIFIED  FOR  SSOD  WILL 

BE  LISTED  BELOW  BY  REACH.  TO  CHANGE  A  VALUE 
RESPOND  TO  THE  TRAILING  'T  WITH  THE  NEW 
VALUE . 


1  ) 

0.7 

.6 

2) 

0 . 7 

3) 

0.7 

4) 

0 . 7 

■  6 

5) 

0.7 

.  0 

6) 

0.7 

.  6 

7) 

0. 7 

.  3 

8) 

0 . 7 

REVIEW  ANOTHER  SOURCE/SINK  TERM  (ANS:  YES  OR  NO)?  n 
CONTINUE  FOR  THIS  TIME  PERIOD  (ANS:  YES  OR  NO)?  n 
REVIEW  ANOTHER  TIME  PERIOD  (ANS:  YES  OR  NO)  ?  y 
WHICH  TIME  PERIOD  DO  YOU  WANT  TO  REVIEW  ?  2 
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INDICATE  WHICH  OF  TH2  FOLLOWING  PARAMETERS  YOU 

WANT  TO  REVIEW  (RESPOND  WITH  THE  PROPER  NUMBER ) 
’  )  HYDRAULIC  PARAMETERS 

2)  INITIAL  CONDITIONS  IN  TRIBUTARIES 

3)  INITIAL  COM  D I T IONS  IN  EFFLUENTS 

4 )  KINETIC  parameters 

5)  distributed  source/sink  parameters 

6)  BIOLOGICAL  parameters 


INDICATE  THE  INDEX  NUMBER  OF  THE  DISTRIBUTED 
SOURCE/SINK  TERM  YOU  WISH  TO  REVIEW 


1 ) 

SL 

2) 

SMH  3 

3) 

SN0  3 

4  ) 

S02 

5) 

SP04 

b) 

SSOD 

7) 

sss 

b) 

ST 

9) 

SAT  1 

10) 

SAT2 

11  ) 

SAT  3 

12) 

SAT4 

13) 

SAT5 

THE  VALUES  CURRENTLT  SPECIFIED  FOR  SSOD  WILL 

BE  LISTED  BELOW  81  REACH.  TO  CHANGE  A  VALUE 
RESPOND  TO  THE  TRAILING  WITH  THE  NEW 

VALUE. 


1  ) 

0.7 

2) 

0.7 

3) 

0.7 

4) 

0.7 

.6 

5) 

0 . 7 

~T 

b) 

0.7 

.  6 

7) 

0.  7 

TT 

B) 

0 . 7 

r-i 

REVIEW  ANOTHER  SOURCE/SINK  TERM  (ANS:  YES  OR  NO)?_£ 
CONTINUE  FOR  THIS  TIME  PERIOD  (ANS:  XES  OR  NO)?_n_ 
REVIEW  ANOTHER  TIME  PERIOD  (ANS:  YES  OR  NO)  ?  n 


YOUR  DATA  SET  IS  STORED  IN  THE  FILE  'TAPE9'. 
REMEMBER  TO  SAVE  OR  REPLACE  IT  IF  YOU  WANT  TO 
USE  IT  IN  A  LATER  SESSION. 

NOTE.  I  f 

/replace . tape9 
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Part  3 


RTV  Test  Using  Data  File  Created  and  Revised 


SIMULA  HONS  HA/E  BEEN  STARTED 
JSI.'IG  THE  FILE  *  T  A  PE  9  ’  AS 
CONTROL  INPUT. 

SIMULATED  WATER  QUALITY  PROFILES  HAVE  BEEN 
OUTPUT  TO  *  TA  PE  3  3 "  FOR  2  TIME 
PERIOD(S)  FOR  THE  FOLLOWING  WATER 
QUALITY  ATTRIBUTES: 


1 ) 

TEMP. 

2) 

BO  05 

3) 

TS3 

“*) 

NH  3 

5) 

N02 

6) 

NO  3 

7) 

P0« 

i ) 

D.O. 

9) 

HARD  . 

10) 

PH 

1  1  ) 

ZN 

12) 

CL2 

13) 

COLI 

WATER  QUALITY  SIMULATIONS  COMPLETE. 

YOU  MAI  NOW  PROCEED  TO  STV  TESTING. 
REMEMBER  TO  SAVE  OR  REPLACE  ’TA PE 3 3  ’  IF 
YOU  PLAN  TO  USE  IT  IN  LATER  SESSONS. 

NOTE.  !  ! 

/replace, tape33 

/-wq rt  v 


LOADER  INFORMATION. 

MAP  OPTIONS  =  OFF 
GLOBAL  LIBRARY  SET  IS  - 
GCSALPH 


♦  THIS  RTV  ROUTINE  TESTS  FOR  VIOLATIONS  OF  AMBIENT  ♦ 

♦  WATER  QUALITY  STANDARDS  AND  QUANTIFIES  THE  ♦ 

♦  SPACIAL  EXTENT  OF  THESE  VIOLATIONS  ♦ 

♦!*s*:+!+:*:*:t:*:t!*s*J*!*!*:*:*:*!ti***tntt*:*i*l*:* 


EXAMPLE  NO.  1 


I C  A  T  £  WHICH  WATER  QUALITY  ATTR I 8UTE ( S )  ARE 
TO  BE  ANALYZED. 


1  ) 

TEMP  . 

5) 

N02 

2) 

3005 

b) 

NO  3 

3) 

TSS 

7) 

PQ4 

N  ) 

NH  3 

0) 

D.O 

•  9) 

HARD  . 

10)  PH 

11)  ZN 

12)  CL2 

13)  COL  I . 

RESPOND  WITH  THE  TOTAL  NUMBER  OF  ATTRIBUTES 

FOLLOWED  BY  THE  APPROPRIATE  INDEX  NUMBERS 
?5,2 .3.M.3 , ’ 1 
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Section  A:  Violation  Test 


INPUT  LOCAL  AMBIENT  WATER  QUALITY  STANDARDS 


-ME 

PERIOD 

NO .  1 

BOD5 

UPPER 

LEVEL 

STANDARD 

? 

i£ 

LOWER 

LEVEL 

STANDARD 

7 

TSS 

UPPER 

LEVEL 

STANDARD 

7 

12 

LOWER 

LEVEL 

STANDARD 

? 

N  H  3 

UPPER 

LEVEL 

STANDARD 

? 

1 . 5 

LOWER 

LEVEL 

STANDARD 

7 

D.O  . 

UPPER 

LEVEL 

STANDARD 

7 

LOWER 

LEVEL 

STANDARD 

7 

ZN 

UPPER 

LEVEL 

STANDARD 

7 

1 . 0 

LOWER 

LEVEL 

STANDARD 

7 

STANDARDS  CONSTANT  OVER  TIME  (X  OR  N)  ?  y 


REPORT  ON  WATER  QUALITY  VIOLATIONS 


1)  TOTAL  RIVER  MILES  IN  VIOLATION  OF  STANDARDS. 

TIME  PERIOD 

ATTRIBUTE  123456 


B305 

33.5 

13-5 

TSS 

33.5 

33-5 

NH3 

16.7 

10.3 

D.O. 

0.0 

0.0 

ZN 

3.2 

0.0 

DO  YOU  WANT  A  GRAPHICAL  OUTPUT  FOR  ANY  ATTRIBUTE 
PROFILES  CANS:  YES  OR  NO)?  y 
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ATTRIBUTE  NUMBER  ?  2_ 
TIME  PERIOD  ?  1 


SOD 5  :  TIHE  PERIOD  NO. 


1  20 

C110 
0  100 
N  90 
C  80 
.  70 

(  60 
M  50 
0  40 
/  30 
L  20 
) 

10 


00 

0  00 


00 

000  0  00 
0000  000 


00000 


00 

0 

0 

0 

0 

0 

000000000000000 


0  000000 **♦*♦♦♦*♦*♦♦♦♦♦*♦♦■*♦♦♦♦♦♦■»***♦•*♦♦ ♦♦♦♦♦♦♦*♦*♦♦*♦♦♦♦**♦+* 
0  4  a  12  16  20  24  28  32  36  40 

2  6  10  14  18  22  26  30  34  38 

DISTANCE  DOWNSTREAM  (MILES) 

PLOT  8005  FOR  ANOTHER  TIME  PERIOD  ?  _y_ 

TIME  PERIOD  ?  2 

B0D5  :  TIME  PERIOD  NO.  2 


000000 


55  *- 
♦  ♦ 

C50  ♦ 

045  ♦♦ 

N  40  ♦♦ 

C  *♦ 

.35  * 

30  *♦ 

(25  ♦♦ 

M  ♦♦ 

020  ♦ 

/ 1  5  ♦♦ 

L’O  . s 

)  ♦.  0  0 
5  ♦  00  0  0 

0  00000000000 ♦♦♦♦■.♦♦♦♦♦♦♦♦»» ♦♦♦*00000000000000000000000000000000 
0  4  8  12  16  20  24  28  32  36  40 

2  6  10  14  id  22  26  30  3“)  38 

DISTANCE  DOWNSTREAM  (MILES) 


0 

0000 
0  0 
0  0 
0 
0 
0 
0 
0 
0 
0 


000 
0  0 
0  0 
0 
0 
0 
0 
0 


00000000 

00 


47 


PLOT  3005  FOR  ANOTHER  TIME  PERIOD  ? 

PLOT  ANOTHER  HATER  QUALITY  ATTRIBUTE  ?  y_ 

ATTRIBUTE  NUMBER  ?  3_ 

TIMS  PERIOD  ?  1 

“  TSS  :  TIME  PERIOD  NO.  1 


100  ♦ 

95  ♦  0000 

C  90  *  0  0000000000 

0  85  t  0  0 

S  75  *  0000  0  0 

C  70  ♦  0  000000  0 

.  65  ♦  0  0 

55  ♦  0  0 

(  50  +  0  0 

H  40  «  0  00000 

o  35  *  o  00000000000000000000000000 

/  30  ♦  0 

L  20  ♦  0 

)  15  . S 

5  ♦  000000 

Q  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦♦ 

o  4  a  12  16  20  24  28  32  36  40 

2  6  10  14  1 8  22  26  30  34  38 

DISTANCE  DOWNSTREAM  (MILES) 

PLOT  TSS  FOR  ANOTHER  TIME  PERIOD  ?  y 

TIME  PERIOD  7  _2_ 

TSS  TIME  PERIOD  NO.  2 


50.  ♦ 

47.5  *  0000 

C  45.  t  00 

042.5  ♦  0000  0  0 

N37.5*  00  00 

C  35.  ♦  0  000000  0 

.32.5  ♦  0  0000000000 

27.5  *  0  0 

(225.  ♦  0  000000000000000000000000000000 

M  20.  ♦  0 

017.5  *  0 

/I  15.  ♦  0 

L  10 . S 

)  7.5  ♦ 

2.5  ♦ 

0  4  8  12  16  20  24  28  32  36  40 

2  6  10  14  18  22  26  30  34  38 

DISTANCE  DOWNSTREAM  (MILE? 
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PLOT  TSS  FOR  ANOTHER  TIME  PERIOD  ?  n 

plot  another  water  gualiti  attribute  ?x~ 

ATTRI3UTS  NUMBER  ?  jj_ 

TIME  PERIOD  ?  1 

NH 3  :  TIME  PERIOD  NO.  1 


IN. 

C  *  3 . 
012. 
Nil. 
CIO. 
.  9. 
8  . 
(  7. 
M  6. 
0  5. 
/  4. 

L  3. 
)  2. 
1  . 
0. 


0 

0 

00 
-  0 
0  0 

0  0 

0  00 
0  00 
0  00 

0  0 
0  000 
0  00000 


. s 

*00  0 
*♦♦*♦000  000  ♦♦♦♦♦♦♦♦♦♦♦*♦**♦♦**00  00  00000000000  0  000000000000000 
0  4  8  12  16  20  24  28  32  36  40 

2  6  10  14  18  22  26  30  34  38 

DISTANCE  DOWNSTREAM  (MILES) 


PLOT  NH3  FOR  ANOTHER  TIME  PERIOD  ? 

TIME  PERIOD  ?  2 

~  NH3  :  TIME  PERIOD  NO.  2 


6.5 

C  6. 
05.5 
N  5. 
C  4 . 5 
.  4 . 

3-5 
(  3. 
M2 . 5 
G  2. 
/ 

61.5 
)  1. 
.5 
0  . 


0 

00 

0  00 
0  0 
0  00 


♦  0  000 
♦  0  00 
♦  0  00 

♦  00 

*00 
0000  0  0 


. s 

♦00  ooooooooo 

♦  00  0 
♦♦♦♦♦000000*******************0000000000000000000000000000000 
0  4  8  12  16  20  24  28  32  36  40 

2  6  10  14  18  22  26  30  34  38 

.  STANCE  DOWNSTREAM  (MILES) 
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PLOT  NH3  FOR  ANOTHER  TIME  PERIOD  ?  n 

PLOT  ANOTHER  WATER  OUALITT  ATTRIBUTE  ?  /  “ 

ATTRIBUTE  NUMBER  ? 

TIME  PERIOD  ?  1 

”  D.O.  :  TIME  PERIOD  NO.  1 


a. 5  * 

8.  +00000 

CT.5  +0  OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOQOOOOQOOOO 
0  7.  +0  0  0 

N  6 . 5  +0 


C  o . 

.5.5  +0 
4.5  0 
(  4.  0 

M3  -5  0 
0  3-  0 
/2.5  0 
LI  .  5  0 
)  1  ■  ♦ 

.5  ♦ 

0 .  ♦♦♦♦ 
0 

2 


S 


6  10  14  18  22  26  30  34  38 

DISTANCE  DOWNSTREAM  (MILES) 


0.0. 


FOR  ANOTHER  TIME  PER  CO  ?  v 


TIME  PERIOD  ?  _2 

D.O.  :  TIME  PERIOD  NO.  2 


9.5 
9. 

C8.5 
-1  3  . 
M  1' .  5 
C6.5 
.  6. 

5.5 
(4.5 

M  4, 

03-5 
/2.5 
L  2. 
>1.5 
.5 
0. 


♦0000000000 

♦0  00  0  000000 

-0  0  ooooocoo  oooo-' c'"’"  "^mooooooooc"'" 

^  GOO 

0 

. . 


2 


8  ’2  16  20  24  28  32  36  40 

6  10  14  13  22  26  30  34  38 

DISTANCE  DOWNSTREAM  (MILES) 
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PLOT  D.O.  FOR  ANOTHER  T  IMS  PERIOD  7  fl 

J,_jr  a .'io Mia  : : 7 ; 7  j r -r  ?  j_  “ 

ATTRIBUTE  NUMBER  7  ’ 

TIME  PERIOD  ? 

zn  time  period  wo.  1 


i 

\ 


\ 


I 

1 


2.  ♦ 

1-9  000000 
C  1  .  8  *  0 

01.7*  0 

N1 .5  *  0 

Cl  .  »  ♦  0 

.  1  .  3  *  0 

1.1.*  0 


(  I . s 

M  .  8  *  0 

G  .7  ♦  0 

/  .6  *  0  0000000000000 

L  .4  ♦  000000  000000  0 

)  -3  ♦  0000  0 

.  ♦  0000000000000000000000000000000 
0. 

0  4  8  12  16  20  24  28  32  36  40 

2  6  10  14  ’8  22  26  30  34  38 

DISTANCE  DOWNSTREAM  (MILES) 

plot  ZN  FOR  ANOTHER  TIME  PERIOD  7  y 


TIME  PERIOD  ?  2_ 

ZN  :  TIME  PERIOD  NO.  2 


| 

i 


n 

i 


1  . 

C  .  9 
0 . 8 
N 

C.7 

.  .6 


0. 


000000 
*♦  0 
*♦  0 

0  0000 
0  0000000  0 
000000000  000000000 

0 

^^.^^^^^♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦+00000000000000000000000000000000 

4  8  12  16  20  24  28  32  36  40 

2  6  10  14  18  22  26  30  3»  38 

DISTANCE  DOWNSTREAM  (MILES) 


PLOT  ZN  FOR  ANOTHER  TIME  PERIOD  7  _n 

PLOT  ANOTHER  WATER  QUALITI  ATTRIBUTE  7  n_ 


THIS  CONCLUDES  ’WORTV”.  TOU  MAI  EXECUTE  MORE 
R  TV  ROUTINES  NOW.  BEGIN  A  MITIGATION 
LOOP  OR  SIJNOFF. 

iiliitiaiiitliiiitltlliiiiltltiMiiiiitttiitlii 

2.993  CP  SECONDS  EXECUTION  TIME 
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Section  B:  Sapvobio  Index  Analysis 


\ 

/-Si Pt  V 


SAPBOBIC  INDEX  ANALYSIS 

Foa 

EXAMPLE  NO.  1 


WATSB  QUALITY  B I V  ER  MILES  IN  TIME  PERIOD 

DESIGNATION  12  3^56 


PUREST  WATER 

0.0 

0 . 0 

CLEAN  WATER 

3-2 

26.5 

MILD  POLLUTION 

2.5 

0.0 

POLLUTED 

6.8 

0.0 

HEAVILY  POLLUTED 

’  «  .  0 

in  J 

RAW  SEWAGE 

13.5 

2.7 

SEPTIC  CONDITION 

0.0 

0.0 

DO  YOU  WANT  FURTHER  QUANTIFICATION  OF  THIS 
(ANS:  YES  OR  NO)  7  y 


INPUT  TIME  PERIOD  OF  INTEREST  ?  1 

THIS  SECTION  ISN'T  OPERATIONAL  YET,  BUT  THE 
OUTPUT  WILL  BE  LOCATIONS  OF  ZONES  IN  EACH 
WATER  QUALITY  DESIGNATION  FOR  THE  SPECIFIED 
TIME  PERIOD. 

.038  CP  SECONDS  EXECUTION  TIME 
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thu'tion  (' 


Itov  ironmcn  La  l  Toxicrity 


/ -t  urt  v 

LOADER  I N  F  0  R  M  A  T 1 0 .1 . 

MAP  OPTIONS  =  OFF 
GLOBAL  LIBRARY  SET  IS  - 
GCSALPH 


THIS  RT V  ROUTINE  TESTS  ENVIRONMENTAL  TOXICITY 


EXAMPLE 

NO.  " 

..„ICAT 

S  WHICH  WATER 
TO  BE  ANALYZED 

QUALITY  ATTRI3UTECS)  ARE 
FOR  THEIR  TOXIC  EFFECTS 

1  ) 

TEMP. 

5) 

N02 

2) 

BOD5 

6) 

N03 

3) 

TSS 

7) 

PO  4 

A) 

NH3 

8) 

D.O. 

9) 

HARD. 

10) 

PH 

1  1  ) 

ZN 

7  2) 

CL2 

13) 

COLI 

RESPOND  WITH  THE 
FOLLOWED  BY 
7  3,«,H.12 


TOTAL  NUMBER  OF 
THE  APPROPRIATE 


TOXICANTS 
INDEX  NUMBERS 


SPECIFY  TARGET  SPECIES: 

DESIGNATE  STREAM  TYPE  (WsWARM  WATER,C*COLD  WATER)  7 _w 

REPRESENTATIVE  SPECIES  LIST: 

1  )  FATHEAD  MINNOW 

2)  CARP 

3)  BLUEGILL 

«)  CHANNEL  CAT 
5)  LARGEMOUTH  BASS 

RESPOND  WITH  NUMBER  OF  TARGET  SPECIES  DESIRED  AND 
WITH  THEIR  APPROPRIATE  INDEX  NUMBER ( S ) • 

?  7  ,  3 

INPUT  THE  96  HOUR  LC50 ’S  FOR  THE  FOLLOWING  SPECIES 
AND  POTENTIONAL  TOXICANTS: 


BLUEGILL 

NH3 

ZN 

CL2 


7  .6 

7  J~ 

7  .2 


SPECIFY  TIME  PERIOD  OF  INTEREST  7  1 
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MAXIMUM  AND  MEAN  (IN  PARENTHESES)  TOXICITY  UNITS 


TARGET  SPECIES 
BLUEGILL 


TOXICANT 

TOTAL  N  H  3  ZN  CL2 

28.136  22.238  .250  7. "77 

(  4.272)  (  <4. 272  )  (  .  054  )  (  2.903  )  ( 


DO  XOU  WANT  A  GRAPHICAL  OUTPUT  FOR  TOXICITX  UNITS 
VS.  LOCATION  DOWNSTREAM  (ANS:  YES  OR  NO)  ?  y 


INPUT  TOXICANT  NUMBER  ? 

1  )  NH  3 

2)  ZN 

3)  CL2 

INPUT  TARGET  SPECIES  INDEX  7  ' 


BLUEGILL  ;  NH  3  I  TIME  PERIOD  1 

( T  =  TOTAL  T.U.'S;  ‘sT.U.'S  FROM  SPECIFIED  ATT.) 


30 

28 

T 

26 

TT 

T24 

TT 

022 

T  TT 

X20 

T  ‘  TT 

.  18 

T*  *  TT 

16 

T  •  TT 

ui  4 

T. 

T  ••  T 

N  1  2 

.  7  ? 

7  t  •  7 

fff 

I  <J 

. . TT  TTT 

r 

TTTTTTT 

T  6 

. .  TT  TTTTT  •• 

•••  TTT 

S  4 

• 

* 

2 

0 

• 

• 

• ‘TTTTT TTTTTTT TTTT TTTTT TTTTTTT TTT TTT 

0  4 

8  12 

16  20  24  28  32  36  40 

2  6  10  14  18  22  26  30  34  38 

DISTANCE  DOWNSTREAM  (MILES) 


PLOT  TOXICITY  IMPACTS  FROM  NH  3 

FOR  ANOTHER  TARGET  SPECIES  7  _n 
PLOT  FOR  ANOTHER  TOXICANT  7  y 
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INPUT  TOXICANT  NUMBER  ?  3. 

INPUT  TARGET  SPECIES  INDEX  ?  1 


3LUSGILL  ;  EM  ;  TIME  PERIOD  1 

( T  =  TOT  AL  T.U.'S;  •*T.U.'S  FROM  SPECIFIED  ATT.) 


30 

28 

T 

26 

TT 

T2  4 

TT 

022 

T  TT 

X20 

T  TT 

.  18 

T  TT 

1  0 

T  TT 

U  1  4 

T. 

T  T 

M  1  2 

.  TT 

T  TT 

TT 

i  a 

.  .  TT  TTT 

T 

TTTTTTT 

T  6 

TT  TTTTT 

TTT 

S  4 

T 

2 

TTTTT TTTTTTTTT TTTTTTT TTTTTTTTTT TTT 

0  . . »••»•••••••■»<•••••••• . * . lit . . 

OR  a  12  16  20  24  28  32  36  40 

2  6  10  14  18  22  26  30  34  38 

DISTANCE  DOWNSTREAM  (MILES) 


PLOT  TOXICITX  IMPACTS  FROM  ZN 

FOR  ANOTHER  TARGET  SPECIES  ?  _n 
PLOT  FOR  ANOTHER  TOXICANT  ?  y 

INPUT  TOXICANT  NUMBER  ? 

INPUT  TA AGET  SPECIES  INDEX  ?  1 


BLUEGILL  ;  CL2  ;  TIME  PERIOD  1 


(T=TOTAL  t 

.U.'S;  »«T.O.'S  FROM  SPECIFIED  ATT.) 

30 

28 

T 

26 

TT 

T2  4 

TT 

022 

T  TT 

X20 

T  TT 

.  18 

r 

TT 

16 

T 

TT 

U1  4 

T. 

T 

T 

N  1  2 

.  TT 

T 

TTTT 

I  8 

.  .TT 

TTT  T 

TTTTTTT 

T  6 

TT*TTTTT# * • • 

•••••  TTT 

S  4 

« 

2 

• 

TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT 

0 

4  8 

12  16  20  24  28  32  36  40 

2 

6 

10  14  18  22  26  30  34  38 

DISTANCE  DOWNSTREAM  (MILES) 
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PLOT  TOXICITY  IMPACTS  FROM  CL2 

FOS  ANOTHER  TARGET  SPECIES  ?  n 
PLOT  FOR  ANOTHER  TOXICANT  ?  n  '' 


DO  YOU  WISH  TO  CONTINUE  ANALYSES  FOR  ANOTHER  TIME 
PERIOD  CANS:  Y  OR  N)  ?  y 

SPECIFY  TIMS  PERIOD  OF  INTEREST  ?  2 


REPORT  ON  TOXICITY  IMPACTS  IN  TIME  PERIOD  2 


MAXIMUM  AND  MEAN  (IN  PARENTHESES)  TOXICITY  UNITS 


TARGET  SPECIES 
BLUEGILL 


TOXICANT 
TOTAL  NH3 

HI. 276  10.394 

I  -.114)  (  2.114)  ( 


ZN 

.060 

.015) 


CL2 
2 . 727 
.826)  ( 


DO  YOU  WANT  A  GRAPHICAL  OUTPUT  FOR 
VS.  LOCATION  DOWNSTREAM  (ANS: 


TOXICITY  UNITS 
YES  OR  NO)  ?  y 


INPUT  TOXICANT  NUMBER  V 

INPUT  TARGET  SPECIES  INDEX  ?  1 


BLUEGILL 
( T  =  TOTAL  T . U .  'S  ; 


;  NH3  ;  TIME  PERIOD  2 
•sT.U.'S  FROM  SPECIFIED  ATT.) 


42.5 
40. 

37-5 
T  35. 
032-5 
X27  •  5 
.  25. 

22.5 
U  1  20 . 
N  15. 
112.5 
T  10. 
S  7.  5 

2-5 
0  . 


T 

TTT 
T  TT 

.  T  T 

T  TT 
T  TTT 

T  TT 

T  TTTT 

TT  T  T 

.  TT  T  TTTT 

TTTTTTT • * • *  TTTTTT 

•  •••••#•••«  T 

•  ••in  •  •«•«««•  r 

. *••••• . ••TTTTTTTTTTTTTTTTrTTTTTTTTTT  TTTTTT 

04  8  12  16  20  24  28  32  36  40 

2  6  10  14  18  22  26  30  34  38 


DISTANCE  DOWNSTREAM  (MILES) 


PLOT  TOXICITY  IMPACTS  FROM  NH  3 

FOR  ANOTHER  TARGET  SPECIES  ?  n 
PLOT  FOR  ANOTHER  TOXICANT  7  y 


I  :i I*  I  r  i'OX  10ANT  NUMOKA  ?  ^ 

I  N  P  U  T  TARGET  SPECIES  INDEX  7  1 


BLUEGILL  ;  Z.N  ;  TIME  PERIOD  2 

(TsTOTAL  T.U.'S;  •sT.U.'S  FROM  SPECIFIED  ATT.) 


42.5 
40  . 

37.5 
T  35- 
032.5 
X27.5 
.  25. 

22.5 
U120. 
N  15. 
112.5 
T  10. 
S  7.5 

2.5 
0  . 


T 

TTT 
T  TT 


TT 


TT 


T 
T 
T 
T 
T 
T 
T 

TTTTTTT 


TT 

TTT 


TT 

TTTT 

T 

TTTT 


TTTTTT 

T 

T 

••••••TtTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT 


12  16  20  24  28 

10  14  ia  22  26  30 

DISTANCE  DOWNSTREAM  (MILES) 


32  36  40 

34  38 


PLOT  TOXICITY  IMPACTS  FROM  ZN 

FOR  ANOTHER  TARGET  SPECIES  ?  n 
PLOT  FOR  ANOTHER  TOXICANT  ?  y  *” 

INPUT  TOXICANT  NUM9SR  7 

INPUT  TARGET  SPECIES  INDEX  ?  1 


BLUEGILL  ;  CL2  ;  TIME  PERIOD  2 

(TrTOTAL  T.U.'S;  »=T.U.'S  FROM  SPECIFIED  ATT.) 


42.  5  - 

40 .  .  T 

37.5  .  TTT 

T  35.  -  T  TT 

032.5  .  T  T 

X27.5  .  T  T T 

.  25.  .  T  TTT 

22.5  .  T  TT 

U 1 20 .  .  T  TTTT 

N  15.  TT  T  T 

112.5  •  TT  T  TTTT 

T  10.  .  TTTTTTT  TTTTTT 

S  7.5  .  T 

2.5  .  . .  T 

0.  •••••• . •••••••••TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT 

o  4  8  12  16  20  24  26  }2  36  v0 

2  6  10  14  18  22  26  30  34  38 

DISTANCE  DOWNSTREAM  (MILES) 


PLOT  TOXICITY  IMPACTS  FROM  CL2 

FOR  ANOTHER-  TARGET  SPECIES  ?  n 
PLOT  FOR  ANOTHER  TOXICANT  7  n 
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DO  YOU  WISH  TO  CONTINUE  ANALYSES  FOR  ANOTHER  TIME 
PERIOD  (ANS:  Y  OR  N)  ?  n 


THIS  CONCLUDES  'TURTV YOU  MAY  EXECUTE  MORE 
RTV  ROUTINES  NOW,  BEGIN  A  MITIGATION 
LOOP  OR  SIGNOFF. 

2.004  CP  SECONDS  EXECUTION  TIME 


aye 

3KVN9RS  COSTS:  690.570  SRUS  AT  $.0068  »  $4.70 
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APPENDIX  C: 

RIAS  SOURCE  PROGRAMS 


type , rl as/g 


NOTE . ! 

a  a  * 

a  a  a 

NOTE . ! 

a 

NOTE . ! 

• 

WELCOME  TO  THE 

NOTE. ! 

• 

RATIONAL  IMPACT  EVALUATION  SYSTEM 

NOTE. 1 

a 

NOTE. ! 

a 

U.S.A.C.E. 

NOTE. i 

a 

CONSTRUCTION  ENGINEERING  RESEARCH 

NOTE. ! 

a 

LABORATORY 

NOTE. 1 

a 

NOTE . ! 

a 

VERSION  1 . 2 

A  J  r  £  .  ! 

a 

NOVEMBER  ’979 

NOTE . 1 

a 

NOTE . ! 

a  a  a 

a  a  a 

NOTE . ! 1 

NOTE. ! t 

NOTE. ! ! 

NOTE. ! 

THIS  COMPUTER  BASED  IMPACT  ASSESSMENT  SYSTEM 

NOTE. ! 

CAN  BE  ACTIVATED  3Y  EXECUTING  A  SERIES  OF 

NOTE. ! 

PROCEDURE  FILES  WHICH  ARE  DESCRIBED  BELOW: 

NOTE. ! ! 

NOTE. ! 

i  > 

-FQUES 

NOTE. ! 1 

NOTE. ; 

THIS  PROC  FILE  EXECUTES  A  PROGRAM  WHICH  ASKES 

NOTE. ! 

A  SERIES  OF  QUESTIONS  ABOUT  THE  ENVIRONMENTAL 

NOTE. 1 

SETTING  FOR  A  SPECIFIC  PROJECT  AND  ORGANIZES 

NOTE. ! 

THE  INFORMATION  OBTAINED  INTO  A  DATA  FILE 

NOTE. ! 

WHICH  IS  USED  AS  INPUT  FOR  .IMULATIONS  AND 

NOTE. ! 

OTHER  EVALUTION  PROTOCOLS.  FQUES  ’  CAN  BE 

NOTE. ! 

USED  TO  SETUP  A  NEW  DATA  FILE  OR  TO  REVISE 

NOTE . ! 

AN  EXISTING  DATA  FILE.  OUTPUT  FROM  'FQUES' 

NOTE. ! 

IS  TO  A  MASS  STORAGE  FILE  CALLED  ' T APE9  * . 

NOTE. 1 ! 

NOTE. 1 

2) 

-SIMWQ 

NOTE . ! ! 

NOTE . ! 

THIS  PROC  FILE  EXECUTES  A  WATER  QUALITY 

NOTE  .  1 

SIMULATION  MODEL  WHICH  PREDICTS  PRIMARY 

NOTE . ! 

PHYSICAL  AND  CHEMICAL  IMPACTS  IN  THE  AQUATIC 

NOTE. ! 

ENVIRONMENT  BASED  OF  THE  ENVIRONMENTAL 

NOTE. 1 

SETTING  AND  PROJECT  SPECIFICATION  INFORMA- 

NOTE .  ! 

TION  STORED  IN  THE  OUTPUT  FROM  'FQUES'. 

NOTE.  ! 

OUTPUT  FROM  'SIMWQ'  IS  TO  A  MASS  STORAGE 

NOTE.  ! 

FILE  CALLED  '  TA  P  E  3  3  '  • 

NOTE. ! ! 

NOTE.  1 

3) 

-RTVS 

NOTE.  1  1 

NOTE. 1 

THIS  PROC  CALLS  A  LISTING  OF  A  SERIES  OF 

NOTE. ! 

PROGRAM 3  WHICH  QUANTIFY  ENVIRONMENTAL  IMPACTS 

NOTE . ! 

IN  THE  AQUATIC  ECOSYSTEM.  THE  UJF.R  CAN 

NOTE . 1 

THEN  CHOOSE  WHICH  RTV  MODELS  ARE  APPLICABLE 

NOTE.  ! 

AND  CAN  EXECUTE  EACH  WITH  THE  APPROPRIATE 

NOTE.  ! 

PROC  CALL.  ALL  RTV  MODELS  REQUIRE  INPUT 

NOTE.  1 

FROM  ' T  A  PE  9  '  AND/OR  '  TA  PE  3  3  ‘  • 

NOTE,  t  ! 

NOTE.  ! 

a  a  a  i 

NOTE.  !  t 
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type  , 

- t'quea / g 

GET,  BFQUES/IDORI  AA756P  . 

RWF  . 

B FQUES 

NOTE. ! 

j 

NOTE. ! 

i 

NOTE.  1 

! 

NOTE.  ! 

! 

NOTE.  ! 

YOUR  DATA  SET  IS 

NOTE. ! 

REMEMBER  TO  SAVE 

NOTE. ! 

USE  IT  IN  A  LATER 

NOTE. ! 
/ 

r 

STORED  IN  THE  FILE  ' TA  PE  9  *  . 
OR  REPLACE  IT  IF  YOU  WANT  TO 
SESSION . 


type , slrawq/g 

NOTE.  |  1 

NOTE  .  !  ! 

NOTE. 1 1 

NOTE . 1 

SIMULATIONS  HAVE  BEEN  STARTED 

NOTE.  ! 

USING  THE  FILE  ' TA  PE  9  '  AS 

NOTE.  ! 

CONTROL  INPUT. 

GET, 33 IMWQ/ I  0=39’  447562. 


RWF. 
BSIMWQ . 
NOTE . ! ! 
NOTE .  !  ! 
NOTE  .  !  ! 
NOTE  .  !  ! 
NOTE, ! 
NOTE.  1 
NOTE. ! 
NOTE. ! 
NOTE. I ! 
/ 


WATER  QUALITY  SIMULATIONS  COMPLETE. 

YOU  MAY  NOW  PROCEED  TO  R TV  TESTING. 
REMEMBER  TO  SAVE  OR  REPLACE  'TA  PE  3  3  *  IF 
YOU  PLAN  TO  USE  IT  IN  LATER  SESSONS . 


type , rt va/g 

M0TE. !  THE  FOLLOWING  RATIONAL  THRESHOLD  EVALUATIONS 

NOTE. !  MAY  BE  INITIATED  AT  THIS  TIME  BY  EXECUTION  OF 

NOTE. I  OF  THE  APPROPRIATE  'PROC  FILES: 

NOTE. ! ! 


NOTE . ! 

-WQRTV 

NOTE.  1 

ASSESSMENT  OF 

OVERALL 

WATER  QUALITY 

IMPACTS  . 

NOTE.  ! 

-3IRTV 

NOTE.  1 

ASSESSMENT  OF 

IMPACTS 

OF  ORGANIC  ENR 

ICHMEN  T . 

NOTE. ! 

-TURTV 

NOTE. 1 

ASSESSMENT  OF 

ENVIRONMENTAL  TOXICITY 

IMPACTS 

NOTE.  1  ! 
NOTE.  1 

REMEMBER  TO  HAVE  TAPE9'  AND  ' 

’  TA  PE  3  3  ’ 

AVAILABLE 

NOTE.  I  AS  LOCAL  FILES  BEFORE  EXECUTION. 

NOTE .  I  I 
/ 
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type , turt v/ g 


GET,BTURrV/ID«31,11t‘t7562. 
RWF  . 

GRAB . GC3ALPH/F . 

Btuarv . 


ytpe  •BE- • 
t ype , at  urt  v/g 


t  ype , wqrt  v/g 

GET, BWQR TV/ IDs  3  Ml  4  47  562. 
GRAB , GC3ALPH/F . 

RWF. 

BWQRTV . 

/ 


type , slrtv/g 


GET,BSIRTV/IDs341447562. 
RWF  . 

BSIRTV. 

/ 


SOU..CE  LISTING  FOR  FQUeJS 


PROGRAM  FILTSRQ( INPUT, OUTPUT , T A P E 1 : I N PUT , TA PE2=OUTPUT , TAPE9 ) 
COMMON /NAMES /TITLE t  d ) , CN AME ( 20 ) , SNAME ( 20 ) , KNAMP ( 20 ) 

COMMON /CNTRL /NR , NTP , NWQC ,NB,NIT,NI£,IT(15),IE(15),HC0DE, 

♦  MCODE ( 8 ) , NCWQC , NNCWQC , NTWOC , NS , NK 

COMMON / 1 WQC/TWQ( 15, 20) ,  £WQ(  1 5 , 20  )  , TO ( 1 5 ) , EQ ( 1 5  > 

COMMON /KDATA/K20(  20  , 20) 

COMMON /S DATA/S20C 20, 20) 

COMMON /MDATA/LR (20 ) , 0A( 20 ) , DEPTH £20) , VEL(20) , BCODE( 15). 

♦  H A  V ( 20 )  , HB  V ( 2  0 ) , H  A  D ( 20) ,HBD(20) 

C0MM0N/8DATA/A ( 20 ) 

INTEGER  TITLE.CN AME. HCODE, KNAME. 'ME 
INTEGER  DINDEXl  125) ,CNTRLS(8) 

REAL  LR.K20 

DATA  ( C  N  A  M  E ( I ),I*1,3)/"TSMP.","BOD5". "TSS" , "NH3","N02","N03" 

♦  "P04" , "D.O. "/ 

DATA  (  KNAME  ( I  )  ,  1=  1 , 20  )  /  '•  KANH  3  "  ,  "  K  ANO  3  *  , “KAPOA" , "KL" ,  ”KLS" , 

♦  "KNDET" , "KNH3" , "KN02" , "KN03" , "KPDET" , "KP04" , "KR" , "KSS" , 

♦  "KSSS  "  ,  "KT"  ,  "KTSS "  ,  "KNCA  1  ”  ,  "KNC  A2  "  .  "  .<  N  C  A  3  "  ,  "KNCAa"  ' 

DATA  (SHAME( I), 1*1,20)/ "SL", " S  N  H  3 " , "SNO 3 " , "S02" , "SPO A " , 

♦  "SSOO" , "SSS" , "ST" , "SAT  1 " , "SAT2" , "SAT  3” , "SATA" , "SAT5" , "SAT  6" 

♦  "SAT7", "SAT8", "SAT9", "SAT10", "SAT1 1" , "SAT12"/ 

PRINT*. "  » 

PRINT*,"  " 

PRINT*,"  " 

PRINT*," 

PRINT*,"  t  #• 

PRINT*,"  t  THIS  PROGRAM  ALLOWS  #" 

PRINT*,"  #  THE  USER  TO  BUILD  UP  A  NEW  #« 

PRINT*,"  #  DATA  FILE  OR  TO  REVISE  AN  OLD  DATA  #" 

PRINT*,"  t  FILE  FOR  SUBSEQUENT  CONTROL  AND  #• 

PRINT*,"  #  INPUT  FOR  WATER  QUALITY  SIM-  P" 

PRINT*,"  P  ULATONS  UNDER  'RIAS'  P" 

PRINT*,"  P  P" 

PRINT*,"  PPFPPPPPPFPPPPPPPPPFPPPPPPPPPFPPPFPPPPPPP" 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  " 

CALL  0PENMS(9,DINDEX, 125,0) 

PRINT*,"  DO  YOU  WISH  TO  BEGIN  CREATING  A  NEW  DATA  FILE" 

PRINT*,"  (ANS:  YES  OR  NO)  ", 

READt 1,901)  I ANS 1 
IF ( EOF ( 1 ) )  10,11 

IFdANSI  ,NE,  »T«.  AND,  IANS1  ,NE.  "N"  )  GO  TO  10 
IFUANS1  ,EQ."N")  GO  TO  500 
PRINT*,"  " 

PRINT*," 

PRINT*,"  TIME  INVARIANT  PARAMETERS” 

PRINT*," 

PRINT*,"  " 

PRINT*,"  " 

PRINT*."  I)  TYPE  IN  THE  NAME  OF  THE  DATA  SET" 


TIME  INVARIANT  PARAMETERS" 

szssazssassszrsstsssssasszssssxzzsn 


I)  TYPE  IN  THE  NAME  OF  THE  DATA  SET" 


RE A D (  1,911)  (TITLEt I  )  . I»  1  , 8) 
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I.-UOfUJ)  20,2  1 
21  PRIST*,"  " 

PRIST*,"  II)  ANSWER  THE  FOLLOWING  QUESTIONS  ABOUT  SIMULATION 

PRINT*,"  CONTROL  PARAMETERS." 

PRIST*. "  " 

30  PRINT*,"  HOW  MANY  REACHES  (ANS:1-20)  ", 

READ ( 1 , •)  NR 

IF  C  EOF  C  1 ) )  30,31 

31  PRINT*,"  HOW  MANY  TIME  PERIODS  (ANS: 1-12)  ", 

READ! 1 , * )  NTP 

IF  C  EOF  Cl))  31,32 

32  PRINT*, "  HOW  MANY  ADDITIONAL  CONSERVATIVE  WATER  QUALITY" 

PRIST*,"  ATTRIBUTES  BEYOND  THE  CORE  ATTRIBUTES” 

PRIST*,"  (ANS:  0-12)  ", 

R  E A  D (  1 , * )  NCWQC 
IFCEOF(I))  32,33 

33  PRIST*,"  HOW  MANY  ADDITIONAL  SOMCONSER V  AT  I V E  ATTRIBUTES" 

PRIST*,"  (AMS:  0-9)  ”, 

READl 1 , • )  NNCWQC 
IF  (  EOF  (  1  )  )  3  3,  3 ** 

39  NTWQC  =  di-NCWQC»NNCWQC 
NWQC=8 

MK  = 1 6  +  NNCWQC 
NS=NTWQC 

IF ( NTWQC . L£ . 6 )  00  TO  70 

PRIST*,"  INPUT  THE  NAMES  OF  THESE  ADDITIONAL  PARAMETERS. 

IF (NCWQC . LT . 1  )  GO  TO  51 

PRINT*, "  CONSERVATIVE  ATTRIBUTES:" 

DO  50  1=1 .NCWQC 

99  PRINT* ATT.  NO.  " , 1*8 , "  )  • , 

READ(1,913)  CNAME(8*I) 

IF ( EOF (  1  )  )  R9.50 

50  CONTINUE 

IF( NNCWQC .LT. 1 )  GO  TO  70 

51  PRINT*,"  (fONCONSERVATI VE  ATTRIBUTES: " 

DO  60  1=1, NNCWQC 

59  PRIST*,"  ATT.  NO.  " , I*8  +  NCWQC , ” ) " , 

REAO(1,913)  CNAME(8*NCWQC*I) 

If ( EOF (  1  ) )  59.60 

60  CONTINUE 

70  PRINT*,"  INPUT  THE  MODEL  I NC  CODES  FOR  THE  EIGHT  CORE" 

PRINT*,"  ATTRIBUTES  (ANS:  0  OR  1)  " 

PRIST*," 

READ(1,»)  (MCODE(I), 1=1,8) 

IF(E0r(1))  70,71 

71  PRIST*,"  " 

PRINT*,"  III)  INPUT  THE  PARAMETERS  DESCRIBING  THE  WATERSHED” 

PRINT*,"  STRUCTURE  FOR  THIS  SIMULATION." 

PRINT*."  " 

100  PRINT*,"  HOW  MANY  TRIBUTARY  INPUTS  ", 

REA0(1,»)  NIT 

IF ( EOF ( 1 ) )  100,101 

101  IF ( NIT . LT .  1  )  GO  TO  120 

PRINT*,"  INPUT  CODES:” 

DO  It'  1=1 ,NIT 

110  PRINT*,"  ",!,")", 

R  E  A  D  (  1  ,  •  )  I  T  (  I  ) 

IF ( EOF ( 1 ) )  110,111 

111  CONTINUE 

120  PRINT*,"  HOW  MANY  POINT  SOURCE  DISCHARGES  ", 

READl 1 , •)  NIE 

IF (EOFl 1 ) )  120, 121 

121  if(nie.lt.i)  co  ro  i«i 

PRINT*. "  INPUT  CODES: " 
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t  30 

1  40 
1  4  1 

1  42 


1  45 


150 
151 
1  60 


1  6  1 


170 


1  80 

1  8  T 

189 

1  90 
300 


210 


2 


1 


00  140  t  = 1  . N I  £ 

PRINT*,"  ".I,")", 

READ ( 1 , •  )  IS(  I  ) 

IF ( EOF (1))  130,140 
CONTINUE 


BIFURCATION 


INPUT  THE  LENGTH  OF  EACH  PEACH  AND  THE" 
DRAINAGE  AREA  UPSTREAM  FROM  THE  " 
TOP  OF  EACH  REACH  CANS:  MILES  AND" 
SQUARE  MILES) . " 

" , IR , *  )" , 


INDICATE  HOW  YOU  WISH  TO  MODEL  STREAM  " 
HYDRAULICS . « 

0)  MEAN  DEPTH  AND  VELOCITY  SPECIFIED" 

FOR  EACH  REACH  AND  TIME  PERIOD." 

1)  HYDRAULIC  SATING  PARAMETERS  USED  FOR" 
EACH  REACH.” 

(ANS:  USE  EITHER  0  OR  1 )”, 


PRINT*. "  HOW  MANY  MAIN  BRANCH  BIFURCATIONS 

READ ( 1 , * )  NB 
I F ( EOF ( 1 ) )  141,142 

IF(NB.LT. 1)  GO  TO  151 

PRINT*,"  BIFURCATION  CODES:" 

DO  150  Ij  1  , N8 
PRINT*, " 

READ ( 1,915)  BCODE(I) 

IF ( EOF  < 1 ) )  145,150 
CONTINUE 
PRINT*,"  " 

PRINT*, " 

PRINT* , " 

PRINT* , " 

PRINT*, " 

DO  170  I  R  =  1  ,  N  R 
PRINT* , " 

R  EAD ( 1  , •)  LR ( IR ) , DA ( IR ) 

IF ( EOF ( 1 )  )  161,170 
CONTINUE 
PRINT*,"  ” 

PRINT*,"  IV) 

PRINT* , » 

PRINT*, » 

PRINT*, » 

PRINT*, " 

PRINT*, " 

PRINT*, " 

REA0O.314)  HCOOE 
IF ( EOF ( 1  )  )  180,181 
IFCHCOOE. LT. T)  GO  TO  300 

PRINT*,"  INPUT  THE  HYDRAULIC  RATING  PARAMETERS  FOR  EACH" 

PRINT*,"  REACH;  RESPOND  WITH  FOUR  PARAMETERS  IN  " 

PRINT*,"  THE  FOLLOWING  ORDER:  A V , BV , AD , BD . * 

DO  190  1=1, NR 

PRINT*,"  ",I,">", 

REA  D (  1  , • )  HAV(I) ,HBV(I ) ,HAD< I ) ,HBD(I) 

IF ( EOF ( 1 )  )  189,190 
CONTINUE 
PRINT*,"  ■ 

PRINT*. *  * 

PRINT*,"  " 

CALL  LOADTIP(ITP) 

DO  450  ITP* 1 , NTP 
PRINT*,"  " 

PRINT*, "======================*===============================" 

PRINT*,"  TIME  VARIANT  PARAMETERS  FOR  TIME  PERIOD  NO.  ",ITP 
PRINT*, "=====.=============================«=•=================" 

PRINT*,"  " 

PRINT*,"  " 

IF(HCODE.GT.O)  GO  TO  221 

PRINT*,"  I)  HYDRAULIC  PARAMTERS." 

PRINT*,"  " 

PRINT*,"  INPUT  MEAN  VELOCITIES  FOR  EACH  REACH." 

PRINT*,"  ", 

READ ( 1 , •  )  < VEL( II ) , 11= 1 ,NR) 

IF ( EOF ( 1  )  )  210,211 
PRINT*, "  " 

PRINT*,"  INPUT  MEAN  DEPTHS  FOR  EACH  REACH." 
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2  20  PRINT*,"  ", 

READC 1 , • )  ( DEPTH! II ) , II* 1 , NR  ) 

IF ( EOF ( 1 ) )  220,22' 

221  PRINT*,"  " 

PRINT*,  "  II)  BOUNDARY  CONDITIONS  AT  TRIBUTARIES." 

PRINT*,"  " 

PRINT*,"  _  .-T  AMBIENT  WATER  QUALITI  CONDITIONS  FOR" 

PRINT*,"  TRIBUTARIES  1  THROUGH  ",UIT 

DO  22 5  IC*  1  , NTWOC 

224  PRINT*,"  " , CNAHE( I C ) , " : " , 

READ( 1 , • )  ( TWQ ( IIT, IC ) , IIT* 1 , NIT) 

IFvSOFd))  224,225 

225  CONTINUE 
PRINT*,"  " 

2  2  o  PRINT*,"  MEAN  DISCHARGE:", 

REACH,*)  ;  T2(  IIT  )  ,  IITs  1  ,  NIT) 

If ( EOF ( 1 ) )  22 0 , 227 
2  2  7  ?  R I N  T • , "  " 

PRINT*,"  III)  30  UN  DA  R I  CONDITIONS  AT  EFFLUENT  DISCHARGES. 

PRINT*,"  " 

PRINT*,"  . .. . U  T  AMBIENT  WATER  3UALIT  i  CONDITIONS  FOR" 

PRINT*,"  EFFLUENTS  1  THROUGH  ",NIE 

DO  235  IC* 1 , NTWQC 

234  PRINT*,"  " , C  N  A  M  E ( I C ) , " : " , 

R£AD(1,»)  (EWQ( 115, IC) , IIE*  t ,MIE) 

IF (EOF ( 1 ) )  234,235 

235  CONTINUE 
PRINT*,"  " 

236  PRINT*,"  MEAN  DISCHARGE:", 

READ (1 . * )  (EQ( IIE) , IIE* 1 ,NI£) 

IF ( EOF ( 1 ) )  235,237 

237  PRINT*,"  " 

PRINT*,"  IV)  REACTION  RATE  COEFFICIENTS." 

PRINT*,"  " 

PRINT*,"  INPUT  THE  INDICATED  RATE  COEFFICIENT  FOR" 

PRINT*,"  REACHES  1  THROUGH  " , NR 

PRINT*,"  " 

DO  250  IK*  1  ,  NIC 

249  PRINT*,”  ".XNAME(IK), 

READ!!,*)  (K20( IR, IK) , IH* 1 ,Nfl) 

IF  t  EOF (  1  )  )  249,250 

250  CONTINUE 
PRINT*,"  " 

PRINT*,"  V)  DISTRIHUTED  SOURCE/SINK  PARAMETERS." 

PRINT*,"  " 

PRINT*,"  INPUT  THE  INDICATED  SOURCE/SINK  TERM  FOR" 

PRINT*,"  REACHES  1  THROUGH  ",NR 

PRINT*, ”  " 

DO  260  13*1, NS 

259  PRINT*,"  ",5NAME(IS), 

R  EAD (  1 ,  * )  (S20( IH, IS) , IH* 1 ,NB) 

IF ( EOF (  1  ) )  259,260 

260  CONTINUE 
PRINT*,"  " 

PRINT*."  VI)  BIOLOGICAL  PARAMETERS." 

PRINT*,"  " 

PRINT*,"  INPUT  ESTIMATED  ALGAL  BIOMASS  FOR  REACHES" 

PRINT*,"  1  THROUGH  » , NR 

271  PRINT*,"  ", 

R  EA  D (  1 . • )  (A( IR ) , IR* 1 ,Nfl ) 

IF (EOF (  1  )  )  27  1,450 
450  CALL  LOADVPS(ITP) 

CALL  CLOSMS ( 9 ) 


500 


P 


PRINT*," 

" 

PRINT* , " 

It 

PRINT*," 

ft 

PRINT* , " 

It 

PRINT* , " 
PRINT* , « 

n 

DO  XOU  WISH  TO 

REVIEW  AND/OR  REVISE  THE  EXISTING 

PRINT* , * 

DATA  FILE 

n 

PRINT*, * 

(ANS;  i ES 

OR  NO)", 

READ (1,901)  IANS2 
IF ( EOF (  1  )  )  5  10,511 

511  IF  (  IANS2.  WE.  "T"  .  AND  .  IANS2  .  NE.  "N"  )  GO  TO  510 
IF ( IAN32. EQ . "N”  )  00  TO  098 
IF l IAN31 . EQ . "N"  )  00  TO  519 
CALL  0PENMS(9,DINDEX, 125,0) 

519  CALL  UNLOAD(ITP) 

520  PRINT*, "  " 

PRINT*,"  " 

600  CALL  DSUM 

90  1  FORMAT (  1  A  1  ) 

902  FORMAT (111) 

903  F0RMAT(3A10) 

911  FORMAT(SAIO) 

913  FORMAT! 1 A5 ) 

919  FORMAT (111) 

915  FORMAT ( 1F5. 2) 

923  FORMAT ( 8F 10 . 3 ) 

898  CONTINUE 

999  STOP 


INDEX  FOR  ALL  VARIABLES  AND  PARAMETERS  USED  IN  'SIMWQ ' 


A  ( I  R  ) 

ALPHA  C IC  ) 
ATT 

3C0DE  C I B ) 
SETA ( IC , I  ) 
3Q(IB) 

3WQUB,  IC) 
*CNAME( IC  ) 
CNTRLS ( 8 ) 


DA(IR) 
OELTAdC,  I) 
DEPTH(IR) 

D INDEX ( I ) 
EPS  I L ( IC, I ) 
EQ(IE) 

EWQ ( IE , IC ) 
GAMMA ( IC, I ) 
H A D ( IR  ) 

H A V ( IR  ) 

HBD ( IR  ) 

H3V ( IR ) 
HCODE 
IA  TT 
13 
IC 

lE(IIE) 

INUM 


ALGAL  CONCENTRATIONS  IN  REACH  'IR'  (MG/L) 

PARAMETER  FOR  ANALTTICAL  W Q  SOLUTION  (MG/L) 

KEK  FOR  MASS  STORAGE  OF  WQ  VECTOR 
BIFURACTION  CODE  DEFINING  WATERSHED  BRANCHES 
PARAMETER  FOR  ANALTTICAL  WQ  SOLUTION  (MG/L) 

DISCHARGE  AT  MINOR  TRIBUTART  AT  IB'TH 
BIFURCATION  (FT**3/SEC) 

WQ  VECTOR  AT  MINOR  TRIBUTART  AT  IB'TH  BIFURCATION 
ALPHABETIC  NAMES  OF  WATER  QUALITT  ATTIBUTES 
ARRAT  HOLDING  SIMULATION  CONTROL  PARAMETERS  FOR 
TRANSFER  TO  AND  FROM  MASS  STORAGE 
(NH,NTP,NTWQC,£TC. ) 

DRAINAGE  AREA  UPSTREAM  FROM  IR'TH  REACH  (MILES**2) 
PARAMETER  FOR  ANALTTICAL  WQ  SOLUTION 
MEAN  DEPTH  OF  IR'TH  REACH  (FEET) 

MASTER  INDEX  FOR  MASS  STORAGE  OF  SIMULATION  DATA 
PARAMETER  FOR  ANALTTICAL  WQ  SOLUTION 

DISCHARGE  AT  IE'TH  POINT  SOURCE  EFFLUENT  INPUT  (FT*3/SEC) 

WATER  QUALITT  VECTOR  AT  IE'TH  EFFLUENT  INPUT 

PARAMETER  FOR  ANALTTICAL  WQ  SOLTIONS 

HTDRAJLIC  RATING  PARAMETER  FOR  MEAN  DEPTHS 

HTDRAULIC  RATING  PARAMETER  FOR  MEAN  VELOCITIES 

HTDHAUL IC  RATING  PARAMETER  FOR  MEAN  DEPTHS 

HTDRAULIC  RATING  PARAMETER  FOR  MEAN  VELOCITIES 

SWITCH  FOR  SPECIFTIWG  MODELING  CHOICE  FOR  DEPTH  VEL 

ATTRIBUTE  INDEX  (1  TO  20) 

BIFURACTION  INDEX  (1  TO  5) 

WQ  CONSTITUENT  INDEX  (>  TO  20);  SAME  AS  I A  TT 
EFFLUENT  INPUT  INDEX  ('  TO  15) 

INDEX  FOR  SPECIFICATION  OF  RECORD  KET  FOR  MASS  STORAGE 
OF  WQ  VECTORS 


SAME  AS  IPP 


IP 


IPP 

la 

i  a  e  c 
irscip 

ITCIT) 

X  TO  P 
ITP 

I  ago 

•<T 

.<20 

La  ua ) 

M C 0 D E  X  IC  ) 

NAATTS 

M3 

MCWQC 

Nit 
N  XT 
NX 

NNCWQC 

HP? 

N.a 

NS 

n  r? 

NTWQC 

NWQC 

PH 

PP 

g(ipp) 

QB 

QE 

0 1 N  0  E  X  (  I  ) 

OT 

go 

QIOK(IX) 

gios(is) 
ST( IR , IS) 
320 ( IH ,  IS ) 
TEMP 

T ITLE ( 3 ) 
TLR 

T0( IT) 

TTWR 

TWg( IT, IC) 
7  E  L  (  I  R  ) 

wq( ipp) 
wqb( is , ic ) 
WQ£(IE,IC) 
WQT  ( IT . IC  ) 

wqo ( ic ) 

X ( IPP, 3  ) 

X  I B 
XINC 
xi  a 

X  L  A  3  i 

XWR 

SETA 


INDEX  OF  PROFILE  POINT  IN  SIMULATION  RESULTS  (1  TO  200) 
INDEX  OF  REACH  (1  TO  20) 

RECORD  SET  INDEX  FOR  MASS  STORAGE  OUTPUT  OF  WO  VECTOR 
INDEX  OF  RECIEVING  REACH  (DOWNSTREAM)  OF  I3'TH  BIFUCATION 
INDEX  OF  INPUT  LOCATION  OF  TRIBUTARIES 
INDEX  SPECIFYING  TO?  OF  SUPERIOR  CHANNEL  BRANCHES 
INDEX  OF  TIME  PERIOD  (1  TO  12) 

INDEX  OF  WATER  DUALITY  CONSTITUENTS  (1  TO  20) 

TEMPERATURE  ADJUSTED  KINETIC  RATE  COEFFICIENT  (I /DA!.' 
KINETIC  RATE  COEFFICIENT  AT  20  DEGREES  C  (l/DA I) 

LENGTH  OF  IR'TH  REACH  (MILES) 

MODELING  CODE  FOR  IC'TH  ATTRIBUTE 
NUMBER  OF  WQ  ATTRIBUTES 

NUMBER  OF  BIFURCATIONS  OF  STREAM  CHANNEL  (0<N8<5) 

NUMBER  OF  CONSERVATION  WO  CONSIT’JENTS  (  0  <N  C  WQC  <  1  2- N  NC  WQ  C  ) 
NUM3ER  OF  POINT  SOURCE  EFFLUENT  INPUTS  (0<NIE<15) 

NUMBER  OF  TRIBUTARY  INPUTS  (1<«IT<I5) 

NUMBER  OF  RATE  COEFFICIENTS  NECESSARY  FOR  SIMULATION 
NUMBER  OF  NON-CONSERVATION  WQ  CONSTITUENTS 
NUMBER  OF  PROFILE  POINTS  IN  SIMULATION  OF  CURRENT 
NUMBER  OF  REACHES  IN  CURRENT  SIMULATION  (1<NR<20) 

NUMBER  OF  SOURCE/SINK  TERMS  IN  SIMULATION  (0<NS)20, 

NUMBER  OF  TIME  PERIODS  IN  SIMULATION  SUN  (1<NT?<12; 

NUMBER  OF  TOTAL  Wq  CONSTITUENTS  IN  SIMULATION  (NTWQCOO) 
NUMBER  OF  CORE  WATER  OUALTIY  CONSTITUENTS  IN  SIMULATION 
PERCENT  OF  NITROGEN  IN  SUSPENDED  SOLIDS  (1/100: 

PERCENT  OF  PHOSPHORUS  IN  SUSPSENDED  SOLIDS  ',1/100) 

STREAM  DISCHARGE  AT  IPP'TH  POINT  OF  OUTPUT  PROFILE  ( f  T • • 3 / 
DISCHARGE  INPUT  FROM  BIFURCATION  TO  CURRENT  REACH  CFG) 
DISCHARGE  INPUT  FROM  TRIBUTARIES  TO  CURRENT  REACH  iCF3) 
MASTER  INDEX  FOR  MASS  STORAGE  OF  SIMULATED  WO  PROFILES 
DISCHARGE  INPUT  FROM  TRIBUTARIES  TO  CURRENT  REACH  (CFS) 
DISCHARGE  AT  TOP  OF  CURRENT  REACH  AFTER  DILUTIONS  (CFS) 
TEMPERATURE  CORRECTION  FACTOR  FOR  IX'TH  RATE  COEFFICIENT 
TEMPERATURE  CORRECTION  FACTOR  FOR  IS'TH  SOURCE/SINK  TERM 
TEMPERATURE  CORRECTED  30URCE/SINX  TERM  FOR  IR'TH  REACH 
SOURCE/SINK  TERM  AT  20  DEGREES  C  (MG/L/DAY) 

WATER  TEMPERATURE  (DEGREES  C) 

SIMULATION  NAME 

TOTAL  LENGTH  OF  STREAM  TO  BE  SIMULATED  (MILES) 

DISCHARGE  IN  IT'TH  TRIBUTARY  (CFS) 

TOTAL  TIME  OF  TRAVEL  WITHIN  CURRENT  REACH  (DAYS) 

WATER  QUALITY  VECTOR  FOR  IT'TH  TRIBUTARY  INPUT 
MEAN  VELOCITY  IN  IR'TH  REACH  (MILES/DA i  ) 

WATER  QUALITY  VECTOR  AT  IPP'TH  POINT  IN  SIMULATION 
WATER  QUALITY  VECTOR  FOR  IB'TH  BIFURCATION  INPUT 
WATER  QUi  ITY  VECTOR  FOR  IS'TH  EFFLUENT  INPUT 
WATER  QUALITY  VECTOR  FOR  IT'TH  TRIBUTARY  INPJT 
WATER  QUALITY  VECTOR  AT  TOP  OF  CURRENT  REACH  AFTER 
DISTANCE  ARRAY  SPEC  IF  INC  DISTANCE  DOWNSTREAM.  REACH 
INDEX  OF  LOCATION  ON  WATERSHED  BRANCHES 
INCREMENT  BETWEEN  PROFILE  POINTS  IN  SIMULATION  OUTPUT 
INDEX  OF  LOCATION  IN  TERMS  Or  STREAM  REACH 
DISTANCE  AT  TOP  OF  CURRENT  REACH 
DISTANCE  FROM  TOP  OF  CURRENT  REACH 
PARAMETER  FOR  ANALYTICAL  WQ  SOLUTION 


SUBROUTINE  DSUM 


COMMON  /  N  AMES/TI  rLEi  3  )  ,  CNAMEt  20)  ,  SHAME  (  2  0  )  ,  KHAME  (  2  0  ) 

COMMON /CN TRL/NR , NT? ,  NWOC , Nfl , NIT . NIS . IT( 1 5 ) . IE( 1 5 ) . HCODE , 

♦  MCODE(d) , NCWQC , NHCWQC , NTWQC , NS , NK 

COMMON/ I WQC/TWQ( 15,20), £W0( 1 5 , 20  )  , TQ (  1 5  )  , EQ ( 1 5  ) 
C0MM0N/KDATA/K2G(20,20) 

C0MMCN/SDATA/S20( 20, 20) 

COMMON /MDATA/LR( 20 ) , DA ( 20 ) , DE PTH ( 20 ) , VEL ( 20 ) , BCODE (  1 5 )  , 

♦  HAV(20 ) ,H3V<  20 ), HAD ( 2C) ,K3D( 20) 

COMMON /BDATA/A( 20) 

REAL  K20.LR 

INTEGER  TITLE.CN AM £, HCODE, KN A ME, 3NAME 
INTEGER  T IN  PUT, EIN PUT .START 
PRINT*,"  " 

PlIilT1,'::::::::: . . . . . . 

PRINT*,"  REVIEW  Of  CONTENTS  OF  CURRENT  DATA  FILE" 

FlLiT','::::::::::::::::: . . . . 

PRINT*,"  " 

PRINT*,"  THE  CURRENT  CONTENTS  OF  IOUR  DATA  FILE  HAVE  THE" 

PRINT*,"  TITLE:" 

PRINT*,"  " 

PRINT*,"  " , ( TITLE ( I ) , Is  1 , 5 ) 

PRINT*,"  ", (TITLE(I),I=6,d) 

PRINT*,"  " 

PRINT*,”  THIS  DATA  SET  SPECIFIES  SIMULATION  OF  THE" 

PRINT*,"  FOLLOWING  WATER  QUALITY  ATTRIBUTES:" 

PRINT*,"  ", (CNAME( I) , 1=1 ,5) 

PRINT* , "  ( CNAME ( I ) , 1=6 , 1 0) 

IF ( NTWQC . LT . 1 1 )  GO  TO  10 
PRINT*,"  ", (CNAME(I) , 1=1 1 . 15) 

IF ( NTWQC . LT . 16)  GO  TO  10 
PRINT*,"  ", (CNAME(I) , Is  16,20) 

10  PRINT*,"  » 

PRINT*,"  SIMULATIONS  WILL  3E  RUN  FOR  " , N T P , "  TIME  PERIODS" 

PRINT*,"  Fn  "  '  TOTAL  OF  ",NR,"  STREAM  REACHES." 

CALL  WAIT 

PRINT*,"  THE  > 1  *  i  3  ? .  WATERSHED  STRUCTURE  IS  AS  FOLLOWS:" 

PRINT*,"  " 

PRINT*,"  R L r  . LENGTH  DRAINAGE  INPUTS" 

PRINT*,"  NO.  (MI.)  (SO. MI.)  EFF .  TRIB." 

PRINT*,"  " 

DO  30  I R  = 1 , N  R 
TINPUT="  " 

EINPUTs"  " 

DO  23  I I T= 1 , N IT 

IF(IT(IIT).EQ.IR)  TINPUTsIT(IIT) 

23  CONTINUE 

DO  26  IIEs 1 , N I E 

IF(IE(II£).£Q.IR)  EINPUTsIECIIE) 

26  CONTINUE 

30  PRINK  1 , 900  )  IR , LR ( IR ) , DA( IR ) , E IN  PUT , T I N PUT 
900  FORMAT ( 7X, 1I2,3X,2(1F!0.3,2X),7X,1I2,6X,!I2) 

PRINT*,"  " 

PRINT*,"  THE  NUMBER  OF  MAJOR  BIFURCATIONS  OF  THE  MAIN” 

PRINT*,"  CHANNEL  OF  THIS  RECIEVING  STREAM  IS  " , NB 

IF(NB.LE.O)  GO  TO  40 
STARTs  1 
PRINT*,"  " 

DO  39  INBsl.NB 

IEND=BC0DE(INB)*100-tFIX(8CODE(IN8))*100 

PRINT*,"  BRANCH  NO.  ",IN3,"  INCLUDES  REACHES  ".START,"  TO  ", 

*  I  END 
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39 


STARIdEND+l 
PRINT*. "  " 

PRINT*,"  THE  MAIN  CHANNEL  INCLUDES  REACHES  ".START,"  TO  " , NR 

CALL  WAIT 
40  IF(HCODE.GT.O)  GO  TO  50 

PRINT*,”  HYDRAULIC  MODELING  WILL  BE  DONE  USING  MEAN  * 

PRINT*,"  VELOCITIES  AND  DEPTHS  f'OR  EACH  REACH  AND  TIME" 

PRINT*,"  PERIOD." 

GO  TO  60 

50  PRINT*,"  HYDRAULIC  MODELLING  WILL  BE  DONE  USING  HYDRAULIC  " 

PRINT*,"  RATING  PARAMETERS  RELATING  MEAN  VELOCITY  AND  " 

PRINT*,”  DEPTHS  TO  DISCHARGE  IN  A  REACH." 

o 0  CONTINUE 

PRINT*,"  « 

PRINT*,"  NONE  OF  THE  ABOVE  PARAMTERS  CAN  BE  ALTERED  WITHOUT" 

PRINT*,"  CREATING  A  TOTALLY  NEW  DATA  SET  (I.E.,  BY" 

PRINT*,"  STARTING  OVER  WITH  ' F Q U E S ' ) . " 

PRINT*. "  " 

70  PRINT*,"  DO  YOU  WANT  TO  CONTINUE  (ANS:  YES  OR  NO)", 

„^AD(1,90n  IANS 

901  FORMAT ( 1  A  1  ) 

IF ( EOF ( 1 ) )  70,71 

71  IFdANS.NE.  "Y". AND. IANS. NE.  "N")  GO  TO  70 
IF  (  IANS.  Eq  .  "Y"  )  GO  TO  100 

RETURN 

100  PRINT*,"  " 

PRINT*,"  WHICH  TIME  PERIOD  DO  YOU  WANT  TO  REVIEW", 

R  E A  D (  1  , • )  ITP 

IF  C  EOF (  1  ) )  100,  101 

It (I7P.LT. 0. OR. ITP.GT.NTP)  GO  TO  100 

101  CALL  UNLTVPS(ITP) 

PRINT*,"  " 

110  PRINT*,"  INDICATE  WHICH  OF  THE  FOLLOWING  PARAMETERS  YOU" 

PRINT*,"  WANT  TO  REVIEW  (RESPOND  WITH  THE  PROPER  NUMBER)" 

113  PRINT*,"  1)  HYDRAULIC  PARAMETERS" 

PRINT*,"  2)  INITIAL  CONDITIONS  IN  TRIBUTARIES" 

PRINT*,"  3)  INITIAL  CONDITIONS  IN  EFFLUENTS" 

PRINT*,"  4)  KINETIC  PARAMETERS" 

PRINT*,"  5)  DISTRIBUTED  SOURCE/SINK  PARAMETERS" 

PRINT*,"  6)  BIOLOGICAL  PARAMETERS" 

1 1 1  PRINT*, "  ", 

R  E  A  D (  1  , • )  IPARAM 

I F ( EOF ( 1 ) )  111,112 

112  IFdPARAM.LT.  1  .  OR  .  I  PAR  All .  GT  .  6  )  GO  TO  111 
GO  TO  ( 150 . 200, 250, 300, 350, 400) , IPARAM 

150  PRINT*,"  " 

IF ( HCODE . NE . 1 )  GO  TO  15! 

PRINT*,"  THIS  DATA  SET  USES  HYDRAULIC  RATING  PARAMETERS" 

PRINT*,"  THEREFORE  NO  MEAN  VALUES  FOR  VELOCITY" 

PRINT*,"  OR  DEPTH  ARE  NEEDED." 

PRINT*,"  " 

GO  TO  500 

151  PRINT*,"  THE  CURRENT  VALUES  FOR  MEAN  VELOCITIES  AND  DEPTHS" 

PRINT*,"  FOR  EACH  REACH  DURING  TIME  PERIOD  " , ITP , "  ARE" 

PRINT*,"  LISTED  BELOW.  TO  CHANGE  THE  CURRENT  VALUES" 

PRINT*,"  RESPOND  TO  TRAILING  '?'  BY  THE  NEW  VALUES.” 

PRINT*,"  (NO.  REACH, VELOCITY, DEPTH)" 

DO  170  IR=1,NR 

PRINT*,"  " ,  I  R , "  )  " , VEL( IR ) , DEPTH( IR ) , 

R  E  A  D (  1  , • )  VEL( IR ) , DEPTH( IR ) 

IF ( EOF (  1  )  )  169,'70 

169  CONTINUE 

170  CONTINUE 
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CO  TO  500 
200  PRINT*,"  “ 

PRINT*,"  INDICATE  WHICH  NAT  E  R  QUALITY  ATTRIBUTE  YOU  ARE" 

PRINT*,"  INTERESTED  IN:“ 

209  DO  210  I  C  :  1  , NTWQC 

210  PRINT*,"  ",IC,")  " , CNAME ( IC  ) 

PRINT*,"  ",NTWQC*1.")  DISCHARGE" 

PRINT*,"  ", 

21  1  READ (  1  , •)  I WQC 

I? ( EOF ( 1 ) )  211,212 
212  IF(IWQC.E3.NTWQC*1 )  GO  TO  230 

PRINT*,"  THE  VALUES  FOR  " , CNAMEC IWQC ) . "  FOR  TIME  PERIOD  " 

PRINT*,"  " , I T  P , "  ARE  LISTED  BELOW  BY  TRIBUTARY. " 

PRINT*,"  TO  CHANGE  A  VALUE  RESPOND  TO  THE  TRAILING 

PRINT*,"  WITH  THE  NEW  VALUE." 

DO  220  I  IT: 1  , NIT 

PRINT*,"  "  ,  I  IT  ,  "  )  " , TWQ ( I IT , IWQC ) , 

RSADC 1 , • )  TWQ( IIT, IWQC ) 

IF ( EOF (  1  )  )  220 ,219 

219  CONTINUE 

220  CONTINUE 
PRINT*,"  " 

GO  TO  241 

230  PRINT*,"  THE  VALUES  FOR  TRIBUTARY  DISCHARGE  FOR  TIME  PERIOD" 

PRINT*,"  "  ,  IT  P , "  ARE  LISTED  BY  TRIBUTARY  BELOW.  TO" 

PRINT*,"  CHANGE  A  VALUE  RESPOND  TO  THE  TRAILING  '?'  BY" 

PRINT*. "  THE  NEW  VALUE." 

DO  240  IIT: 1 , NIT 

PRINT*,"  " , II T , ”  )  " , TQ ( I IT ) , 

READd,*)  TO  (  IIT) 

IF  (  EOF  (  1  )  )  I 9  ,  -  4U 

239  CONTINUE 

240  CONTINUE 

241  PRINT*,"  " 

PRINT*,"  REVIEW  ANOTHER  TRIBUTARY  ATTRIBUTE  (ANS:  YES  OR  NO)" 

-ADC  1 , 901 )  I A  N  S  4 
IF (EOF  (  1  )  )  24  1  ,  242 

242  IFCIANS4. HE. "H" . AND. IANS4.NE. "Y" )  GO  TO  241 
IFCIANS4.EQ." N")  GO  TO  500 

243  PRINT*, "  WHICH  ATTRIBUTE", 

READ ( 1 , • )  IWQC 

IF ( EOF ( 1 ) )  243,244 

244  IFCIWQC.LT. 1 .0R.IWQC.GT.NTWQC+1)  GO  TO  209 
GO  TO  212 

250  PRINT*,"  " 

PRINT*,"  INDICATE  WHICH  WATER  QUALITY  ATTRIBUTE  YOU  ARE" 

PRINT*,"  INTERESTED  IN:” 

259  DO  260  IC : 1 , NTWQC 

260  PRINT*,"  ",IC,">  ", CNAME (IC) 

PRINT*,"  " , NTWQC* 1 , " )  DISCHARGE" 

PRINT*,"  ", 

261  READ  < 1 , • )  IWQC 

IF ( EOF ( 1 ) )  261,262 

262  IF ( IWQC . EQ . NTWQC* 1 )  GO  TO  280 

PRINT*,"  THE  VALUES  FOR  ", CNAME C IWQC ), "  FOR  TIME  PERIOD  " 

PRINT*,"  " , IT  P , "  ARE  LISTED  BELOW  BY  EFFLUENT." 

PRINT*,"  TO  CHANGE  A  VALUE  RESPOND  TO  THE  TRAILING 

PRINT*,"  WITH  THE  NEW  VALUE." 

DO  270  1  1 E : 1 , N I  £ 

PRINT*,"  " , 1 1 E , "  )  ", EWQC IIE, IWQC) , 

READ ( 1 , *  )  EWQ (IIE, IWQC ) 

Ir ( EOF (  1  )  )  270,269 

269  CONTINUE 

270  CONTINUE 
PRINT*,"  " 

GO  TO  29' 
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280  PRINT*,"  THE  VALUES  FOR  EFFLUENT  DISCHARGE  FOR  TIME  PERIOD" 

PRINT*,"  " , I T  P , "  ARE  LISTED  3*  EFFLUENT  BELOW.  TO" 

PRINT*."  CHANGE  A  VALUE  RESPOND  TO  THE  TRAILING  B  X  " 

PRINT*,"  THE  MEW  VALUE." 

DO  290  HE*  1  ,  NIE 

PRINT*,"  ".HE,")  "  ,  £Q  ( 1 1 E  )  , 

READd,*)  EQ(IIE) 

IF ( EOF (  1  )  )  289,290 

289  CONTINUE 

290  CONTINUE 

291  PRINT*,"  " 

PRINT*,"  REVIEW  ANOTHER  EFFLUENT  ATTRIBUTE  (ANS:  XES  OR  NO)", 

READ (1,901)  IANS4 
IF ( EOF (  1  )  )  291 ,292 

292  IF  ( IANS4  .  ME  .  "N"  .  AND  .  IAN54.  NE  .  "X"  )  GO  TO  291 
IF  (  IANS4  .  EQ  .  "N"  )  GO  TO  500 

293  PRINT*,"  WHICH  ATTRIBUTE", 

READd,*)  I WQC 

I F ( EOF ( 1 ) )  293,294 

294  IF ( iwqc . LT.  I  . OR . IWQC .GT .NTWQC* 1  )  GO  TO  259 
GO  TO  262 

300  PRINT*,"  " 

PRINT*,"  INDICATE  THE  INDEX  NUMBER  OF  THE  KINETIC  PARAMETER" 

PRINT*,"  XOU  ARE  INTERESTED  IN  REVIEWING" 

309  DO  310  IKsl.NK 

310  PRINT*,"  * , IK , » )  * , KNAME ( IK ) 

311  PRINT*,"  ", 

READ  C  1 , • )  IK 

IF ( EOF (  1  ))  311,312 

312  IF ( IK.LT. 1 .OR.IK.GT.NK)  GO  TO  311 
PRINT*,"  " 

313  PRINT*,"  THE  VALUES  CURRENTLX  SPECIFIED  FOR  ", KNAME ( IK ), "WILL " 

PRINT*,"  BE  LISTED  BELOW  BX  REACH.  TO  CHANGE  A  VALUE* 

PRINT*,"  RESPOND  TO  THE  TRAILING  *?*  WITH  THE  NEW" 

PRINT*,"  VALUE." 

DO  320  IRsl.NR 

PRINT*. "  " , IR , " >  " , K20 ( IR , IK ) , 

READd,*)  K20  ( I R  ,  IK  ) 

IF ( EOF ( 1  )  )  319,320 

319  CONTINUE 

320  CONTINUE 

321  PRINT*. "  REVIEW  ANOTHER  RATE  COEFFICIENT  (ANS:  XES  OR  NO)", 

READd, 901)  IANS5 

IF (EOF (  1  ) )  321  ,  322 

322  IF( IANS5 . NE. "X" . A. '.IANS5 . NE. "N" )  GO  TO  321 
IF ( IANS5 . EQ . "N" )  GO  TO  500 

330  PRINT*,"  WHICH  ONE", 

R  EA  D (  1  , •)  IK 

IF ( EOF (  1  ) )  330,331 

331  IF ( IK ■ LT . 1. OR.IK.GT.NK)  GO  TO  309 
GO  TO  313 

350  PRINT*,"  " 

PRINT*,"  INDICATE  THE  INDEX  NUMBER  OF  THE  DISTRIBUTED" 

PRINT*,"  SOURCE/SINK  TERM  XOU  WISH  TO  REVIEW" 

359  DO  360  IS*  1  , NS 

360  PRINT*. "  ",IS,">  " , SNAME ( IS ) 

361  PRINT*,"  ", 

READd,*)  IS 

IF ( EOF (  1  )  )  361,362 

362  IFUS.LT.O.OP.  IS.GT.NS)  GO  TO  361 

363  PRINT*."  " 

PRINT*,"  THE  VALUES  CURRENTLX  SPECIFIED  FOR  " , SN AME ( I S ) , " W I LL " 

PRINT*. "  BE  LISTED  BELOW  BX  REACH.  TO  CHANGE  A  VALUE" 

PRINT*,"  RESPOND  TO  THE  TRAILING  *7*  WITH  THE  NEW" 

PRINT*,"  VALUE." 
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DO  370  IRsl.NR 

PRINT*,"  " , I R  ,  "  )  ",S20(IR,IS), 

READO,*)  S20CIR.IS) 

IF ( EOF (  1  ) )  369,370 

369  CONTINUE 

370  CONTINUE 

371  PRINT*,"  REVIEW  ANOTHER  SOURCE/SINK  TERM  (ANS:  IES  OR  NO)", 

READO, 901)  IANS6 

IF ( EOF ( 1 ) )  371,372 

372  IF(IANS6. NE . "N" . AND . IANS6 .NE . "I" )  GO  TO  371 
IF(IANS6.£Q."N")  GO  TO  500 

330  PRINT*,"  WHICH  ONE", 

READO,*)  IS 

IF ( EOF  Cl))  380,381 

381  IF < IS. LT.1 .OR. IS.GT.NS)  GO  TO  359 
GO  TO  363 
400  PRINT*,"  " 

PRINT*,"  THE  CURRENT  VALUES  FOR  MEAN  ALGAL  CONCENTRATIONS" 
PRINT*,"  ARE  LISTED  BELOW  BI  REACH  AND  FOR  TIME  PERIOD" 

PRINT*,"  " , ITP , " .  TO  CHANGE  A  VALUE  RESPOND  TO  THE  " 

PRINT*,"  TRAILING  B*  THE  NEW  VALUE." 

DO  410  IR.1.NR 

PRINT*,"  " , I R , " )  " , A ( IR ) , 

READO,*)  A  ( IR  ) 

IF ( EOF (  1  ) )  409,410 

409  CONTINUE 

410  CONTINUE 
PRINT*,"  " 

500  PRINT*,"  " 

PRINT*,"  CONTINUE  FOR  THIS  TIME  PERIOD  (ANS:  XES  OR  NO)", 

501  READO, 901)  IANS  1 
IF ( EOF (  1  )  )  501,502 

502  IF(IANS1 ,NE. "T" . AND. IANS1 .NE. "N")  GO  TO  501 
IF(IANS1.EQ."N")  GO  TO  510 

PRINT*,"  " 

508  PRINT*,"  INDICATE  PARAMETER  TT  PE  ", 

READ*, IPARAM 

IF (EOF (  1  ) )  508,509 

509  IFdPARAM.LT.O.OR.  IPARAM.  GT.  6)  GO  TO  113 
GO  TO  (  150,200,250, 300, 350, 40Q) , IPARAM 

510  PRINT*,"  " 

PRINT*,"  REVIEW  ANOTHER  TIME  PERIOD  (ANS:  YES  OR  NO)", 

511  READ(1,901)  IANS2 
IF ( EOF (  1  ) )  511,512 

512  IF(IANS2.NE."7".AND.IANS2.NE."N")  GO  TO  511 
CALL  LOADV PS ( ITP ) 

IF(IANS2.EQ."X")  GO  TO  100 

RETURN 

END 


72 


570 


580 


COMMON /N AMES /T ITLE(3 ) .CNAME (20) ,SNAME( 20) , KN AME ( 20 ) 
COMMON/CNTRL/NR,NTP,NUQC,NB,NIT,NIE, IT( 15) , IE ( 1 5 ) , HCODE , 

■  MCODE ( 8 ) , NCWQC , NNCWQC , NT'AQC , NS , NK 
COMMON /IHQC/TWQ(  15,20) ,£VQ(  1 5 , 20 ) , TQ ( 1 5 ) , EQ( 1 5 ) 

COMMON/ KDATA/K20( 20, 20) 

COMMON /SDATA/S 20 ( 20 ,20) 

COMMON /MDATA/LR( 20) , DA (20 ) .DEPTH (20) , V  EL ( 20 ) , BCODE ( 15) , 

•  HAV(20) ,HBV(20) .HAD (20) ,HBD(20) 

C0MMON/BDATA/A(20 ) 

REAL  LR.K20 

INTEGER  TITLE, HCODE, CNAME 
INTEGER  DIND£X( 1 25 ) , CNTRLSt 8 ) 

DATA  ( CNAME ( I ) , 1  =  1  , 8  >/"TEMP. " , "B0D5" , "TSS" , "NH  3 " , "N02" , "NO  3" 

•  "P04" , "D.O. "/ 

•>  LOAD  DATA  INTO  MASS  STORAGE  FILE 
CALL  WRITM3( 9, TITLE, 3, 1 ,-1 ) 

CNTRLS ( 1 )*NR 
CNTR LS ( 2  )  i NT  P 
CNTRLS( 3 ) =NCWQC 
CNTRLS(4)=NNCWQC 
CNTRLS(5)sHC0DE 
CNTRLS ( 6 )  =  NB 
CNTRLS(7)=NIT 
CNTRLS ( 8 ) sNI E 

CALL  HRITMS(9, CNTRLS, 8, 2,-1 ) 

CALL  HRITMS(9, MCODE, 8, 3.-1) 

NAATTSsNCHQC*NNCHQC 
IF(NAATTS.LE.O)  GO  TO  570 
CALL  WRITMS(9, CNAME, NTWQC, 4,-1) 

CALL  WRITM3(9,BCODE,NB,5,-1  ) 

CALL  WRITMS(9,IT,NIT,6,-1) 

CALL  NRITMSt  9 , IE . TIE, 7,  -  ■  ' 

CALL  WRITMS(9, LR, NH , 8, -1 ) 

CALL  WRITMS(9,DA,NR,9,-1) 

IF(HCODE.EQ.O)  GO  TO  580 
CALL  NRITMS(9, HAV, NR, 10,-1 ) 

CALL  WRITMS(9.HBV, NR, 11,-1) 

CALL  WRITMS( 9, HAD , NR , 12, -1 ) 

V.lL  WRI TMS  (9.HBD.NR,  13,-1) 

RETURN 

t  NTR  K  LOADVPS 

/ TWQCsU+NCWQC+NNCWQC 

NK= 16+NNCWQC 

NSsNTWQC 

II=13*(ITP-1 ) *9 

IF (HCODE. EG. 1 )  GO  TO  535 

INUMs I !♦ 1 

CALL  WRITMS(9, VEL.NR, INUM, -1 ) 

I NUM= I I *2 

CALL  NR  I TMS (9, DEPTH, NR, INUM, -1) 

NUM=NIT»NTNQC 

INUM=II*3 

CALL  WRITMS(9,TWQ, 300, INUM, -1) 

INUM=II*4 

CALL  WRITMSt 9,TQ,NIT, INUM, -1 ) 

N  UM*N I E •HTNQC 
I  NUM  » 1 1 *5 
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CALL  WRITMS(  9,  E'NO,  300  ,  INUM,  -  1  > 

XNUM»XI»6 

CALL  WRITMS(9,EQ,NIE,INUM,-1 ) 

NUMsNK*NR 

INUMsII*7 

CALL  tf RI TMS (9, K20 , 400, INUM , - 1 ) 

NUMsNS»NR 

INUMsIIfB 

CALL  WRITMS(9,S20, 900, INUM, -1) 

INUM*II*9 

CALL  WRITMS(9, A,NR,INUM,-1  ) 

RETURN 

END 

SUBROUTINE  UNLOAD(ITP) 

C3ssi33itiis3sim3isi*3ui:ii«3s33s:iss:i3i:stsii33saii:<3ss3s 

COMMON/NAM£S/TITLE(8),CNAME(20),SNAM£(20>,KNAME(20) 

COMMON /CNTRL /N  R ,NTP,NWQC,NB,NIT,NIE,IT( 1 5 ) , I E (  15),HC0DE, 

*  MC0DE(8),NCWQC,NNCWQC,NTWQC,NS,NK 

COMMON / 1 NQC /TWO ( 15,20) ,EWQ( 1 5 , 2 0 ) , TO ( 1 5 ) , EQ ( 1 5  ) 
COMMON/XDATA/K20(20,20) 

COMMON /SDATA/S20( 2 0,20) 

COMMON/MDATA/LR(20) , DA (20) .DEPTH (20) , VEL(20) ,BCODE( 15) , 

♦  HA V ( 20 ) ,HBV(20) , HAD (20) ,HBD(20) 

C0MM0N/BDATA/A(20 ) 

REAL  LR.X20 

INTEGER  TITLE, HCODS , C3AME 
INTEGER  DIND£X{ 125) ,CNTRLS(8) 

C - >  UNLOAD  DATA  FROM  MASS  STORAGE  FILE 

CALL  READMS(9, TITLE, 8, 1 ) 

CALL  READMS(9,CNTRLS,8,2) 

NR*CNTHLS( 1 ) 

NTP*CNTRLS ( 2 ) 

NCWQC*CNTRLS(3) 

NNCWQC*CNTRLS(U) 

HC0DEsCNTRLS(5) 

NBsCNTRLS ( 6 ) 

NIT*CNTRLS(7) 

N I EsCNTRLS ( 8 ) 

NUQC  >  8 

NIC* 1 6+NNCWQC 
NS*8-*NNCHQC*NCWQC 
NTWQC«8»NCWQC*NNCWQC 
CALL  READMS (9,MCODE,8,3) 

NAATTS*NCWQC+NNCHQC 
IF(NAATTS.LE.O)  GO  TO  570 
CALL  READ '■*S(9,C NAME, NT WQC,  9  ) 

570  CALL  REA0MS(9,BC00E,N8, 5) 

CALL  READMS(9,IT,NIT,6) 

CALL  READMS (9,IE,NIE,7) 

CALL  READMS(9,LR,NR,8) 

CALL  READMS (9, DA, NR, 9) 

IP(HCODE. EQ. 0)  GO  TO  580 
CALL  READMS( 9, HAV , NR , 10) 

CALL  REA0MS(9,HBV,NR,  1)  ) 

CALL  REA0M3(9,HA0,NR, 12) 

CALL  READMS(9,HBD,NR, 1 3) 

580  RETURN 
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X 


ENTRT  UNLTVPS 

11=1  3*( ITP-1 ) *9 

IF ( HCODE . EQ . 1 )  GO  TO  585 

INUM=II*1 

CALL  READMS (9,VEL,NR,INUM) 

INUM=II*2 

CALL  REAOMS(9, DEPTH, NS, INUM) 

N  UM  =  N IT  *NTWGC 
INUM=II+3 

CALL  READMS(9,TWQ, 300, INUM) 

INUM  =  II**t 

CALL  READMS(9,TQ,NlTtINUM) 

NUM=NIE*NTVQC 

INUM=II*5 

CALL  READMS(9,SUa, 300.INUM) 

INUM=II*6 

CALL  READMS(9,SQ,NI£,INUM) 

NUM=NK«NH 
I MUM =11*7 

CALL  READMS(9 , K20 , 400 , INUM) 

NUM=NS»NH 
INUM= I 1*8 

CALL  READMS (9, 320, 400, INUM) 

INUM=II*9 

CALL  READMS (9, A, NR, INUM) 

RETURN 

END 

:3333=======5===33=33=3==3==333==3====3====333= 

SUBROUTINE  WAIT 

[33X3333333333 

READ* , DUMMI 
IF  { EOF ( 1  )  )  1 ,2 
GO  TO  2 
DO  3  1=1,5 
PRINT*,"  " 

RETURN 
ENO 
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o  n  o  ooo  o  o  o  fg  ooo 


SOURCE  LISTING  Of  3IUWQ 


PROGRAM  WQMAIN (OUTPUT. TAPER =0UTPUT,TAPE9 .TAPE  33) 

COMMON /NAMES /TITLE (d>  .CNAME (20) 

COMMON /CMTRL/NR , N  T  P  ,  N  W  Q  C , NB , NIT , NIE , IT ( 1 5 ) , IE ( 1 5 ) ,  .(CODE, 

♦  MCOOE ( 8 ) . I R . 13 , IPP , NCWQC , NNCWQC . NTWQC.NPP 
COMMON /IWQC/TWQ( 15,20), EWQ ( 1 5 , 20 ) , TQ ( 1 5 ) , EQ ( 1 5 ) . 

♦  8WQ ( 5 , 20 ) , BO ( 5 ) 

COMMON /EDA  TA/K20(  20, 20)  ,KT(  2  0  ) /QD  A  T  A/3  1  QiC  (  20  )  ,Q1 03(20) 

COMMON /SO ATA/320( 2 0,20) ,ST(20) 

COMMON /MDATA/LR( 20) . DA ( 20 ) . DEPTH ( 20 ) , VEL ( 20 ) , BCODE( 5 ) , 

+  HAV(20) ,H8V(20) , H AO ( 20 ) , HBD ( 20 ) 

COMMON/BDATA/PN , P? , A ( 20 ) 

COMMON /A  PA  RAM/ ALPHA (20) .BETA (20, 5) .GAMMA (20, 5) .DELTA (20) , 

♦  EPSIL(20,5) , ZETA ( 20 ) 

DIMENSION  HQ(200,20) , 0(200) , X(200, 3) 

REAL  K20 , KT , LR 

INTEGER  TITLE, HCODE , CNAME , D INDEX ( 125) 

DATA  Q10K/20*1.0/,Q10S/20»1.0/.PP/0.05/,PN/0.10/ 

DATA  (CNAME  (I),I*1,8)/"TEMP."  ,  "TSS"  , "BOD" , "NH3" , "N02" , "N03" , 

♦  "POP" , "D .0. " / 

->  INPUT  DATA  FROM  'TAPE9  ' 

CALL  0PENMS(9, DINDEX, 125,0) 

CALL  LOADDAT ( ITP ) 

->  BEGIN  SIMULATIONS  FOR  DISCRETE  TIME  PERIODS 


DO  100  ITP=1,HTP 
CALL  LOADTVP ( ITP ) 

- >  WATER  QUALITY  SIMULATIONS 


1  00 
C 


200 


CALL  SIMWQ ( WQ , Q , X ) 

->  OUTPUT  RESULTS  TO  'WQFLS ' 

CALL  WQOUTMS(WQ,Q,X,ITP) 
CONTINUE 


PRINT*,"  » 

SIMULATED  WATER  QUALITY  PROFILES  HAVE 
OUTPUT  TO  'TAPE33'  t'0R  ".IT?-’." 
PERIOD(S)  FOR  THE  FOLLOWING  WATER 
QUALITY  ATTRIBUTES : " 


PRINT*, " 

PRINT* , » 

PRINT*, " 

PRINT* , " 

PRINT*,"  " 

DO  200  I  ATT  *  1 , NTWQC 
PRINT*,"  "  ,  I  ATT , " ) 

PRINT*,"  " 

PRINT*, "*♦*♦*♦»♦*♦*♦*♦=* 
CALL  CLOSMS ( 9 ) 

STOP 


BEEN” 

*  ?  «  «  C*  M 


,  CNAME ( I ATT ) 
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INDEX  f OR  ALL  VARIABLES  AND 


ITERS  USED  IN  '31.TWQ' 


A  (  I  R  ) 

ALPHA ( IC ) 
ATT 

3C0DEC IS ) 
3ETA( IC , I ) 
3Q( 13 ) 

3WQ( 13 , IC  ) 
CNAMSUC) 
CNTRL3 ( 3 ) 


DA (  IR ) 
DELTA( IC, I) 
3  E  ?  T  H  ( IR) 
DINDEXU) 
EPSILi IC , I) 
EQ( IE) 

E WQC IE . IC ) 
GAMMA  (  IC  .  I  ) 
H  A  D  (  I  it  ) 

H  A  V  (  I R  ) 

H3D( IR  ) 

H8V ( IS ) 

HCODE 

IATT 

IB 

IC 

IE(IIS) 

INUM 

IP 

IPP 

IR 

I R  EC 

IRECIP 

IT(IIT) 

ITO? 

ITP 
I WQC 
KT 
K20 

LR( IR  ) 
MCODE( IC  ) 
NAATTS 
NB 

NCWQC 

NIE 

NIT 

NX 

NNCWQC 

NPP 

NR 

NS 

MTP 

MTWQC 

NWQC 

PN 

PP 


ALGAL  CONCENTRATIONS  REACH  'IR'  (tl'3/L) 

PARAMETER  rOH  ANALYTI-.  *  WQ  SOLUTION  (MG/L) 

KEY  FOR  MASS  STORAGE  OF  WQ  VECTOR 

BIFUR ACTION  CODE  DEFINING  WATERSHED  BRANCHES 

PARAMETER  FOR  ANALYTICAL  NO  SOLUTION  (MG/L) 

DISCHARGE  AT  MINOR  TRIBUTARY  AT  I3'TH 
BIFURCATION  <FT*«3 /SEC) 

WQ  VECTOR  AT  MINOR  TRIBUTARY  AT  IB'TH  BIFURCATION 
ALPHABETIC  NAMES  OF  WATER  QUALITY  ATTIBUTE3 
ARRAY  HOLDING  SIMULATION  CONTROL  PARAMETERS  FOR 
TRANSFER  TO  AND  FROM  MASS  STORAGE 
(NR, NTP , N  TWQC , ETC .  ) 

DRAINAGE  AREA  UPSTREAM  FROM  IR'TH  REACH  (MILES*»2) 
PARAMETER  FOR  ANALYTICAL  WQ  SOLUTION 
MEAN  DEPTH  OF  IR'TH  REACH  (FEET) 

MASTER  INDEX  FOR  MASS  STORAGE  OF  SIMULATION  DATA 
PARAMETER  FOR  ANALYTICAL  WQ  SOLUTION 

DISCHARGE  AT  IE'TH  POINT  SOURCE  EFFLUENT  INPUT  (FT*3/3EC) 
WATER  QUALITY  VECTOR  AT  IE'TH  EFFLUENT  INPUT 
PARAMETER  FOR  ANALYTICAL  WQ  30LTI0NS 
HYDRAULIC  RATING  PARAMETER  FOR  MEAN  DEPTHS 
HYDRAULIC  RATING  PARAMETER  FOR  MEAN  VELOCITIES 
HYDRAULIC  RATING  PARAMETER  FOR  MEAN  DEPTHS 
HYDRAULIC  RATING  PARAMETER  FOR  MEAN  VELOCITIES 
SWITCH  FOR  SPECIFYING  MODELING  CHOICE  FOR  DEPTH  VEL 
ATTRIBUTE  INDEX  (I  TO  20) 

-IFURACTION  INDEX  (1  TO  5 ) 

WQ  CONSTITUENT  INDEX  (1  TO  20);  SAME  AS  IATT 
EFFLUENT  INPUT  INDEX  (1  TO  15) 

INDEX  FOR  SPECIFICATION  OF  RECORD  KEY  FOR  MASS  STORAGE 
OF  WQ  VECTORS 
SAME  AS  IPP 

INDEX  Of  PROFILE  POINT  IN  SIMULATION  RESULTS  (1  TO  200) 
INDEX  OF  REACH  (1  TO  20) 

RECORD  KEY  INDEX  FOR  MASS  STORAGE  OUTPUT  OF  WQ  VECTOR 
INDEX  OF  RECIEVING  REACH  (DOWNSTREAM)  OF  IB'TH  BIFUCAIION 
INDEX  OF  INPUT  LOCATION  OF  TRIBUTARIES 
INDEX  SPECIFYING  TOP  OF  SUPERIOR  CHANNEL  BRANCHES 
INDEX  OF  TIME  PERIOD  (1  TO  12) 

INDEX  OF  WATER  QUALITY  CONSTITUENTS  (1  TO  20) 

TEMPERATURE  ADJUSTED  KINETIC  RATE  COEFFICIENT  (1/DAY) 
KINETIC  RATE  COEFFICIENT  AT  20  DEGREES  C  (1/DAY) 

LENGTH  OF  IR'TH  REACH  (MILES) 

MODELING  CODE  FOR  IC'TH  ATTRIBUTE 
NUMBER  OF  WQ  ATTRIBUTES 

NUMBER  OF  BIFURCATIONS  OF  STREAM  CHANNEL  (0<NB<5) 

NUMBER  OF  CONSERVATION  WQ  CONSITUENTS  ( 0 <NC WQC < 1 2- N NC WQC ) 
NUMBER  OF  POINT  SOURCE  EFFLUENT  INPUTS  (0<NIE<15) 

NUMBER  OF  TRIBUTARY  INPUTS  { 1<MIT<15) 

NUMBER  OF  RATE  COEFFICIENTS  NECESSARY  FOR  SIMULATION 
NUMBER  OF  NON-CONSERVATION  WQ  CONSTITUENTS 
NUMBER  OF  PROFILE  POINTS  IN  SIMULATION  OF  CURRENT 
NUMBER  OF  REACHES  IN  CURRENT  SIMULATION  (1<NH<20) 

NUMBER  OF  SOURCE/SINK  TERMS  IN  SIMULATION  (0<TNS<20) 

NUMBER  OF  TIME  PERIODS  IN  SIMULATION  RUN  (  1 <N TP < 1 2 ) 

NUMBER  OF  TOTAL  WQ  CONSTITUENTS  IN  SIMULATION  (NTWQC<20) 
NUMBER  OF  CORE  WATER  QUALTIY  CONSTITUENTS  IN  SIMULATION 
PERCENT  OF  NITROGEN  IN  SUSPENDED  SOLIDS  (*/100> 

PERCENT  OF  PHOSPHORUS  IN  SUSPSENDED  SOLIDS  ( * / 1 0 0 ) 
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o  o  o 


Q(I?P)  STREAM  DISCHARGE  AT  IPP'TH  POINT  OF  OUTPUT  PROFILE  (FT»»3/SEC 

OB  DISCHARGE  INPUT  FROM  BIFURCATION  TO  CURRENT  REACH  (CFS) 

OE  DISCHARGE  INPUT  FROM  TRIBUTARIES  TO  CURRENT  REACH  (CFS) 

GINDSX(I)  MASTER  INDEX  FOR  MASS  STORAGE  OF  SIMULATED  W Q  PROFILES 

OT  DISCHARGE  INPUT  FROM  TRIBUTARIES  TO  CURRENT  REACH  (CFS) 

QO  DISCHARGE  AT  TOP  OF  CURRENT  REACH  AFTER  DILUTIONS  (CFS) 

QJOXdi C)  TEMPERATURE  CORRECTION  FACTOR  FOR  IK'TH  RATE  COEFFICIENT 

QIOS(IS)  TEMPERATURE  CORRECTION  FACTOR  FOR  IS 'TH  SOURCE/SINK  TERM 

ST(IR.IS)  TEMPERATURE  CORRECTED  SOURCE/SINK  TERM  FOR  I R  ' TH  REACH 

520(IR,I5)  SOURCE/SINK  TERM  AT  20  DECREES  C  (MG/L/DA X) 

TEMP  WATER  TEMPERATURE  (DEGREES  C) 

T I T L E ( 8 )  SIMULATION  NAME 

TLR  TOTAL  LENGTH  OF  STREAM  TO  BE  SIMULATED  (MILES) 

TQ(IT)  DISCHARGE  IN  IT'TH  TRIBUTARY  (CFS) 

TTWR  TOTAL  TIME  OF  TRAVEL  WITHIN  CURRENT  REACH  (DAYS) 

TWQ(IT,IC)  WATER  QUALITY  VECTOR  FOR  IT'TH  TRIBUTARY  INPUT 

VEL(IR)  MEAN  VELOCITY  IN  IR'TH  REACH  (MILES/DAY) 

WQdPP)  WATER  QUALITY  VECTOR  AT  IPP'TH  POINT  IN  SIMULATION 

WQBdB.I C)  WATER  QUALITY  VECTOR  FOR  IB'TH  BIFURCATION  INPUT 

WQEdS.IC)  WATER  QUALITY  VECTOR  FOR  IE 'TH  EFFLUENT  INPUT  1 

WQT(IT.IC)  -  WATER  QJALITY  VECTOR  FOR  IT'TH  TRIBUTARY  INPUT 

WQO(IC)  WATER  QUALITY  VECTOR  AT  TOP  OF  CURRENT  REACH  AFTER  < 

XdPP.3)  DISTANCE  ARRAY  SPECIFING  DISTANCE  DOWNSTREAM,  REACH  ' 

XIB  INDEX  OF  LOCATION  ON  WATERSHED  BRANCHES  ' 

XINC  INCREMENT  BETWEEN  PROFILE  POINTS  IN  SIMULATION  OUTPUT  < 

XIR  INDEX  OF  LOCATION  IN  TERMS  OF  STREAM  REACH  < 

XLAST  DISTANCE  AT  TOP  OF  CURRENT  REACH  ' 

XWR  DISTANCE  FROM  TOP  OF  CURRENT  REACH  « 

ZETA  PARAMETER  FOR  ANALYTICAL  WQ  SOLUTION  I 

< 

END  1 


SUBROUTINE  ASOLN ( TTWR , tfQ ) 

. . . . 

COMMON/CNTRL/NR.NTP.NWQC ,NB,NIT,NIE, IT ( 1 5 ) , IE ( 1 5 ) , HCODE , 

+  MCODE ( 8 ) , IR, IB, I P P , NCWQC , NNCWQC , NTWQC.NPP 
COMMON /APAR AM /ALPHA (20), BETA (20, 5). GAMMA (20, 5), DELTA (20), 

♦  EPSIL(20,5),2ETA(20) 

DIMENSION  NE(20) , WQ ( 200 , 205 

DATA  (NE(I),I*1,20)/3*',2,2*3,2,5,12*1/ 

D0SAT(TEMP)=1«. 652-0.  A  1 022»TEMP*0 . 00799  !  »TEMP  *  »2 

*  -0.0000777A»TEMP*»3 

DOWNSTREAM  PROFILE  CALCULATION  OF  WATER  QUALITY 
DO  899  J *  1 , NTWQC 

WQdPP, J)»ALPHA( J )*DELTA( J)*TTWR *ZETA(J) »TTWR»«2 
IIsNE(J) 

DO  799  K=.1,II 

WQ(IPP, J) »WQ(IPP, J)*BETA( J,K)»EX P( -GAMMA (J,K)*TTWR) 
WQdPP, J)*WQ( I PP, J)*EP3IL( J,K )«TTWR* EXP (-GAMMA ( J,K) »TTWR ) 
799  CONTINUE 
899  CONTINUE 

WQdPP,8)rOOSAT(  WQdPP,  1  )  ) -WQ(IPP,  8) 

RETURN 
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SUBROUTINE  DILUT£(WQ,WQ0,QG) 

common/cntrl/nh,  ntp  , nwqc , mb, nit, nie, i~( 15) . isl 15) . hcooe, 

♦  :ICODE(  d) ,  IR,  13.  IPP.NCUQC.NNCWQC,  NTWQC.  WPP 
COMMON / 1  HOC /TWO ( 15,20),EWQ(15.20).TQ(l5)tEQ(l$). 

♦  BNG(5,20),9Q(5) 

C0MMGN/MDArA/LR(20) ,DA(20),DEPTH(20),VEL(20),BCODE(5), 

♦  HAV(20 ) , HBV(20 ) , HAD< 20) , H8Q(20) 

DIMENSION  WOE ( 20 ) ,UQT(20) , WQB I  20 ) , HQ( 200 , 20 ! , HQ0C20 ) ,00(20) 
REAL  LR 


C - >  ZERO  TRIBUTARY  AND  EFFLUENT  INPUTS 

C 

101  QT  =  0 . 0 
QE  =  0  . 0 
QBsO  . 0 

DO  102  Is  1  ,  N  TWQC 
HQ3(I)sO.O 
HQE( I ) sO  .  0 

102  W  Q  T  ( I  )  s  0 . 0 
C 

c - >  SELECT  TRIBUTARY  INPUTS  FOR  CURRENT  REACH 

C 

DO  110  Is  1 , N IT 

IF ( I T ( I ) . NE. IR )  GO  TO  110 

QTsTQ(I) 

DO  105  Jsl.NTWOC 
105  WQT( J)*TWQ(I, J) 

GO  TO  111 

110  CONTINUE 
C 

o - >  SELECT  POINT  SOURCE  EFFLUENT  INPUTS  FOR  CURRENT  REACH 

C 

111  DO  120  Is  1, HIE 

IF ( I E ( I ) .NE.IR)  GO  TO  120 
OEsEQ(I) 

DO  115  J.'  .  MTNOC 
115  NQE( J)=£WQ( I, J) 

GO  TO  121 

120  CONTINUE 
C 

c - >  SELECT  INPUTS  FROM  BIFURCATIONS 

C 

121  DO  130  Isl.NB 
iRECIPsIHT(BCODEd)  ) 

IFdRECIP. NE.IR)  GO  TO  130 
QB=3Q( I ) 

DO  125  Jsl.NTWQC 
125  WQB ( J ) sBWQ ( I , J ) 

GO  TO  131 

130  CONTINUE 
C 

C - >  DILUTE  HATER  QUALITY  VARIABLES 

C 

131  IF(IR . EQ .  1 )  GO  TO  1  NO 
DO  135  Is  1 , NB 

ITO?s3COOE< I )*100.-INT(BCQDE£I))*!00*1 
IF(ITOP.EQ.IR)  GO  TO  140 
135  CONTINUE 

GO  TO  150 
1  A  0  Q0( IR ) sQE*QT 

DO  145  Ks I, NTWQC 

145  WQ0(K>=(QE»WQE(K ) *QT • WQT ( K )  ) /QU ( IR ) 

GO  TO  2  10 

150  DO  200  Ksl.NTWQC 

200  WQO  (  iC  )  *  (  QE  *HQ£  (  K  )  *QT  •WQT(  K  )  *QB  *WQB  (  K  )  *Q0(  IR  -  1 )» 

*  WQ( IPP-1 . X) )/Q0( IR ) 

210  CONTINUE 
RETURN 
END 
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SUBROUTINE  LOADDAT ( ITP  ) 

HEAL  LR.K20 

CDMMON/NAMES/TITLEC  d ) , CNAME ( £0 ) 

COMHON/CNTRL/NR,  NTP  ,  UUQC  ,  NB  ,  HI  T  .  NIE  ,  IT  (  1  5  )  .  IS  (  !  5  )  .  HCODE  , 

♦  HCODE  (  8  ),IR,I3,IPP,  MC1J3C,  MNCWGC,  HTKQC,  NPP 
COMMON /I XQC/TWQ( 1 5 , 20 > , EWQ ( 1 5 , 20 ) , TQ ( 1 5 ) , EQ ( 1 5 ) , 

♦  a«G(5.20),BQ(5) 

COMMON /KDATA/K20( 20, 20) , XT ( 20 ) /GDA TA / Q 1  OK  (  2 0 ) ,0103(  20) 
COMMON/S  DA  TA/S20(20, 20) , ST (20) 

COMMON/MDATA/LR( 20) , DA ( 20 ) , DEPTH ( 20 ) , VEL ( 20 ) , BCODE ( 5 ) , 

♦  HAV( 20) , HBV ( 20 ) , HAD ( 20 ) , HBD( 20) 

COMMON /BDATA/PN , PP, A ( 20) 

INTEGER  TITLE , HCODE , CNAME , DINDEX ( 1 2 5 ) , CNTRLS ( 3 ) 

C - >  UNLOAD  DATA  PROM  MASS  STORAGE  FILE 

CALL  RE ADMS (9, TITLE, 8,  1  ) 

CALL  READMS( 9 .CNTRLS, 3, 2 ) 

NR*  CN  TRLS (  1 ) 

NTP=CNTRLS( 2 ) 

NCWQC*CNTRLS(3) 

NNC’XQC  =  CNTRL5(  4  ) 

HC0DE*CNTRLS(5 ) 

NB  =  CN TRLS ( 6  ) 

NI?*CNTRLS(7) 

NIS*CMTRLS(8 ) 

NWQC  » 3 

NX* 1 6*NNCWQC 
NS*3*NNCWQC*NCWQC 
NTWQC*8*NCWQC*NNCWQC 
CALL  READMS(9,MCODE,8,3) 

UAATTSsNCWQC+NNCWQC 
IF(tlAATTS.LE.O)  GO  TO  570 
CALL  READMS ( 9 , CNAME , NTWGC , 4 ) 

570  CALL  R£ADMS(9,8CODE,NB,5) 

CALL  READMS(9,IT,Nir,6) 

CALL  REA0MS(9,IE,NIE. 7) 

CALL  R£ADMS<9,LR.NR,8) 

CALL  REAOHS ( 9 , DA , NR , 9 ) 

IF(HCOOS.EQ.O)  GO  TO  580 
CALL  R£ADMS(9,HAV,NR, 10) 

CALL  READMS<  9,HBV,NR,  1 1  ) 

CALL  READMS ( 9 ■ HAD , N R . 12) 

CALL  READMS ( 9 , HBD , NR , 13) 

580  RETURN 

ENTRX  LOADTVP 
IIs13*(ITP-1)*9 
IF (HCODE . EQ . 1 )  GO  TO  585 
....M*II*1 

CALL  READMS(9,VEL,NR,IHUM) 

INUM *11*2 

CALL  READMS(9, DEPTH, NR, INUM) 

585  INUM. 11*3 

CALL  READMS(9,TWQ, 300, INUM) 

INUM* I 1*4 

CALL  READMS(9,TQ, NIT, INUM) 

INUM* 1 1*5 

CALL  READMS  (  9'.  EWQ,  300  ,  INUM) 

INUM *11*6 

CALL  READMS ( 9 , EQ , N IE , INUM) 

INUM* 1 1*7 

CALL  READMS(9,K2l), 400, INUM) 

INUM. 11*8 

CALL  READMS(9.S20, 400 , INUM) 

INUMsII+9 

CALL  REA  DM3 (9, A, NR, INUM) 

CONTINUE 

RETURN 

END 
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subroutine  pcalc;nqo) 


C3MM0N/APARAM/ALPHA ( 20 ) . 3ETA ( 20 , 5 ) , GAMMAl 20 . 5 ) . DELTA( 20 ) , 

*  EPSIL ( 20 , 5 ) , 2£TA( 20 ) 

COMMON /CNT3L/NH  ,  MTP,  N’rfQC  ,  M3  ,  NIT.NIS,IT(15),I£(15>.  HCODE  , 

*  HCODE ( d ) , IS , IB , IPP . NCtfQC , NNC'NQC , NTafQC , NPP 
C0MMON/.<DATA/K20(  20, 20)  ,  XT( 20 ) /QDATA/Q  10X(20),Q10S(20) 
COMMON/SDATA/S20( 20 , 20 ) , 3T ( 20 ) 

COMMON /3DAT A / PN  ,  PP  ,  A(  20) 

COMMON /MDATA/LR (20),DA(20).DEPTH(20),VEL (20), 3C0D£( 5 ) . 

*  HAV( 20) , HBY( 20) ,HAD( 20) ,H3D( 20) 

DIMENSION  200(20) 

BEAL  XT, K20, LR 

DATA  (DELTA ( I ) , 1= 1 ,20 ) /20»0 . 0/ 

DATA  t(EPSIL(I,J),I=1,20),J=1,5)/100«0.0/ 

DATA  (ZETA( I ) , 1=1 ,20) /20*0. 0/ 

DGSAT(rEMP)  =  14.o52-0.4  1G22#TEMP-*0. 007 9 9 1  ‘TEMP  * *2 
+  -0.00007774»TEM?'*3 
C 

C----  TEMPEBAT'JBE  PARAMETERS 

C  - 

100  GAMMA (  1 .  1  )  =  XT ( 15) 

IF  ( GAMMA (1,’).SQ.O.O)  GO  TO  150 
ALPHA ( 1 ) *ST( 8 ) /GAMMA (1,1) 

8ETA( 1 , 1 )=WQO( 1 ) -ALPHA ( 1 ) 

GO  TO  200 

150  DELTA( 1 ) =ST( 8) 

ALPHA ( 1 ) =  NQO( 1 ) 

BETA ( 1 , 1 >=0. 0 
C 

C-  — -  BOD  PARAMETERS 
C  - 

200  GAMMA12, 1 ) = XT ( 4 ) .XT ( 5 ) 

IF (GAMMA (2,  1  ) .£0.0.0)  GO  TO  220 
AL?HA(2)=ST( 1 )/GAMMA(2, 1 ) 

B  ET  A ( 2 ,  1 ) =WQ0(2)-ALPHA(2 ) 

GO  TO  300 

220  DELTA(2)=ST( 1 ) 

ALPHA(2)=WqO(2) 

BETA ( 2 , 1 )=0.0 
C 

C----  TOTAL  SUSPENDED  SOLIDS  PA  AMETERS 

C  . . 

300  GAMMA ( 3 ,  1  ) =  XT ( 1 3 ) 

Ir"(0E?Trf<  IR)  .GT.  0 . 0  )  GAMMA  (  3,  1  )  =  GAMMA  (  3,  ’  )  *XT(  1  4  )  /  DE  PTH  (  I  R  ) 
IF ( GAMMA ( 3 ,  1  )  • EQ- 0. 0  )  GO  TO  340 
310  ALPHA(3)-ST(7)/GAMMA(3,  ') 

B  ET  A ( 3,1) =WQO( 3)-ALPHA( 3) 

GO  TO  400 

340  DELTA(3)=ST(7) 

ALPHA ( 3 )  =  WQO ( 3 ) 

BETA ( 3 , 1 )=0.0 
C 

C----  AMMONIA  PARAMETERS 
C  - 

400  GAMMA ( 4, 1 )=GAMMA( 3, 1 ) 

GAMMA (4.2 ) =XT( 7 ) 

IF (GAMMA ( 4 ,  1  )  . EQ . 0 . 0 )  GO  TO  430 
(GAMMA ( 4 , 2 ) . EQ . 0 . 0 )  GO  TO  460 
ALPHA ( 4 )  =  (KT(0 ) *PN ‘ALPHA ( 3 )-3T(2  >-XT( 1  ) «A(  IR  )  ) /GAMMA ( 4 ,2) 

IF ( GAMMA ( 4 , 1 ) . EQ . GAMMA (4,2))  GO  TO  410 

BETA  (  4  ,  1  )  =  (K7(  6)  »PN»BETA  (  4  ,  1  )  .'/(GAMMA  (4,2)  -GAMMA  (4,1)) 

GO  TO  490 

410  BETA ( 4 , 1 ) =0 . 0 

EPSIL(4,')=XT(6)»PN»ALPHA(3) 

GO  TO  490 
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430  IF  (GAMMA  l  4  ,  1  )  .  tu  .  oAMMA  l  4 ,  <L  I  )  GO  1U  hou 

ALPHA!  4  )  a  (  KT!o)*P;i*ALPHA(  3)*>T!2)-Kf(  1  )  *A  (  IR)  )  /GAMMA  (4,2) 
BETA ( 4 ,  1  )  :0 . 0 
ZETA(4)*<T(0)»PN«ST(7)/2 
GO  TO  490 

460  BETA  (  4  ,  l)  =  i!T(  6  )  »PN» ALPHA!  3) /GAMMA (4  ,  1  ) 

ALPHA ( 4 ) =0. 0 

DELTA!  4  )  =oT  <  2  )-KT(  1  )  «A(  IR  )  «-  ( iC  T  (  6  )  »PN»S:(  7  )  ) /GAMMA! 3,  I  ) 

00  TO  500 

460  DELTA ( 4 ) =ST( 2 ) 

ALPHA ( 4 ) =WQO( 4 ) 
a£TA(4, 1 )=0.0 
oETA(4,2)=0.0 
GO  TO  500 

4  90  BETA! 4 , 2 ) =  WQ0( 4 ) -ALPH  A( 4 ) -BETA (4 , 1  ) 

C 

C----  NITRITE  PARAMETERS 

C  - 

500  GAMMA ( 5, 1 ) .GAMMA ( 3, 1 ) 

GAMMA(5,2)=GAMMA(4.2) 

GAMMA ( 5 . 3 ) s  KT ( 6 ) 

IF ( GAMMA ( 5 . 3 ) . FQ . 0 . 0 )  GAMMA! 5, 3)  = . 00001 
AL?HA(5)=GAMMA(4.2)»ALPHA, 4 ) /GAMMA (5,3) 

IF  ( GAMMA  (  3  .  3  )  .  E<3  .  GAMMA  (  5  ,  )))  Go  TO  510 

BETA!  5  ,  1  )  .GAMMA!  4, 2  )  »BETA(  4  ,  1  )  /  (  GAMMA!  5,  3  )  -CAMflA!  5  ,  1  )  ) 

502  IF (GAMMA! 5 , 3 ) . £Q . GAMMA ( 5 , 2) )  GO  TO  520 

BETA! 5 , 2 ) .GAMMA (5,2)*BETA(4,2) / (GAMMA ( 5 , 3 ) -GAMMA (5,2)) 

506  BETA ( 5 , 3)=WQ0(5)-ALPHA(5)-BETA(5. 1 ) -BETA (5,2) 

GO  TO  600 

510  BETA ( 5  ,  1  )  =  0.0 

EPS IL ( 5 , ' ) .GAMMA (4,2)*BETA(4,1) 

GO  TO  502 

520  B£TA(5,2)=0.0 

EPSIL(5,2).GAMMA(5,2)#SETA!4,2) 

GO  TO  506 

C 

C----  NITRATE  PARAMETERS 


600  GAMMA! 6, 1 ) .GAMMA! 3, 1 ) 

IF (GAMMA! 6, 1 ) . EQ . 0 . 0 )  GAMMA! 6, 1 )=. 00001 
GAMMA (6, 2) .GAMMA (4, 2) 

IF ( GAMMA (6,2).EQ.O.O)  GAMMA ( 6 , 2 ) = . 0000 1 
GAMMA (6,3) .GAMMA ( 5 , 3) 

ALPHA(6).WQ0(6) 

BETA (6, 1 )= GAMMA! 5, 3 ) ‘BETA! 5, 1  ) /GAMMA! 6, 1 ) 

BETA (6 , 2 ) .GAMMA (5 , 3 ) #BETA( 5 , 2 ) /GAMMA (6,2) 
BETA(6,3)=BETA(5,3) 

DELTA(6).KT(a)»ALPHA(5)*ST(3)-KT(2)»A(IR) 

C 

c----  PHOSPHATE  parameters 

c  - 

700  GAMMA(7, 1 ). GAMMA! 3, 1 ) 

GAMMA ( 7 , 2 ) »KT ( 1 1 ) 

IF(GAMMA(7, 1 ) .EQ.0.0)  GO  TO  730 
IF(GAMMA(7,2) .EQ.0.0)  GO  TO  760 

ALPHA { 7 ) = ( XT ( 10) •??• ALPHA! 3 ) ♦ST ( 5 > -KT ( 3 ) *A ( I R  )  ) / KT ( 1  1  ) 
IF (GAMMA ( 7 , 1 ) . EQ • GAMMA ( 7,2))  GO  TO  710 

BETA!  7,  1  ).(KT'<  10)  »PP»BETA!  3.  1  )  )  /  (  GAMMA  (  7 , 2  )  -GAMMA  (  3,  1  )  ) 
GO  TO  790 

710  BETA ( 7,11=0.0 

£?3Ik-v7,O.KIv10/*  F*AuPHA  v  3/ 

GO  TO  790 

730  IF (GAMMA ( 7 , 1 ' • EQ . GAMMA! 7,2) )  GO  TO  760 

ALPHA(7).(KT( 1 0 ) • P P • A L PH A ( 3 ) *3T ( 5 ) -K T t 3 ) * A ! I R  )  ) /KT (  1  1) 
BETA ( 7, 1 ) =0 . 0 
2ETA(7)=KT( 10)*PP*ST(7)/2 
GO  TO  790 


/  u  u 


7  30 


790 

C 

C - 

L 

aoo 


302 
309 
80  a 
808 


810 


320 


330 


8«0 


C 

c - 

c 

900 


930 

C 

c - 

c 

950 


980 


970 

C 

C 

C 

999 


»  £,  4  H  i  <  •  •  ,/SfW  v  '  '>  J  “rr-*wfnnv  3  /  '  »J  * .  i  >•»  *  v  i  ,  *  y 

9  l  ?  H  A  (  n  5  o . 

DELTA (  n =8  7(  )  )-< 1 (  j  i  »A  C IR ) *(K  T(  10 ) *PP»ST( 7 )  ) /GAMMA (  1 .  1) 
00  70  JOO 

del  "a  .  7  •  =ot;  'i 

A  L  P  ii  A  (  '  )  -  0  0  J  \  7  i 
BETA  C  7 .  1  )  =  0  .  J 
BETA l 7 , 2 ) =0 . 0 
00  TO  300 

BETA l 7,2)=WQ0( 7)-ALPHA(7)-3ETA(7. 1 ) 

DISSOLVED  OXTGEN  PARAMETERS 


G  A  M  M  A  (  3  ,  1  )  =GA!1MA  (  3  ,  1  > 

GAMMA(8,2)=GAMMA(5, 3) 

GAMMA (3,3 ) s GAMMA (9,2) 

GAMMA (3,9) sGAMMA (2,1) 

GAMMA ( 3 , 5 ) =KT( 12) 

ALPHA(8)  =  (<T(9)»ALPHA(2)*3-93»KT( 7 >  »A LPHA ( 9 ) *  1 .  1 9 »K T ( 3 ) »AL PH A ( 5 ) - 
+  ST(5)*ST(6))/KT( 12) 

If (GAMMA ( 3 , 5 ) . EQ . GAMMA ( 3 , 1 ) )  GO  TO  810 

3ETA( 3, 1 ) *( 3 • 93*3AMMA( 9 , 2 ) »8ETA( 9 , 1 )*1 .  1 9»GAMMA ( 5 . 3 ) »B ET A ( 5 ,  1  )  )  / 

*  (GAMMA (8,5) -GAMMA ( 3  t ' ) ) 

IF (GAMMA ( 8 , 5 ) ■ £3 • GAMMA (5,3))  GO  TO  320 

BETA ( 8 , 2  )  =  (  1  .  1 9»GAMMA( 5, 3) ‘BETA (5 , 3) )/( GAMMA (8,5) -GAMMA (5,3)) 

IF ( GAMMA ( 3 , 5 ) - EQ • GAMMA (9,2))  GO  TO  330 

BETA ( 8 , 3  )  =  ( 3 . 93»GAMMA( 9 , 2 ) »3ETA (5 , 2 ) > / ( GAMMA (3,5) -GAMMA (9,2)) 

IF (GAMMA (8, 5) .£Q.GAMMA(2, 1 ) )  GO  TO  890 

BETA ( 8 , 9 ) s (GAMMA ( 2 , 1 ) -KT ( 5 ) ) *3ET A ( 2 , 1 )/ (GAMMA ( 8, 5 ) -GAMMA (2,1)) 
BETA ( 8 , 5 ) = ( DOSAT (  WQO  ( 1 ) )-WQ0( 8 ) ) -ALPHA ( 3 ) -BETA ( 8 , 1 ) -BETA ( 8 , 2 ) 

♦  -BETA ( 8 , 3 )-0ETA( 8 . 9 ) 

GO  TO  900 

BETA ( 3 , 1 )=0. 0 

EPSIL (8,1 )=3 . 93»GAMMA(9 , 2) ‘BETA ( 9, 1 ) ♦ 1 . 19»GAMMA( 5, 3)*BETA(5, 1 ) 

GO  TO  809 
BET  A ( 8 , 2 ) s  0 . 0 

EPSIL (8,2) si  .  19»GAMMA(5, 3) ‘BETA (5, 1 ) 

GO  TO  306 
BETA ( 3 , 3>s0.0 

epsil ( 8, 3 )  =  3- 9 3 *gamma < « , 2) »beta  ( u ,  i  ) 

GO  TO  303 
BETA ( 3, 9  )  =0 . 0 

EPSIL(8,9)*uAMMA(2, 1 ) -KT ( 5 > *8ET A ( 2 , 1 ) 

GO  TO  308 

CONSERVATIVE  WQ  CONSTITUENTS 


IF ( NCW5C . LT . 1 )  GO  TO  950 
DO  930  1=1, NCWQC 
ALPHA ( 8*1 )s'WQO(  3*1 ) 

BETA ( 8*1 ,  1 ) =  0 , 0 
GAMMA ( 8*1 ,  1  )  =  0.0 
0ELTA(8+I)=ST(8*I) 

NONCONSERVATIVE  HQ  CONSTITUENTS 


IF ( NNCWQC . LT . 1 )  GO  TO  999 
DO  970  Isl.NNCHQC 
GAMMA(3*NCWQC*I, 1 ) =KT( 8*NCWQC*I ) 

IF (GAMMA ( 3*NCH3C*I , 1 ). EQ . 0 . 0 )  GO  TO  960 

ALPHA ( 8*NCWQC*I )= ST ( 8* NCWQC* I ) /GAMMA ( 8*NCWQC*I , 1 ) 

BETA( 0*NCWQC*I ) = WQG ( 3 *MCWQC*I )- ALPHA ( 8*NCWQC.I ) 

GO  TO  970 

DELTA l J*NCWQC»I ) =3T ( 3*NCWQC ♦ I ) 

ALPHA ( 3*NCWQC*I ) *WQ0(8* NCWQC*!) 

BETA ( 3*NCWQC*I , 1 ) =0 . 0 
CONTINUE 


RETURN 

ENO 
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SUBROUTINE  RCALCCSMP) 


COMMON  /CiJTRL/NR  ,  NIP  ,  NWOC  ,  N  B  ,  .V I  P  ,  M I  £  ,  1 7  (  1  5  )  ,  1 2  (  15),  HCODE  , 
MCODE ( 3 ) , I R , IB , IPP ,  NCWQC  ,  N NCWQC  .  NTWQC  ,  HPP 
COMMON  /KOA  r  A /K20  (  2  0  ,  2  0  )  .  .<T  (  20)/QDA  PA /Q 1  OK ( 20 ) , Q 1  OS ( 2 0 ) 
COMMON /SDATA/S20(20,20),3T(20) 

REAL  <20. K T 


u 

C - >  Q10  CONVERSIONS  OF  RATE  CONSTANTS 

C 

NK  *  1  6«-NNCWQC 
NSsd»NC'NQC*NNCWqC 
DO  100  I s  1  , NK 

100  KT(I)sK20(IR,I)*Q10K(I)**(TEM?-2Q) 

DO  200  1= 1 ,NS 

200  ST(I)=S20(IR, I)*Q10S(I)»*(TEMP-20) 

C 

RETURN 

END 

£=--=232333333333333333333333333333333335333333333333333333333333333333333 

SUBROUTINE  S IMWQ ( WQ , Q , X ) 

COMMON /CMTRL /NR, NTP.NWQC, W3, WIT, NIE, IT( 15) , IE ( 15) ,HCODE , 


♦  MCODE  (  3)  ,  I  a  ,  13  ,  IPP.NCWQC, NNCWQC  ,  NTWQC  ,  NPP 

COMMON /MDA  TA / LR ( 20 ) , D A ( 20 ) , DEPTH ( 20 ) , VEL ( 20 ) , BCODE ( 5  ) , 

♦  HAV( 20) , HBVI 20) , HAD( 20) ,HBD! 20) 

COMMON /  IWQC /TWQ  (  15,20),  EWQ(  1  5 , 20  )  ,  TQ  (  1  5  )  ,  EQ<  1  5  )  , 

♦  BWQ ( 5,20) , BQ ( 5  ) 


DIMENSION  WQO I  20 ) , WQ( 200 , 20 ) , Q0( 20) , Q<200 ) , X(200, 3 ) 

REAL  LR 
INTEOER  HCODE 
C 

C - >  CALCULATE  DISTANCE  BETWEEN  PROFILE  POINTS  FOR  OUTPUT 

C 

TLRxO.O 
DO  10  IR.I.NR 
10  TLRxTLR  +  LRUR) 

XI NC  a  0 .  1 
DO  20  Is  1  , 5 
NPPsTLR/XINC*2»NR 
IF(NPP.LE.200)  00  TO  25 
20  XINCsI *0 . 5 

C 

C - >  CALCULATE  WATER  QUALITY  PROFILE  FOR  A  FIXED  TIME  PERIOD 

C 

2  3  Z ? ’ s  0 


X 1 3  s  N  3 
XLASTsO.O 

100  DO  500  I R s 1 , NR 
XIRsIR 
XWRsO.O 
IPPxIPP+1 
C 

CALL  DILUTE(WQ,WQO,QO) 

C 

DO  110  IWQCsI.NTWQC 
110  WQ(IPP,IWQC)sWQO(IWQC) 
Q(IPP)*QO(IH) 

KIPP,  1  )sXLA3T 
X(IPP,2)»XIR 
XlIPP, 3 ) 3  XIB 
C 

IF ( HCODE . LT . 1 )  00  TO  150 
VEL(IR)*HAV(IR)»QO(IR)»*HBV( IR> 
DEPTH ( IR)*HAD(IR)»QO(IR)«»HBD(IR) 
C 

150  CALL  RCALC ( WQO ( 1 )  ) 

C 
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CALL  PCALCiWGOI 


2  0  0  IPP=IPP*1 

XW3  =  XWR*XI:iO 

XrU'Jil.^.UdB))  GO  TO  300 
X{ IPP ,  1  )  =XLA3T*XWR 
X  (  I  ?  P  .  2  )  =  X  IP 
X ( IP? , 3 ) =  XIB 
3(  IPP) =G0( 13) 

TTWRsXWR/VELC 13  ) 

c 

CALL  ASOLN  t  TTWR , WQ ) 

C 

QdPP)sQO(ia) 

00  TO  2 00 
300  XWR  =  LRCR) 

XU??,  1  )rXLAST*XWR 
X( IPP ,2  )  =  XIR 
X(IPP, 3)=XIB 
Q(I??)=QO(IR) 
ttwh=xwr/vel(ir) 
c 

CALL  ASOLN  (  TT«f3  ,  i«Q  ) 

3 

IFCNB.LT. 1 .OR.XIB.EQ.O. )  GO  TO  500 
DO  POO  Isl.NB 

IEND=3C0DE( I  )»100.  -IUTC3C0DEC  I)  )  •  1  00  . 
IFCIEND.NE.IR)  GO  TO  «00 
XI3=XI3-1 
3QC I )  =3( IPP) 

DO  350  IC*1, NTWQC 
350  Bwj(i,rc)saQ(i?p,ic) 

LOO  CONTINUE 

500  XLAoTsXCIPP,  1  ) 

NPPsIPP 

RETURN 

END 

SUBROUTINE  WQ0UTM3 ( NO , Q , X , I T P  ) 

COMMON /CN THL/NR .NTP.NWQC, NB, NIT , NIG, ITi  15  /, IE 
*  MCOOEC  8  l.IR.IB.IPP.  NC’NQC  ,  NNCWQC  ,  NTWQC  ,  HP? 
INTEGER  CNAME , TITLE , rtCODE 
IFCITP.GT.  1  )  GO  TO  10 
CALL  OPEN  IIS  (33iO  INDEX, 277,0) 

10  IREC=23*(ITP-1 )-1 

CALL  WRITMSC 33, NPP, 1 , IREC) 

IREC=23*( IT?  -  1 )*2 

CALL  WH  X  TftS  C  3  3  .  X,  600,  IREC) 

IREC  =  23* (ITP-1 ) *  3 

CALL  WRITMSC 33, Q, NPP, IREC) 

DO  30  ICal.MTWQC 
DO  20  IP= 1 , NPP 
20  ATTUP)*WQ(IP,  IC) 

IRECs23*(ITP-1 >  •*■  I C* 3 
30  CALL  WRITMSC 33. ATT, NPP, IREC) 

IFCITP.LT. NTP)  GO  ro  too 
CALL  CL03M3C33) 

100  RETURN 

END 


::::::::: ii: :: c: :::: ii 


' ' . ATT  C 200  ) 
1  5  /  ,  N  C  G  0  E  , 
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PROGRAM  UQTEST (INPUT, OUTPUT, TAPE3  =  IN  PUT, TAPE4  =  0UTPUT, 

♦  TAPE9.TAPE33) 

COMMON/NAM£S/TITL£(8) , CNAME ( 20 ) 

COMMON /CNTRL /NR ,NTP,NWOC,NB,NIT,NIE,IT( 15) , IE ( 15) , 

♦  MC0DEC8) , NCWQC , NNCWOC , NTWQC , NPP 
COMMON /MOAT A /LR( 20) , OA ( 20  )  , 0CO DE ( 1 5 ) 

DIMENSION  X(200,3),WQ(240),PTS(2>,XWQSV(20, 1 2 ) , UMQS ( 20 , 1 2 ) , 

♦  LWQSC20, 12) , XX(240) , IATTS(20) 

REAL  LR.LHQS.LS 

INTEGER  TITLE, CNAME 

DATA  ( CNAME ( I ) , I s 1 , B ) /"TEMP . * , "BOOS* , "TSS" , "NH3" , "N02 » , "NO 3" 

♦  "P04"  , "D.O. "/ 

DATA  XWQSV/240*0.0/,UWQ3/240»9999./,LWQS/240*0.0/, I0PEN/0/ 
PRINT*, "  " 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  « 

PRINT*  ,  ”  +  3*3+3  +  2  +  3+  3-*S+3+X  +  34-3  +  3  +  3  +  X  +  3  +  2  +  3-»X  +  3+  3+3  +  S  +  Z+  34.3  +  34 
PRINT*, "♦  THIS  RTV  ROUTINE  TESTS  FOR  VIOLATIONS  OF  AMBIENT 
PRINT*, "♦  WATER  QUALITY  STANDARDS  AND  QUANTIFIES  THE 

PRINT*, "♦  SPACIAL  EXTENT  OF  THESE  VIOLATIONS 

PRnT*,l,tX*343«3«3aa*333+S*3  +  3«l4a  +  3  +  S«343*3  +  3>S+3>84S+3»3*}4 
PRINT*,"  " 

C - >  READ  DATA  FROM  'TAPE9  ' 

CALL  DATINMS 

C - >  DESIGNATE  ATTRIBUTES  OF  INTEREST 

PRINT*, " - 

PRINT*,"  ",  (TITLE(I) ,l3l ,8) 

PRINT*, " - 

PRINT*,"  » 

PRINT*,"  INDICATE  WHICH  WATER  QUALITY  ATTRIBUTE(S)  ARE- 

PRINT*, »  TO  BE  ANALYZED." 

DO  40  1 3  1  ,  4 

40  PRINT*,"  " ,  CNAME ( I ) , "  •, !♦*,«)  " , CNAME ( 

IF ( NTWQC . LE . 8  )  GO  TO  45 
DO  42  lag, NTWQC 

42  PRINT*,"  ",I,">  " , CNAME( I ) 

PRINT*,"  RESPOND  WITH  THE  TOTAL  NUMBER  OF  ATTRIBUTES  " 

PRINT*,"  FOLLOWED  BY  THE  APPROPRIATE  INDEX  NUMBERS" 

45  PRINT*,"  ", 

READ ( 3 , * )  NATTS, ( I ATTS ( I ) , I s 1 , N ATTS ) 

IF ( EOF ( 3 ) )  999,46 

46  IFCNATTS.GT. NTWQC)  GO  TO  45 
PRINT*,"  " 

C - >  SPECIFY  EXISTING  WATER  QUALITY  STANDARDS 

PRINT* , " - 

PRINT*,"  INPUT  LOCAL  AMBIENT  WATER  QUALITY  STANDARDS  " 

PRINT*, " - 

PRINT*,"  " 

PRINT*,"  " 

ITPxl 

52  PRINT*,"  " 

PRINT*,"  TIME  PERIOD  NO.  ",ITP 

DO  59  ICsl, NATTS 

IWQCsIATTS(IC) 

PRINT*, »  " 

PRINT*, »  " , CNAM£( IWQC ) 

PRINT*,"  UPPER  LEVEL  STANDARD  -, 

READO.*)  UWQS(IWQC.ITP) 

IF ( EOF ( 3 ) ) 56 , 56 


1*4) 
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50  PRINT*,"  LOWER  LEVEL  STANDARD  ", 

R  EA  0  (  3  •  *  )  L  WQS  (  I  WQC  ,  I  TP  ) 

IF ( EOF ( 3  )  )  59,59 

59  CONTINUE 

IF ( ITP . £0. 1 )  GO  TO  60 
IT?*ITP*1 

IF ( ITP . GT . NTP )  GO  TO  69 
50  TO  52 

60  PRINT*,"  " 

PRINT*,"  STANDARDS  CONSTANT  OVER  TIME  (t  OR  N)», 

READ (3,990)  IANS 
IF  ( EOF ( 3 ) )  60,6) 

61  IFUANS.  EQ.  "N")  GO  TO  52 
DO  68  ITP*2 , NTP 

DO  66  IC*1, NATTS 
IWqC*IATTS(IC) 

UWQS(IWQC,ITP)=UWQS(IWQC, 1 ) 

66  LWQS(IWQC,ITP ) sLWQS ( IWQC , 1 ) 

6d  CONTINUE 

c - )  read  DATA  FROM  'TAPS33' 

69  DO  100  ITP* 1 , NTP 

DO  99  ICs 1 , NATTS 
IWQC=IATTS(IC; 

CALL  INWQMS  C IWQC , ITP , X , WQ , N PP ) 

DO  90  IPPs 1 , NPP 

IF ( IPP , EQ . 1 )  GO  TO  90 

IF  C  WQ  ( I P  P )  .  LT . UWQS ( IWQC ,ITP).AND.WQ(IPP).OT. L WQS ( IWQC, ITP ) ) 

*  GO  TO  90 

XWQSV(IC,ITP)*XWQSV(IC,ITP)*<X<IPP, 1 )-X(IPP-1 ,1)) 

90  CONTINUE 

99  CONTINUE 

100  CONTINUE 
PRINT*,"  " 

PRINT*,"  * 

PRINT*,"  " 

PRINT*,  * . . . . 

PRINT*,"  REPORT  ON  WATER  QUALITY  VIOLATIONS" 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  1)  TOTAL  RIVER  MILES  IN  VIOLATION  OF  STANDARDS." 

PRINT*,"  * 

PRINT*,"  TIME  PERIOD- 

PRINT*,"  ATTRIBUTE  1  2  3  «  5  6" 

PRINT*,"  -  -  -  -  -  -  - 

DO  300  IC=1, NATTS 
IWQC*IATTS(IC) 

NEND*6 

IFCNTP.LT. 6)NEND*NTP 

PRINT ( 4, 930)  CNAME(IWQC) , (XWQSV(IC,J),J»1,NEND) 

300  CONTINUE 

IF (NTP . LE . 6  )  GO  TO  399 


PRINT*,"  " 

CALL  PAUSE 

PRINT*,"  TIME  PERIOD" 

PRINT*,"  ATTRIBUTE  7  8  9  10  11  12- 

PRINT*,"  -  -  -  -  -  -  - « 

DO  310  CC : 1 , NATTS 
IWQCsIATTS(IC) 


PRINTC 1,930)  CNAME ( IWQC ) , ( XWQSV ( IC , J ) , J*6 , NTP ) 
310  CONTINUE 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  " 

CALL  PAUSE 
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c - >  ask  for  graphical  output 

399  PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  00  YOU  WANT  A  GRAPHICAL  OUTPUT  FOR  ANT  ATTRIBUTE" 

400  PRINT*."  PROFILES  (ANS:  TES  OR  NO)", 

R£AD<3.940)  IANS 

IF< IANS. NE. "T" . AND. IANS . NE. "N" )  GO  TO  ROO 
IF ( IANS . EQ . "N * )  GO  TO  899 
CALL  USTART 
RIO  CALL  UBELL 
PRINT*,"  « 

PRINT*,"  ATTRIBUTE  NUMBER  ", 

READ ( 3 , * )  IWQC 
IF ( EOF  C  3 ) )  R10.R11 
R11  IF(IWQC.GT.NTWQC)  GO  TO  RIO 
R20  PRINT*,"  ■ 

PRINT*,"  TIME  PERIOD  ", 

READ ( 3 , * )  ITP 
IF ( EOF ( 3 ) )  420,421 
421  IF(ITP.GT.NTP)  GO  TO  420 

CALL  INHQMS( IWOC , ITP , X , WQ.NPP ) 

DO  430  IC= 1 , NATTS 
IF(IATTS(IC).NE.IWQC)  GO  TO  430 
USxJWQSdWQC.ITP) 

LSaLWQS ( IWQC , ITP ) 

GO  TO  431 

430  CONTINUE 

431  IF ( US . LT. 9999 . )  GO  TO  460 
WQ(NPP*1 )*LS 

WQ  (  N  P  P  +  2  )  s  LS 
GO  TO  470 

460  WQ(NPP*1 )rUS 
WQ(NPP*2)=US 
470  DO  479  Ixl, NPP 

479  XX(I)=X(I,1) 

XX ( NPP+ 1 )*0.0 
XX(NPP>2) -XX(NPP) 

PTS ( 1 JxNPP 

PT3 ( 2 ) s2 

PRINT(4,910)  CNAMEdWQC)  ,ITP 
CALL  G0PL0T1 (XX, WQ, PTS, IWQC) 

PRINT*,"  » 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  PLOT  ".CNAMEdWQC),"  FOR  ANOTHER  TIME  PERIOD  ", 

480  READ(3.940)  IANS 
IF ( EOF ( 3 ) )  480,481 

481  IFdANS.NE."Y". AND. IANS. NE.*N")  GO  TO  480 
IF ( IANS . EQ . "X"  )  GO  TO  420 

PRINT*,"  PLOT  ANOTHER  WATER  QUALITY  ATTRIBUTE  ", 

490  READ(3,940)  IANS 
IF(E0F(3))  490,491 

491  IFdANS.  NE  .  "Y"  .  AND.  IANS  .NE.  "N")  GO  TO  490 
IF(IANS.EQ. "Y")  GO  TO  410 

CALL  UEND 
940  FORMAT ( 1  A  1  ) 

910  F0RMAT(25X,  1A5,":  TIME  PERIOD  NO.  " , 112) 

930  FO RMAT ( 4X ,  1 A5 , 9 X , 6 ( 1 f 5 . 1 , 2X  ) ) 

899  PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  " 

999  PRINT*,"  " 


PRINT*," . . . . . . . . 

PRINT*,"  THIS  CONCLUDES  'WQRTV ' .  YOU  MAY  EXECUTE  MORE" 

PRINT*,"  RTV  ROUTINES  NOW,  PEGIN  A  MITIGATION" 

PRINT*,"  LOOP  OR  SIGNOFF." 


PRINT*,  "iiiiiiitiitilltlltlicrttiiitufitisiitiiiiiiiinitiii" 

STOP 

END 


88 


r 


SUBROUTINE  DATINMS 

i  COMMON /NAMES/ TITLE (6 ) .CNAME (20) 

COMMON/CNTRL/NH  ,  NTP, WWQC.NB,  .VI T.  (DIE,  IT(  1 5  )  , IE( 1 5  )  . 

‘  *  UCOOS(d),IJCWgC,NNCWOC.NTWaC.NPP 

j  COMMON /MO  A  T  A/ LR(  20)  ,  ilk  (  20  )  ,  BCODEi  15) 

I  SEAL  LR 

INTEGER  TI TLE , HCODE, CNAME 

I  INTEGER  DINDEX( 125) .CNTRLS(d) 

C - >  UNLOAD  DATA  FROM  MASS  STORAGE  FILE 

CALL  0PENMS(9,DIN0EX, 125,0) 

CALL  SEADMSC 9 , TITLE, 5 ,  1  ) 

CALL  READMS(9,CNTRLS,a,2) 

NRsCNTRLS (  1  ) 

1  NT?  =  C  NTRLS ( 2 ) 

NCWQC=CNTRLS( 3) 

NNCWOCsCNTRLSC A ) 

HCOO£=CNTRLS( 5) 

NB*CNTRLS(6 ) 

NXTsCNTfllSl  7) 

NIE«CNTRLS( 8) 

NWOCsS 

HKs 1 6*NNCWQC 

MSs  a»NNCWQC»NC'NQC 

NTWQCs8*NCUQC*NNCWQC 

N A ATTS  sNCWQC+NNCUQC 

IF  (NAATTS . LS . 0 )  GO  TO  570 

CALL  SEA0M5(9,Ct?AM£,NTWQC,9> 

570  CALL  REAOMS( 9. BCODE, N8, 5) 

CALL  READMS(9,IT,NIT,6> 

CALL  READMS(9,IE,NIE,7) 

CALL  READMS(9,LR,NR,8) 

CALL  READMS( 9 , DA , HR , 9 ) 

RETURN 

END 

SUBROUTINE  GOPLOT 1 ( X , T , PTS , IWQC ) 

COMMON /NAMES /T ITLE (  8 ) ,CNAME(20) 

INTEGER  TITLE, CNAME, OPTSU) 

DIMENSION  X(290),]((290),PTS(2) 

DATA  OPTS/"LO" , ”DS"/ 

CALL  UDIMEN ( 7 . ,5.25) 

CALL  USET( “EDGEAXES" > 

CALL  UPSETl "CHARACTER" ,"♦" ) 

CALL  USETt "XBOTH") 

CALL  USET("XBOTH") 

call  UPSET< "XLABEL" , "DISTANCE  DOWNSTREAM  (MILES);") 
CALL  UPSET( "TLABEL" , "CONC .  (MG/L);") 

CALL  UBELL 

CALL  UPL0T(X,T,2. , PTS, OPTS) 

CALL  UFLUSH 
CALL  UPAUSE 
CALL  UERA SS 
RETURN 
END 

SUBROUTINE  PAUSE 
PRINT* , "CONTINUE" , 

READC  3 , * )  DUM 
IF ( EOF ( 3  )  )  10,  10 
10  CONTINUE 
RETURN 
END 
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SUBROUTINE  INWQMS ( I ATT , ITP , X , ATT ,  N PP ) 
DIMENSION  WQ(200),X(200,3) 

INTEGER  0INDEX(277) 

IF(IOPEN.GT.O)  GO  TO  10 
CALL  OPSNMSC 33, OINDEX, 277,0) 

10  IREC1*23,(ITP-1 )  ♦  1 

CALL  READMS ( 33 . NPP , 1 , IREC 1  ) 
IREC2*23#(ITP-1 )*2 
CALL  READMS(33,X,600,IREC2) 
IH£C3*23*(ITP-1 )»3+IATT 
CALL  READMS(33.ATT,NPP,IREC3) 

IOPENs  1 

RETURN 

END 

/ 


type , aairt v/g 


90 


f 

t 

I 

i 


20 

21 


50 

90 

100 


150 


160 


PROGRAM  3  IMA  IN ( I N PUT , OUT  PUT , TA P E 1  =  I N PUT , TA PE2 aOUT PUT ,  t  A  P  E  9 , 

♦  T  APE  3  3 ) 

COMMON/NAME3/TITL£(8) ,CNAME(20) 

C0MM0N/CNTHL/NH.NTP,NWQC,N8 , MIT . NIE . IT ( 1 5 ) . IE ( 1 5 ) . 

♦  MCODEl 3) , NCWQC , NNCWGC , MTWQC . NPP 
COMMON/MDATA/BCODEC 15) 

DIMENSION  X(200, 3) , 30D5 ( 200 > , S I ( 200 ,  12 ) , XSI ( 1 0 ,  1 2 ) , SIT ( 7 ) 
INTEGER  TITLS,CNAM£,SNAME(7,2) 

DATA  X3 1/1 20 *0.0/, SIT/-. 5, 1.5, 2. 5, 3. 5, 5, 5. 5, 6. 5/ 

DATA  SNAME/ " PUREST  W" , "CLEAN  WA","MILO  POL" , " POLLUT ED " , 

♦  "HEAVILY  ","RAW  SENA " , "SEPTIC  C","ATER  ", 

♦  "TER  " , "LUTION  "."POLLUTED", 

♦  "GE  " , "ONDITION"/ 


CALL  DATINMS 

DO  100  I  TP  *  1 , NTP 

CALL  INWQMS(2,ITP,X,S0D5,NPP) 

DO  90  IPP=  1 , M P P 
IF(B0D5(ITP).GT.50. )  GO  TO  20 

SI( IPP , ITP)s( 1 . 0747"B0D5( IPP)-0 . 4729) /(O. 90 408+0. 2 1 8*BOD5 ( I PP ) ) 
GO  TO  21 

SI(IPP,ITP)=(O.Ol89"BOD5(IPP)+7.938)/(1 . 8 82 -0 . 002 1 "B0D5 ( I P P ) ) 
IFCIPP-EQ.  I  )  GO  TO  90 
DO  50  1=1,7 

IF  (SHIPP,  ITP)  .GT.SIT(I)  )  GO  TO  50 
XSI(I.ITP)=XSI(I,ITP)*X(IPP,1)-X(IPP-1,1) 

GO  TO  90 
CONTINUE 
CONTINUE 
CONTINUE 
PRINT" , " 

PRINT", " 

PRINT", " 

PRINT", " 

PRINT" , " 

PRINT", " 

PRINT" , " 

PRINT", " 

PRINT", " 

N  END: 6 
IF ( NT? . 

PRINT", 

PRINT", 

PRINT" , 


SAPROSIC  INDEX  ANALXSIS" 
FOR" 

■,(TITLE(I),I=1,4) 

",(TITLE(I),I=5,8) 


333533333333333333333=33333X33333X33333 


333333333333335" 


LT. 6)NEND=NTP 
"  WATER  QUALITY 

"  DESIGNATION 


RIVER  MILES  IN  TIME  PERIOD" 

2  3  4  5  6 


PRINT","  " 

DO  150  1=1,7 

PRINT <2, 901)  (SNAME(I, J) , J= 1 , 2 ) , ( XSI ( I , ITP ) , I TP= 1 , NEND ) 
CONTINUE 

PRINT","  -  - - - - 

ii  _ _ " 


PRINT","  " 

IF ( NTP . LE . 6 )  GO  TO  200 
NENDx  1  2 

IF ( NTP . LT . 12)NENDsNTP 
NEND=6 

PRINT","  WATER  QUALITY  RIVER  MILES  IN  TIME  PERIOD" 

PRINT","  DESIGNATION  7  8  9  10  11  12 

PRINT","  -  - - - - - ", 

PRINT","  " 

DO  160  1=1,7 

PR INT( 2 , 90  1  )  (3NAME(I,J),J=1,2),(XSI(I,ITP)  ITP=1,NEND) 

CONTINUE 

PRINT","  -  - - - - - ", 
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200  CONTINUE 
PRINT*,"  » 

PRINT*,”  DO  YOU  WANT  FURTHER  QUANTIFICATION  OF  THIS" 

201  PRINT*,"  (ANS:  YES  OR  NO)", 

R£AD( 1 ,902)  IANS 

IF ( EOF ( 1 ) )  201,202 

202  IF ( IANS. N£ . ”N" . AND. IANS . NE. "Y")  00  TO  201 
IF( IANS . SQ. "N" )  00  TO  300 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  INPUT  TIME  PERIOD  OF  INTEREST  ", 

210  R  E AD ( 1  , • )  I  TP 

I? ( EOF ( 1 ) )  210,211 

211  IF(ITP.GT.NTP)  00  TO  210 
PRINT*,"  " 

PRINT*,"  THIS  SECTION  ISN'T  OPERATIONAL  YET,  BUT  THE 

PRINT*,"  OUTPUT  WILL  BE  LOCATIONS  OF  20NES  IN  EACH" 

PRINT*,"  WATER  QUALITY  DESIGNATION  FOR  THE  SPECIFIED 

PRINT*,"  TIME  PERIOD." 

300  CONTINUE 

901  F0RMAT(4X, 2A8, 3X,6( IX, 1F5. 1 ) ) 

902  FORMAT ( 1  A  1  ) 

STOP 

END 

SUBROUTINE  DATINMS 

C333333333333323333333322333S33333333S33XX23223X33332X3S333S3X 

COMMON/ NAMES /TITLE! 8) ,CNAM£(20) 

COMMON /CNTRL /NR, NTP, NWQC, N8, NIT, NIE, IT ( 15 ) , IE( 1 5) , 

♦  MC0DE(8),NCWQC,NNCWQC,NTWQC,NPP 
COMMON /MD ATA /LR ( 20) ,DA(20) , BC  ODE ( 15) 

REAL  LR 

INTEGER  TITLE. HCODE.CNAME 
INTEGER  DINOEXC 1 25  ) , CNTRLS < 8 > 

c - >  UNLOAD  DATA  FROM  MASS  STORAGE  FILE 

CALL  0PENMS(9.DINDEX, 125,0) 

CALL  READMS (9, TITLE, 8, 1 ) 

CALL  READMS(9,CNTRLS,a.2) 

NRaCNTRLSC 1 ) 

NTP  aCNTR  LS ( 2 ) 

NCWQCsCNTRLSC  3) 

NNCWQC=CNTRLS( 4 ) 

HC0DExCNTRLS(5) 

NB’CNTRLS ( 6 ) 

NITxCNTRLSI 7 ) 

HlEzCNTRLSi d ) 

NWQC  *  3 

NKs  1  6-»NNCWQC 
NS=8»NNCWQC*NCWQC 
NTWQCx8*NCWQC*NUCWQC 
NAATTSsNCWQC^NHCWQC 
IF(NAATTS.LE.O)  GO  TO  570 
CALL  READMS { 9 , CN AME , NTWQC , 4 ) 

570  CALL  READMS (9,BCODE,NB,5) 

CALL  READMSt 9 , IT , NIT, 6) 

CALL  READMS(9,IE,NIE,7) 

RETURN 

END 
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Cs3::::s338s3:3s3s33s339i3ss:ss:8s8sssi3f9x:stsa«asi3:2s:3ssi23xsst:3ss:s 

SUBROUTINE  INWQMS ( I ATT , IT P , X . ATT , NP P ) 
C:s::>n:s::s:::333:s333ssi(asssfs<:s3si«::iss3i3na(:isia:isi33itiitss3 

DIMENSION  WQ(200),X(200,3) 

INTEGER  QINDEX (277) 

IF ( ITP . GT . 1 )  GO  TO  10 

CALL  OPENMS(33.QINDEXt277,0) 

10  IREC1s23,(ITP-1 )*1 

CALL  R£ADMS(33»NPP, 1.IREC1) 

IREC2i23*(ITP-1 ) *2 

CALL  READMS(33,X,600,IR£C2) 

IREC3*23*(ITP-1 )♦ 3+1 ATT 
CALL  R£ADMS(33,ATT,NPP,IREC3) 

RETURN 

END 


PROGRAM  TUTSST  (  INPUT,  OUTPUT  ,  TAPE 3*IMPUT , TAPE4 sOUTPUT  , 

♦  TAPS'}, TAPE33) 

COMMON /NAM£3/TITLS(  3)  ,  CNAME (  20) 

COMMOM/CNTRL/NR, MTP.NWOC, SB, HIT, HIE, 

♦  MCODS ( 8 ) , NCWQC , NNCWQC , NTWQC , NPP 
C0MM0N/MD.4TA/LR (  20  )  , DA ( 2 0 ) , BCODE { ’5) 

DIMENSION  X (200, 3) , WO (200) , PTS( 2) , TU ( ’ 0 , 5 , 200 ) , MAXTU ( 1 0 , 5 > , 

♦  MEANTU  (  ’  0 , 5  )  ,  SPNAM£(  '  0, 2)  ,  T  N  A  M  £  (  10),XX(400),LC50(70,5), 

♦  ITDX( 10) , ISP( 70) ,C<  400) ,  TOTTUf  5,2  00)  ,  ISOC  7  0) 

REAL  LR , LC50 , MAXTU , MEANTU , MAXTTU ( 5 ) , MEANTTU ( 5 ) 

INTEGER  TITLE, CNAME , TNAME , SPNAME 

DATA  ( CN AMS ( I ) , I s 1 , 8 ) / " TEMP . " , "B0D5 " . “TSS" , "NH3" , "N02" , "N03 " , 


♦  "P04", "D.O. ”/ 

DATA  ( SPNAME ( 7 ,1) , 1  =  1 ,2)/” FATHEAD  MI" ,"NNOW  "/ 
DATA  (SPNAME ( 2, I ) , I = 7 , 2 )/"CARP  "/ 
DATA  (SPNAMEl  3  ,  I )  ,  1  =  ’  ,  2) /"BL'JEGILL  "/ 
DATA  (SPNAMEt 4, I) , 1= 7 , 2) /"CHANNEL  CA",”T  "/ 
DATA  (SPNAMEt 5 , I ) , 1= 7 , 2 ) / "L A RGEMOUTH " , "  BASS  "/ 
DATA  (SPNAME! 6, I ), 1=1 , 2) /"BROOK  TROU",,,T  "/ 
DATA  (SPNAME ( 7 , I ), 1=7 ,2 ) /"RAINBOW  TR","OUT  "/ 
DATA  (SPNAME(S.I) .1=7 ,2)/"C0H0  3ALM0","N  "/ 


DATA  ISO/S*1 ,5*-'/ 

DATA  TOTTU/’000*0.0/ 

PRINT*,"  » 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,  +  +  +  +  + 

PRINT*, "♦  THIS  RTV  ROUTINE  TESTS  ENVIRONMENTAL  TOXICITY  ♦” 

PRINT*, *♦=♦=♦*♦**=♦=**♦=******♦*♦*♦=♦*♦*♦=♦***♦=♦=♦=♦****=♦*♦»*" 
PRINT*,"  " 

C - >  READ  DATA  FROM  'TAPE9' 

CALL  DATIMMS 

C - )  DESIGNATE  ATTRIBUTES  OF  INTEREST 

PRINT*  ,  " - - - " 

PRINT*,"  " , (TITLE! I ) , 1= 7 ,8) 

*>>*■;  ^  *  ■■ _ _ _ _ _ _ _ „ _ _ _ _ _ _ _ _ _ _  _ _ _ _ „ 

PRINT*, "  " 

PRINT*,"  INDICATE  WHICH  WATER  QUALITY  ATTRIBUTE(S)  ARE" 

PRINT*,"  TO  BE  ANALYZED  FOR  THEIR  TOXIC  EFFECTS," 

DO  1 0  1=7,4 

70  PRINT*,"  ",I,")  " , CNAME ( I > , "  ",I*4,")  ",CNAME(I*4) 

IF ( NTWQC . LE -  8 )  GO  TO  15 
DO  12  I =9, NTWQC 

12  PRINT*,"  ",!,*)  " , CNAME ( I ) 

PRINT*,"  RESPOND  WITH  THE  TOTAL  NUMBER  OF  TOXICANTS  " 

PRINT*,"  FOLLOWED  BY  THE  APPROPRIATE  INDEX  NUMBERS" 

15  PRINT*,"  ”, 

READ( 3 , » )  NTOX, ( ITOX( I ) , 1=1 , NTOX) 

IF ( EOF (  3  )  )  999  ,  '6 

16  IF ( NTOX. GT. NTWQC)  GO  TO  15 
DO  ’3  1=1, NTOX 

IF( ITOX( I ) . LE. NTWQC )  GO  TO  '8 
PRINT*,"  BAD  INPUT;  TRY  AGAIN!" 

GO  TO  15 

18  TNAME ( I )=CNAME( ITOX(  I  )  ) 

PRINT*.”  " 

20  PRINT*,"  " 
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22 


25 


26 


C 

30 


32 

33 
36 

69 

TO 

71 


90 


92 

95 

99 


PRINT*  ,  "  SPECIF!  TARGET  SPECIES:" 

e  s  i  .9  r  • , "  " 

PRINT*,"  DESIGNATE  STREAM  T  If  P  S  <W=W  ARM  WATER, C*C0L  0  WATER)  ", 

READ! 3. 9  NO  )  I3TRM 
I F  (  EOF  13))  2’ ,22 

IF  (  ISTRM  .  NE  .  “C"  .  A.9D.  ISTRM.  NE  .  “W"  )  GO  TO  2' 

SPECIFICATION  OF  TARGET  SPECIES 


PRINT*,"  " 

PRINT*,"  REPRESENTATIVE  SPECIES  LIST:” 

DO  25  1=1,6 

IF( ISTRM . EQ . "W" . AND . ISO( I ) . LT . 0 .OR . 

♦  ISTRM. EQ . "C" . AND . ISQ( I ) .GT. 0)  GO  TO  25 
PRINT*,"  ",I,")  ", (SPMAMEd , J) , J* 1 ,2) 

CONTINUE 


PRINT* , " 
PRINT*, ” 
P  R  1 9  T  *  ,  " 


RESPOND  WITH  NUMBER  OF  TARGET  SPECIES  DESIRED  AND” 
WITH  THEIR  APPROPRIATE  INDEX  UUHBSR ( 3 ) . n 

II 


READ  13.*)  NS?,(IS?(I),I  =  ’,NSP) 

IF ( EOF  <  3 ) )  26,30 

SPEC  I F I  LC  5  0 ' S  FOR  TARGET  SPECIES 
PRINT*,"  * 

PRINT*,"  INPUT  THE  96  HOUR  LC50'S  FOR  THE  FOLLOWING  SPECIES" 

PRINT*,"  AND  POTENT  ION AL  TOXICANTS:" 

DO  35  IS=1,NSP 
PRINT*,"  " 

PRINT* , "  ", (SPNAHEt ISP ( IS ) , J 1  J=1 , 2) 

DO  33  I T  = ’ , NTOX 

PRINT*,"  " , TN  AME ( I T ) , 

READ(3.*>  LC50(IT,IS) 

I r ( EOF ( 3 ) )  32,33 
CONTINUE 
CONTINUE 
CALL  USTART 


PRINT*, "  " 

PRINT*,"  SPEC  I F I  TIME  PERIOD  OF  INTEREST  ", 

READ ( 3 ,  * )  I  TP 

Ir (EOF ( 3  )  )  70,71 

IF (ITP.GT.NTP)  GO  TO  70 

PRINT*,"  " 

DO  99  IS= 1 , NSP 
DO  95  IT=1,NT0X 
3'JMTU  =  0  . 0 


SUMTTU  =  0 . 0 
MA  XTU ( IT ,  IS) =0 . 

ICsITOX(IT) 

CALL  INWqMSI IC , ITP , X , WQ, NPP ) 

DO  90  IPPsl.HPP 

TU(IT,IS, IPP)=WQ(IPP)/LC50(IT,IS) 

IF  (MAXTUdT,  IS)  .  LT.TUdT,  IS  ,  IPP)  )MAXTU(IT,  IS  >  =  TU(  IT,  IS,  IPP) 
TOTTU( IS , I P  P ) s  TOTT  U ( IS, IPP ) +TU( IT, IS, IPP ) 

SUMTU  =  SUHTU*TU ( IT , IS, IPP ) 

MEANTU( IT, IS ) =  SUMT  U /N  PP 

MAXTTU(IS)=0.0 

DO  92  I PP= 1 , NPP 

IF(MAXTTUdS)  .LT.  TOTTU  (  IS.  IPP  )  )  MAXTTUCIS  )  =  TOTTU  I  IS  ,  IPP  ) 
3UMTTU*SUMTTU*T0TTUdS,  IPP) 

MEAIITTUdS)  rSUMTTU/NPP 


CONTINUE 
CONTINUE 
PRINT*," 
PRINT*," 
PRINT*," 
P  R  I  N  r  *  .  " 
PRINT*," 
PR  INI*, ■ 


REPORT  ON  TOXICITT 


IMPACTS 


IN  TIME  PERIOD 


,  ITP 


t 
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AD-AU1  9H? 


UNCLASSIFIED 


ILLINOIS  UNIV  AT  URBANA  DEPT  OF  CIVIL  ENGINEERING  F/S  13/2 

QUANTITATIVE  ASSESSMENT  OF  ENVIRONMENTAL  IMPACTS  IN  THE  A8UATIC— ETC(U> 
JAN  82  R  RIG6INS*  E  HERRICKS.  M  J  SALE  DACA0S-78-R-OO6 

CERL-TR-N-119  ML 


PR I  NT*,"  " 

print*.  *  " 

PRINT*, "  MAXIMUM  AND  MEAN  (IN  PARENTHESES )  TOXICITY  UNITS" 

PRINT*."  " 

PRINT*,"  TOXICANT" 

N  EN  D  :  5 

IF(NTOX.LT.5)  NEND=NTOX 

PRINK  3,91  1 )  (TNAME(IT) , IT: 1 , NEND ) 

DO  300  IS i 1 , NS P 

PRINK  4, 930)  ( S  PN AME (ISP(IS).J) ,  J:1 , 2 ) , MAXTTU ( IS  )  , 

♦  (MAXTUCI ,IS) , Is  1 ,NEND) 

PRINK  4, 931  )  MEANTU ( IS ) , ( MEAMTU ( I , IS ) ,1:1, NEND ) 

300  CONTINUE 

PRINT*. ”  " 

PRINT*,"  " 

PRINT*,"  " 

C - >  ASK  FOR  GRAPHICAL  OUTPUT 

3  ?  5  PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  » 

PRINT*. "  DO  YOU  WANT  A  GRAPHICAL  OUTPUT  FOR  TOXICITY  UNITS" 

400  PRINT*,"  VS.  LOCATION  DOWNSTREAM  (ANS:  YES  OR  NO)" 

R£AD(3.940>  IANS 

I?(IANS.NE."Y". AND. IANS. N£."N")  GO  TO  400 
IF(IANS.EQ."N")  GO  TO  399 

410  CALL  UBELL 
PRINT*,"  " 

PRINT*, «  INPUT  TOXICANT  NUMBER  ", 

REAO ( 3 . * )  IT 
IF(£0F(3))  412,411 

411  IF (IT.LE.NTOX)  GO  TO  420 

412  DO  415  I T* 1 , NTOX 

415  PRINT*,"  " , IT , " )  ",TNAME(IT) 

PRINT*."  ", 

READ ( 3 • *  )  IT 
IF (EOF  C  3 )  )  412,411 

420  PRINT*,"  " 

PRINT*,"  INPUT  TARGET  SPECIES  INDEX  ", 

READ ( 3 , * )  IS 
IF (EOF ( 3 ) )  422,421 

421  IF(IS.LE.NSP)  GO  TO  449 

422  DO  425  IS: 1 , NSP 

425  PRINT*, "  ",  IS, ")  ", (SPNAME(ISP( IS) . J) , J:1 ,2) 

PRINT*,"  ", 

READ ( 3 , *  )  IS 

IF ( EOF (  3  )  )  422,426 

426  IF(IS.GT.NSP)  GO  TO  420 

449  DO  450  IPP*1 , NPP 
XX(IPP)=X(IPP, 1) 

C(IPP)»TU(IT,IS,IPP) 

XX(IPP*NPP):X(IPP,  1  ) 

450  C(IPP*NPP3*T0TTUCIS,IPP) 

PTS ( 1 ) *  M  P  P 

PTS( 2 ) *NPP 
PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  " 

PR  I  NT ( 4, 910)  (SPNAME(ISP(IS) . J) , J>1 ,2) ,TNAME( IT) , ITP 

PRINT*,"  (KTOTAL  T.U.'S;  *:T.U.'S  FROM  SPECIFIED  ATT.)" 

CALL  GOPLOT1 (XX.C.PTS) 

PRINT*,"  " 

PRINT*,"  * 

PRINT*,"  " 

PRINT*,"  PLOT  TOXICITY  IMPACTS  FROM  " , TNAME( IT) 

PRINT*,"  FOR  ANOTHER  TARGET  SPECIES  ", 
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READ(3.940>  IANS 

I F  !  £  0  F  (  3  >  )  4d0.4d'  , 

If ( IANS . NE . "i" . AND . IANS . NS . "M" )  00  TO  4d0 
IfdANS.SO.  "if")  00  TO  420 

PRINT*, "  PLOT  FOR  ANOTHER  TOXICANT  ”, 

REAQ( 3.940)  IANS 
If  C  EOF  C  3 )  )  490,491 

IF( IANS. NE . " i " . AND . IAN3.NE . "N”  )  00  TO  490 
If  (  IANS  .  EQ  .  "I"  )  00  TO  4’0 

PRINT*,"  DO  YOU  WISH  TO  CONTINUE  ANALYSES  FOR  ANOTHER  TIME" 

PRINT*,"  PERIOD  CANS:  i  OR  N ) " , 

READ! 3 ■ 940 )  IANS 
If ( EOF ( 3 ) )  500,50' 

If  (IANS.  NE."jf".  AND.  IANS.  NE."N")  GO  TO  500 
PRINT*. '•  " 

IFCIANS . SO. "Y " )  GO  TO  69 
,auL  UEND 
FORMAT! ’12) 

FORMAT! ’5X, 2A10 , "  ;  ",’A5,"  ;  TIME  PERIOD  «,1I2) 

FORMAT!  3 X, "TARGET  S PEC  I ES " , 3X . "TOTAL" , 3X,5(2X,’A5.2X),/, 

*  'X,  20!  "-"  )  ,  2X,  6!” - ",1X)) 

FORMAT! ’X,2A’0,'X,6C1X.1F7.3>1X)) 

FORMAT!  2 2X  ,  5  (  1  X  .  "  (  "  ,  1 F 6  .  3  ,  "  )  "  )  ) 

FORMAT!  U!) 

PRINT*,  "  ■" 

PRINT*. "  " 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  " 

PRINT*,"  THIS  CONCLUDES  '  T  U  R  T  V  ' .  YOU  MAY  EXECUTE  MORE" 

PRINT*,"  RTF  ROUTINES  NOW,  BEGIN  A  MITIGATION" 

PRINT*,"  LOOP  OR  SIGNOFF." 

PRINT*,"!:::::: . . . . 

STOP 

END 

SUBROUTINE  DATINMS 

COMMON /NAMEG/TITLE! 3) ,CNAME(20) 

C0MMJN/CN7RL/NR,NTP,NWqC,N8,NIT, NIE, 

♦  MCODE(d) .NCWQC.NNCWQC.NTWQC.NPP 
COMMON/MDATA/LR! 20 ) , DA (20 ) , BCODE! 15) 

REAL  LR 

INTEGER  TITLE, HCOOE,CNAME 
INTEGER  DINDEX! 125 ) .CNTRLSta) 

-->  UNLOAD  DATA  FROM  MASS  STORAGE  FILE 
CALL  0PENM319, DINDEX, 125,0) 

CALL  READMS(9, TITLE, 8, 1 ) 

CALL  READMS(9,CNTRLS,8,2) 

NRsCNTRLS! 1 ) 

NTP=CNTRLS(2) 

NCWQCsCNTRLS! 3) 

NNCWQCiCNTRLS! 4) 

HC0DE*CNTRLS(5) 

MB  =  CHTRLS ( 6 ) 

N IT  =  CNTRL3 ( 7 ) 

NI£=CMTRLS(3) 

NWQC : 3 

NKs 16*NNCWgC 

NS »  d*NNC WQC  »NC WQC 

NT  WQC >  d*NC WQC  »MMCWQC 

NAATTS»NCWgC*NMCWQC 

If(NAATTS.LE.O)  GO  TO  570 

CALL  READMS! 9.CNAME, NTWQC . 4 ) 

CALL  READMS!), BCODE, MB, 5) 

CALL  READMS(9. IT, NIT, 6) 

CALL  R E ADMS ( 9 ,  IE, Nit. 7) 
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CALL  READMS(9,LR,NR,d) 

CALL  READMS(9,DA,NR,9) 

RETURN 

END 

SUBROUTINE  GOPLOT 1 ( X . T , PTS ) 

COMMON /N  AMES /TI TLE ( 6 ) ,CNAME(20) 

INTEGER  TITLE, CNAME , OPTS ( 2 ) 

DIMENSION  X ( 400  )  . X ( 400 ) .  PTS ( 2 ) 

DATA  OPTS/"L*" , "LT"/ 

CALL  UDIMBN (7.5,5.20) 

CALL  US£T< "EDGEAXES") 

CALL  UPSETt "CHARACTER" ,"." ) 

CALL  USET ( " X80TH "  ) 

CALL  USET ( " X30TH " ) 

CALL  UPSET("XLABEL", "DISTANCE  DOWNSTREAM  (MILES);") 
CALL  UPS£T( "ILABEL" , "TOX.  UNITS;") 

CALL  UBELL 

CALL  UPL0T(X,Y,2. , PTS, OPTS) 

CALL  UFLUSH 
CALL  UPAUSE 
CALL  UERASE 
RETURN 
END 

SUBROUTINE  PAUSE 
PRINT*, "CONTINUE", 

READ ( 3 • *  >  OUM 
IF ( EOF ( 3  ) ) 1 0 ,  10 

RETURN 

END 

SUBROUTINE  INWQMS (I ATT, ITP, X, ATT.NPP) 

DIMENSION  WQ( 200 ) , X ( 200 , 3) 

INTEGER  QINDEX(277) 

IF (IOPSN.GT.O)  GO  TO  10 
CALL  0PENMS(33,QINDEX, 277,0) 

10  IREC1  s23*(ITP-1  )♦’. 

CALL  REAOMS( 33.NPP, 1.IREC1) 

IREC2s23*(ITP-1 )»2 

CALL  READMS(33,X,600,IREC2) 

IH£C3»23*(ITP-1 )*3*IATT 
CALL  READMS( 33, ATT.NPP, IREC3) 

IOPENs 1 

RETURN 

END 
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Ohio  River  46201 
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USA  ARRCON  61299 
ATTN:  0RC1S-RI-1 
ATTN:  DRSAR-IS 

DARCOH  -  Olr.,  lost..  6  Svcs. 

ATTN:  Facilities  Engineer 
ARRADCQH  07801 

Aberdeen  Proving  Ground  21005 
Amy  Hatls.  end  Nechenlcs  Res.  Ctr. 
Corpus  Chrtsti  Amy  Depot  78A19 
Horry  Olaaond  Laboratorlps  20783 
Ougway  Proving  Ground  84022 
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Sharp*  Amy  Oapot  953)1 
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Tobyhann*  Amy  Depot  18468 
Tooele  Amy  Depot  8407 A 
Hetorvl let  Arsenal  12189 
Turn  Proving  Ground  86364 
volte  Sends  Hlsslle  Rouge  88002 

OLA  ATTH:  0LA-V1  22314 

FOASCOH 

FOASCOH  Engineer,  ATTN:  AFEN-FE 
ATTN:  Facilities  Engineer 
Fort  Buchenen  00934 
Fort  8rtgg  28307 
Fort  Caapbel  1  42223 
Fort  Corson  8Qgi3 
Fort  Oovons  014)3 
Fort  Drua  13601 
Fort  Hood  76644 
Fort  lndlontoun  Gap  17003 
Fort  Inrtn  92311 
Fort  Soa  Houston  78234 
Fort  Louis  98433 
Fort  HcCoy  64666 
Fort  HcPhorson  30330 
Fort  Goorgo  G.  made  20 755 
Fort  Ord  93941 
Fort  Polk  7145g 
Fort  Richardson  99505 
Fort  Alloy  68442 
Prypldlo  of  Sen  Frenclsce  94129 
Fort  Shorldon  60031 
Fort  Steuort  3131) 

Fort  Volnurlght  99703 
Toncovvor  Iks.  98660 

HSC 

ATTN:  HSLO-F  ?8234 
ATTN:  Facilities  Engineer 
Ft tJESlMOivs  Amy  Medical  Cmttr  90240 
Ml  Ur  Retd  Amy  Medical  Center  20012 

INSCOM  -  Ch,  Inttl.  Olv. 

ATTN:  Facilities  Engineer 
Arlington  Hall  Station  (2)  22212 
Vint  HIM  Fam*  Station  221SO 


MM 

ATTN:  facilities  Engineer 
Caaeron  Station  22*14 
Fort  Lesley  J.  McNair  20319 
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